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Abstract 


A method is given for the numerical solution of the barotropic vorticity equation 
over a limited area of the earth’s surface. The lack of a natural boundary calls for an 
investigation of the appropriate boundary conditions. These are determined by a 
heuristic argument and are shown to be sufficient in a special case. Approximate 
conditions necessary to insure the mathematical stability of the difference equation 
are derived. The results of a series of four 24-hour forecasts computed from actual 
data at the soo mb level are presented, together with an interpretation and analysis. 
An attempt is made to determine the causes of the forecast errors. These are ascribed 
partly to the use of too large a space increment and partly to the effects of baroclinicity. 
The réle of the latter is investigated in some detail by means of a simple baroclinic model. 


I. Introduction 


Two years ago the Meteorological Research 
Group at the Institute for Advanced Study 
adopted the general plan of attacking the 
problem of numerical weather prediction by 
a step by step investigation of a series of models 
approximating more and more the real state 
of the atmosphere. In accordance with this 
plan the two-dimensional barotropic model 
was chosen as the first object of study. The 
first two publications? dealt with the numerical 
properties of the linearized barotropic equa- 
tions as a preparation for the numerical integra- 
tion of the non-linear equations. Such integra- 
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tions have now been performed and will be 
described in the present article. 

These integrations would not have been 
possible without the use of a high-speed large- 
capacity computing instrument. We should 
like, therefore, to express our warmest thanks 
to the U. S. Army Ordnance Department and 
the administration of the Ballistic Research 
Laboratories in Aberdeen, Maryland for 
having generously given us the use of their 
electronic computing machine (The Eniac 
[compare footnote $]). The request for the 
use of the Eniac was made on our behalf by 
the U.S. Weather Bureau and we should like 
to thank them also for their gratifying in- 
terest and support. 

The reasons for regarding the integration 
of the barotropic equations as an essential 
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first step of the general program are as follows: 
(1) An accumulation of evidence indicates 
first that the effects of baroclinicity do not 
manifest themselves in a steady, widespread 
conversion. of potential into kinetic energy, but 
rather in sporadic and violent local over- 
turnings accompanying what, for want of a 
better term, may be called baroclinic instability, 
and second, that when these effects are not 
predominant, the motion is quasi-barotropic. 
It is hoped therefore that the barotropic pre- 
dictions, by their agreements and disagreements 
with observation, will provide a basis for an 
a priori classification of the large-scale atmos- 
pheric motions. One has the suspicion that 
certain processes which have heretofore been 
classed as baroclinic will be found to have a 
barotropic explanation. (2) If the barotropic 
forecasts are found to be sufficiently accurate 
approximations to the upper flow, it is possible 
that they can be profitably incorporated into 
practical forecast procedure. (3) Just as the 
analysis of the linearized barotropic equations 
served as a pilot study for the integration of 
the non-linear, barotropic equations, so will 
these integrations supply the necessary back- 
ground for the treatment of the three- 
dimensional equations. 

The most easily integrated of the baro- 
tropic equations are the primitive Eulerian 
equations, in which the local time derivatives 
of the field variables are given explicitly in 
terms of their space derivatives. Although 
the virtual unobservability of the geostrophic 
deviation and the horizontal divergence 
renders the initial time derivatives in the 
Eulerian equations highly inaccurate, the in- 
dications are that, contrary to an earlier im- 
pression4, the error will occur only as a small 
amplitude gravitational oscillation about an 
essentially correct large-scale flow — provid- 
ing the space and time increments in the 
finite difference equations are chosen to 
satisfy the Courant-Friedrichs-Lewy condi- 
tion (1928) for the computational stability of 
the equation governing the motion of long 
gravity waves. However two main considera- 
tions led the writers to decide against at- 
tempting the integration of the Eulerian 
equations at this time. First, the computa- 
tional stability condition states that the ratio 


4 CHARNEY (1949, p. 12). 


J. G. CHARNEY, R. FJÖRTOFT, J. VON NEUMANN 


of the space to the time difference must ex- 
ceed the gravitational wave velocity, about 
300 m sec", and demands a time increment 
of some ıs minutes or less. This means that 
a twenty-four hour forecast would require 
nearly 100 time cycles for the integration, a 
formidable number for the machine that was 
available to the writers. Then there was also 
the difficulty that the characteristic property 
of large-scale non-divergent barotropic mo- 
tions, the conservation of absolute vertical 
vorticity, is obscured in the integration of 
the Eulerian equations — which apply as well 
to the divergent gravity motions — and the 
results thereby made extremely difficult to 
analyse and interpret. These difficulties hav- 
ing been regarded by the authors as serious, 
it was decided instead to base the computa- 
tion on the quasi-geostrophic, non-divergent 
vorticity equation, in which the sole depend- 
ent variable is the height z of a fixed iso- 
baric surface. 


a. SS 
Because of the elliptic character in En of the 


non-divergent vorticity equation, zand certain 
of its space derivatives must be specified as 
functions of time at the boundary of a limited 
forecast region. Since these quantities are 
known only initially, one has to fıx their values 
in a more or less arbitrary manner, and it 
becomes necessary to know how rapidly 
influences from the boundary propagate into 
the forecast region. This problem has been 
treated by CHARNEY (1949) who arrived at 
the conclusion that large-scale influences travel 
with a speed not radically different from that 
of the wind, and these influences consequently 
are provided for by integrating over an area 
not much larger than the forecast region. 
Although it is immaterial what values are 
prescribed at the boundary, it is nevertheless 
important that the mathematical form of the 
boundary conditions be known. Experience 
has shown that a violation of these conditions 
may lead to errors adjacent to the boundary 
which propagate into the interior with de- 
structive effect. Since one is dealing with 
boundaries at which the conditions are not 
naturally prescribed by the geometry of the 
motion, as for example at a wall, it is not imme- 
diately obvious what these conditions are. A 
heuristic argument will be advanced to show 
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that the following are probably correct: where 
fluid is entering the region enclosed by the 
boundary both z and the relative vorticity 
must be prescribed, but where fluid is leaving 
the region it is enough to prescribe z. 

The first part of the following discussion, 
the mathematical part, is devoted to the treat- 
ment of the boundary conditions, the method 
of solution of the finite-difference vorticity 
equation, and the computational stability 
criteria. The second part contains a description 
and analysis of the results of four twenty-four 
hour forecasts computed from actual data 
and a final section devoted to an account of a 
baroclinic model which is used to explain 
some of the barotropic forecast discrepancies. 


U. The Vorticity Equation 


We assume that the horizontal winds in the 
large-scale systems vary according to the law 


V (51, 52, p) = A(p) v (sr, $2, po) (1) 


where sr, s2 are orthogonal curvilinear distance 
co-ordinates on the sphere, p is the vertical 
pressure co-ordinate, and po the mean surface 
value of p. If a bar denotes a vertical pres- 
sure average, 1.¢., 


Po 
— I 3 
“= — x dp, (2) 
Po J Ê 


the integrated vorticity equation takes the 
approximate form (CHARNEY [1949, p. 383]) 


| 


—v-Wv (Ke +f); (3) 


where € is the mean relative vertical vorticity 
component, f is the coriolis parameter, and 


K = A2/(A)?. 
If A* is defined by 
Ak = A2/A 
and (3) is multiplied by A*/A we obtain 


Oe 


vH. 


where v* and £* are respectively the wind 


velocity and relative vorticity at the level p* 
defined by A (p*) = A*. Since approximately 


Serpe Bev veep- 
dev, (5) 


p* is the level of non-divergence. It is known 
(Cuarney [loc. cit.]) that p, the level at which 
A = A and v = v, is between 600 and 500 
mb and that A2/A is approximately 1.25. 
Hence we have 


en == ee 1.25 
|v | A A2 

and from the average variation of v with 

height we find that p* is approximately 100 


mb higher than p or between soo and 400 mb. 
We shall take this level to be 500 mb in the 
forecasts. 

Evaluating the vorticity of the geostrophic 
wind 


ye eng. (6) 
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where k is the unit vertical vector and z 
the height of the p* surface, we find 


CX & Fdvve= Am (7) 


and substituting this expression into (4) we 
obtain the quasi-geostrophic vorticity equation: 


an 02 
ISı OS2 


al 
dt 


on dz 
ae fm Js (n, 2), 


(8) 


(4,z) = 


where 7, is the absolute vorticity 


dix +f (0) 


Here A, is the surface spherical Laplacian 
operator, and J, is the Jacobian of n and z 
with respect to sr and s2. Equation (8) is 
taken tobe the basic equation governing the 
large-scale motions in a barotropic atmosphere. 
Its solution may be found iteratively by solving 
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dz . 3 
for == and extrapolating the motion forward 


in time, but for this purpose the boundary 
conditions must first be ascertained. 


UI. The Boundary Conditions 


‘Let R be a region of the earth bounded by 
the simple, closed, rectifiable curve C. If 
Az is known in R and z is given on C, z 

: N : Bu: 
may be obtained in R by solving a Poisson’s 
equation. Hence, if the boundary conditions 
are prescribed in such a way that z is always 
known on C and A, is known in R, the 
solution to (8) for the region R will be deter- 
mined. Suppose now that z is a given func- 
tion of time on C, so that the tangential 
derivative and therefore the normal velocity 
is fixed, and that we know its values initially 
in R. Since according to (4) the absolute 
vorticity is advected with the fluid, after a 
small time öf the distribution of A,z will be 
known everywhere in R except in the part, 
OR, which is penetrated by fluid from outside. 
If in addition A,z had been prescribed on 
that part of C where fluid was entering, we 
should also have been able to say exactly what 
vorticity had entered 6R, since the normal 
velocity on C is known. Hence we may 
assert that the motion is determined by the speci- 
fication of z everywhere on the boundary, and the 
vorticity on that part of the boundary at which 
fluid is entering the interior region. 

To supplement the foregoing heuristic argu- 
ment we now give a demonstration of the 
sufficiency of these boundary conditions in 
the case of a special two-dimensional non- 
divergent flow which, however, appears to 
exhibit the essential mathematical properties 
of the barotropic flow. It will also be shown 
that the specification of only the stream func- 
tion does not determine the motion. 

We consider the two-dimensional incom- 
pressible motion of an inviscid fluid on a 
circular cylinder, in which the force of gravity 
is directed radially inward so that no external 
forces act along its surface. Let y be directed 
along the axis of the cylinder and x be di- 
rected at right angles to y. Because of the 
incompressibility assumption we may intro- 
duce a stream function y and write for the 
vorticity equation 
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an equation corresponding closely to (8). 

Consider the ring shaped domain bounded 
by the two circles y=o and y =a and 
prescribe the boundary conditions y = x on 
each circle for all t > o. Also let y = x in the 
entire domain at tf = o. The motion consists 
initially of a uniform streaming parallel to the 
y-axis. It may easily be verified that there are an 
infinity of solutions of (10) satisfying the initial 
and boundary conditions of the form 


p= x + [BO —B(—a)] À — [8-8 0—y)], 
(11) 


where ß is limited solely by the requirement 
(12) 


It is clear, therefore, that a knowledge of y in 
the boundary does not determine the motion. 
But the motion is determined if Ay is also 
specified on the boundary where fluid is en- 
tering, Le. at y= 0: Let (Ay, = Fi 
then from (11) and (12) 


B (t) = ofort So. 


(Ay),=0 = #8 (0) = QUE 
and 


B(t) =f f Fly) dy dx + At = G(t) + At, 


where A is an arbitrary constant. Substitu- 
tion of this expression into (11) gives 


y=x +[G()—G—a)]*¥—[G)—G¢—a)], 


and y is completely determined. 


IV. The Solution of the Vorticity Equation 


The following method of solving the finite 
difference vorticity equation is well-adapted 
to a variety of high-speed computing machines, 
although it was chosen specifically for use on 
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the Eniac.5 It is not, however, recommended 
for hand computation. 

The spherical earth is first mapped con- 
formally onto a plane. If m is the magnifica- 
tion factor, the Laplacian and Jacobian opera- 
tors transform as follows: 


A, = m4; J, = mJ 


where A and J are the Laplacian and Jacobian 
on the plane. Thus (8) is transformed into 


2 (49 =J {2 | a 


n=hAz+f 
with 


Mae 
i 


The mapping used for the numerical inte- 
grations was the stereographic projection of 
the earth’s surface onto a plane tangent at the 
north pole. In this case we have the following 
relation between the geographical latitude y 
and the distance r from the pole on the map: 


(14) 


cos 9 


I+ sing’ 


where the radius of the equator on the map 
is chosen as the unit of distance. We obtain 


m= — 2 rer le 

= dp en 
f =2Qsing =22 —; (Es) 
Be um (4 Cask Oe 


where 2 is the earth’s angular speed of rotation. 
With the notation 


a AZ 
5 Electronic Numerical Integrator and Computer, 


Ballistic Research Laboratories, Aberdeen Proving 
Ground, Maryland. 
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the system (13) is replaced by 


OL =J (7, 2) ; (16) 


Since it is immaterial what values are 
assigned to z and Az on the boundary, as 
long as it is sufficiently far removed from the 
forecast region, we may prescribe the condi- 
tions 


Oz 

es 0, 

dt 

3 (17) 
noe ofor =, —0, 


where 2, is the tangential derivative of z 
taken in the direction that has the interior 
of the region on the left. 

For simplicity a rectangular area with sides 
L, and L, is chosen, and a rectangular grid of 


F F j Ee 
points is defined by the co-ordinates x = — i, 
) 


ise 
— Io 
y N 


with boundary linsi=o,i=pandj=o, 


ic) 


Ly 
j = 4. The grid intervals = and rn are taken 


to be equal to the common value A s. 


The quantities h and f are independent of t 
and may be determined once and for all from 


(15) and 


(ey Yo)? 


~ 
» 
I 


where x,and y, are the co-ordinates of the pole. 

Using centered space differences and de- 
noting by the subscript ij the value of a 
quantity at the point (i, j), we derive the finite 
difference analogue of (16), 
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Ney, = hy eyes | 
= 
ie) (18) 
de 
where 


Ay (F ost | (5) | (2) 
PNG) mas) Lo NE 
dz Oz 

ie an 7 eae = 


Jun 2) = 
ge 

4 (As)? [+ 

ier u) ehr 


Here 


— zy) 


The solution of (38) for the boundary con- 
dition 


VE 
So (| he 
ot oj 
P= 0: Zoj4+1—~ Zoj—1 | as 
(ks 
HOME 
> (2). 
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N 
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is then given explicitly by 


DE 
AINE TE 


-q) (9) 


à 2 (A 3) p—I q—I p—I g—1 
Gi i PRE > € ; 
—ıI à 
(sin ne + sin? =) (5 
2 q It J nc 
MIT sms . nli . mmf 
sin sin r sin (20) 


) 
The boundary values of = require special 
C 


attention. Here we make 


use of the second 


condition in (17): if fluid is entering the rect- 


dE 
angle, we set — = 


dt 


7) 
o; if fluid is leaving, 98 is 


ot 


determined by the interior values of 7 and z 
In the latter case we shall agree to extrapolate 


Fe 


CAS : 3 3 
Er linearly from the interior. This leads to the 


following scheme: 


4 (= d Ë 
AU" It} yp—2j 


10% mw ln | : es 
a6) —— = 


fae 2 
PNA it ) 

dË (5) 

ot en (&) iq—2 


(21) 
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The corner points are exceptional, but as they 
are not required in the computation they need 
not be considered. 


A 3 dË Oz 2 
Having determined (-——) and{—) from 
ot ji; dt 6 
(18) and (20) we perform the time extra- 
polation by means of the formulas 


Le fi = IE 
fue + At) = és((— At) +2At (5), (t) 
1] 


à Dz 
Zij (É + A t) = Sh (f— At) + SUN L ( ). (t) 
ij 


except at the first step where uncentered time 
differences must be used. The entire process is 
repeated # times if a forecast is desired for the 
time n At. 


V. The Computational Stability of the 
Finite Difference Equations 


It the finite difference solution is to approxi- 
mate closely the continuous solution, As and 
At must be small in comparison to the space 
and time scales of the physically relevant 
motions. But this does not alone insure 
accuracy; the small-scale motions for which 
there is inevitably a large distortion may 
possibly be amplified in the course of com- 
putation to such an extent that they will 
totally obscure the significant, large-scale 
motions. In the following we shall derive the 
criteria which insure that no such amplifica- 
tion will take place. We shall employ a 
heuristic procedure which is, however, pat- 
terned after the rigorous method of COURANT, 
FRIEDRICHS, and Lewy (1928). 

Since we are concerned with the behavior 
of small-scale, high-frequency perturbations, 
we may linearize the vorticity equation and 
assume the coefficients to be constant in a 
region which, though small, is yet sufficiently 
large to contain several of the small perturba- 
tions. If 2’ denotes a small perturbation su- 
perimposed on a smooth large-scale motion, 
the vorticity equation may be written 


(A =) =) (7, z) + h] (A 2,2) + 
+ /(h,2)Az 


I 
dt 
(22) 


and we have to consider the computational 
stability of its finite difference analogue: 


[= (ER Aryl A A ee 


4 h dz ZUR — A; ahs 
oy DNS 
dz À;; Hr 2 — 4j 12’ r 
— HS ( et) Tih, ZAG 2 


in which the coefficients of the terms in z’ 
may be considered constant. If, as before, the 
boundary is rectangular, the perturbation may 
be expanded into a finite Fourier series and 
it is enough to take an individual harmonic of 
the form 


a = elkxruy+n) 
Sa ’ 


where 
al al . 
ba Fe = Elana en 
m am di 
Largan ie ar 


» is the frequency, which may be compJex, and 
L, and Ly are the x and y dimensions. of the 
grid rectangle. Substitution into (23) then 
gives, after some manipulation of terms, 


o—ow t 


= JU (25) 


2 
where 


@ = e®At, (26) 


We a Ws 
— (= sin k As — — sin wh s)— 
As 2 


0 Y OX. 


2) Vig 
I sin kAs— 1 sin ws 
DS it OX dy 


4 PAS in? HAS 


2 


(27) 


Clearly the disturbance will not amplify 
if |o| <ı. To investigate the conditions 
under which this is so, let us first consider (25) 
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without the term J (h, z) At. Put @o, a root of 
this equation, equal to oe'®, where o is positive 
and © real. |wo| = @ <1 is required for sta- 
bility. Since — wo! = — gt e—® is a root 
along with wo, @ = 1 is also required. Hence 
0 = I, @ = e'®, and 


Do Do = eo — e—? = Bi sin ©), 
and the stability condition becomes 
(2) 


© find an upper bound for | a | we consider 
the two terns in (27). The absolute 


led by 


[a =] sin ON = 7. 


value of tue rirst is bes 


| At 
max |h V2 |= V2 — max |mv|. 
i 1 ESS 


As 
Fhe second increases indefinitely in numerical 
value with diminishing k aud «4, but it must 
be ssalled thet 


mn 
iv = EST 
r 7 
Loy Lx 
am TT 
LE = _ > 
i 7o T 


A ¢ 


Since » As is small, we may approximate the 
sine «fms in the numerator and denominator 
by Jueir arguments to obtain the upper bound 


Mere 


Heice the condition t leads to the 


tability criterion 


al 


AS op Las 
er V2 (max Imv|-+ Lee max | V7 ) 


(29) 


Consider now the effect of the term 
J (h, 2) At in equation (25). This equation 
(inc) has no multiple roots except at a = +1. 
If this is excluded, that is, if the stability 
condition is strengthened to | a |< 2, then 
the extra term in question will merely cause 
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a change in w of the order At. Over a 
time f, that is, over t/ Af steps, this will lead 
to an accumulation of errors which does 
not increase indefinitely for Af approaching 
zero. That is, this does not cause an amplitica- 
tion of error that can vitiate the convergence 
of an approximate solution to the true solu- 
tion for At approaching zero. Thus, at most 
an irrelevant modification of the stability 
criterion may be called for; the exclusion of 
the equality sign in the inequality that de- 
fines a lower bound for As/ At. 

Because of the extreme approximative 
character of the derivation of (29) this stability 
criterion can be regarded only as a rough 
directive in the selection of As and Ar. In 
addition, the second parenthetical term on 
the right hand side of (29) originates from a 
term in a that is proportional to At and hence 
is comparable to the term J (h, z) At whose 
neglect was advocated above. Finally the 
size of the L which enters here is debatable: 
To, use the.full:L, (.e:,..L,...L,) seems hs 
called for, since phenomena on the largest 
possible scale are not aimed at by these con- 
siderations. A detailed discussion of these 
factors would take us too far from prin- 
cipal subject. The actual values used for Ax 
and Ar were chosen on the basis of a com- 
bination of the above principles and general 
physical considerations and were ultimately 
justified when the computation was found to 


be stable. 


VI. Presentation and Analysis of Results 


The forecasts were computed for a period 
of 24 hours. The time interval used was at 
first one hour but was increased to two and 
then three hours when it was found that the 
larger intervals gave practically identical 
forecasts and did not lead to computational 
instability. The space interval As was taken 
to be 736 km, or 8 degrees of longitude at 
45 degrees latitude on the map, and the grid 
rectangle consisted of 15 x 18 space intervals. 
An actual grid is shown in figure 1. 

It will be seen that the distances are so great 
that good resolution cannot be expected for 
any but the largest scale motions. Unfortun- 
ately a smaller interval would too greatly 
have reduced the size of the forecast area, for 
the total number of grid points was restricted 
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by the limited internal memory capacity of 
the Eniac. It may be of interest to remark that 
the computation time for a 24-hour forecast 
was about 24 hours, that is, we were just able 
to keep pace with the weather. However, much 
of this time was consumed by manual and 
I.B.M. operations, namely by the reading, 
printing, reproducing, sorting, and interfiling 
of punch cards. In the course of the four 
24 hour forecasts about 100,000 standard 
LB.M. punch cards were produced and 
1,000,000 multiplications and divisions were 
performed. (These figures double if one takes 
account of the preliminary experimentation 
that was carried out.) With a larger capacity 
and higher speed machine, such as is now 
being built at the Institute for Advanced 
Study, the non-arithmetical operations will 
be eliminated and the arithmetical operations 
performed more quickly. It is estimated that 
the total computation time with a grid of 
twice the Eniac-grids density, will be about 
*/, hour, so that one has reason to hope that 
RICHARDSON’s dream (1922) of advancing the 
computation faster than the weather may soon 
be realized, at least for a two-dimensional 
model. Actually we estimate on the basis of 
the experiences acquired in the course of the 
Eniac calculatons, that if a renewed system- 
atic effort with the Eniac were to be made, 
and with a thorough routinization of the 
operations, a 24-hour prediction could be 
made on the Eniac in as little as 12 hours. 

Insofar as possible, weather situations were 
chosen in which the changes of interest 
occurred over North America or Europe, the 
areas with the best data coverage. It must be 
borne in mind, however, that forecasts for the 
western coasts are reduced somewhat in ac- 
curacy by lack of data in the Pacific and 
Atlantic oceans. 

The data were taken from the conventional 
500 mb analyses of the U. S. Weather Bureau 
and were accepted without modification in 
interpolating for the initial values of z at the 
grid points. It was realized that the conven- 
tional analyst pays more attention to wind 
direction than to wind speed and more atten- 

tion to directional smoothness of the height 
contours than to their spacing, but it was 
thought that the more or less random errors 
introduced in this way would be smoothed 
out in the integration. Unhappily this was 


A typical finite-difference grid used in the 
computations. A strip two grid intervals in width at 
the top and side borders and one grid interval in width 
at the lower border is not shown. 


Bier. 


not always so, and it now appears that an 
objective analysis would have been preferable. 

The forecasts were made from the 0300 
GMT soo mb charts for January 5, 30, and 31, 
and February 13, 1949. On these dates the 
weather systems were of so large an amplitude 
that their development could not have been 
explained by small perturbation theory or 
by simple translation. Each forecast is illustrated 
by four diagrams (fig. 2—5): (a) contains the 
initial height contours in units of 100 ft and 
the initial isolines of absolute vorticity in 
units of 1/3 x 1074 sec-!; (b) contains the 
observed height contours and constant absolute 
vorticity lines 24 hours later; in (c) the contours 
of observed height change in hundreds of feet 
are shown as continuous lines and the contovs 
of predicted height change as broken lines; 
(d) shows the predicted height contours and 
constant absolute vorticity lines. 

The spurious boundary influences were 
removed, by excluding from the drawings a 
strip adjacent to the boundaries which was 
two grid intervals in width at the west, east, 
and north boundaries, and one grid interval 
in width at the south boundary where the 
influence velocities were smaller. 

The forecast of January 5, in which the 
principal system was an intense cyclone over 
the United States, was uniformly poor. The 
forecast gave much too small a displacement 
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Fig-02; 
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of the cyclone and also distorted its shape, and 
the predictions of the other motions were 
equally inaccurate. On the other hand, the 
January ‘0 forecast contained a number of 
good features. The displacement and amplifica- 
tion of the trough over the United States at 
about 110° W was well predicted, as was the 
large scale shifting of the wind from NW to 
WSW and the increase in pressure over eastern 
Canada. The displacement of the axis of the 
major trough over the eastern United States 
and Canada was correctly predicted, but the 
strong circulation that developed at its southern 
extremity was not. Proceeding eastwards we 


d 
Forecast of January 5, 1949, 0300 GMT: (a) observed = and n at t = 0; (b) observed 


z and » 


(broken 24-hour height change; (d). 


lines) 


24 hours. The height unit is 100 ft and the unit of vorticity is 1/3 x 10-4 sec, 


find that the amplification of the trough over 
the North Sea together with the characteristic 
breakthrough of the northwesterly winds and 
the corresponding destruction over France of 
the eastern nose of the anticyclone was pre- 
dicted approximately. This is shown by the 
agreement of the predicted with the observed 
height changes over western Europe. On the 
next day the forecast was even better; the 
continued turning of the northwesterly winds 
over the North Sea and their extension into 
southwestern Europe was correctly predicted. 
The major discrepancy was the appearance of a 
sharp anticyclonic ridge south of Newfound- 
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land, which is also indicated by the excessive 
observed pressure rises over Newfoundland. 
We note, however, that the position of the 
center of these rises was correctly predicted. 
On February 13 the major changes occurred at 
the west coast of North America and in the 
Atlantic Ocean and were consequently difficult 
to forecast and to verify. Too great considera- 
tion should not be given to the Atlantic 
forecast, since data for this area were virtually 
non-existent. 

An attempt will now be made to account 
for the errors in the forecasts. The success of 
such an attempt will, of course, depend upon 


Forecast of January 30, 1949, 0300 GMT. (See Fig. 2 for explanation of diagrams.) 


one’s ability to separate the computational and 
analysis errors from those which were due to 
defects in the model. This unfortunately will 
not always be possible. Because of the excessive 
size of the space increments, the computational 
errors were in some instances obviously so 
great that nothing definite could be said 
about the residual errors due to the short- 
comings of the model. The January 5 forecast 
was a case in point. Here the grid interval was 
not at all small in comparison with the scale 
of the systems, and one had to expect a large 
computational error. This was not equally 
true of the remaining forecasts where the scale 


J. G. CHARNEY, R. FJORTOFT, J. VON NEUMANN 


Fig. 4. Forecast of January 31, 1949, 0300 GMT. (See Fig. 


was in general considerably greater. Errors 
in analysis also contributed to the difficulty in 
interpreting the results of the forecasts. These 
will be minimized as far as possible by confin- 
ing the discussion to areas in which the analysis 
was fairly reliable. Ultimately, however, the 
importance of the analysis errors can be ac- 
curately judged only by making a series of 
forecasts with the same data but with varying 
independent analyses, both subjective and 
objective. In view of the above mentioned 
difficulties the following discussion must be 
regarded as highly tentative. 

1. Truncation errors. (Errors due to the re- 
placement of the [strict] differential equation by an 


2 for explanation of diagrams.) 


[approximant] difference equation.) — Examples 
of unsystematic errors due to truncation are 
found by comparing the verification map for 
January 30 (fig. 3 b) with the synoptically iden- 
tical initial map for January 31 (fig. 4a). The 
difference in position of the grids is reflected 
in differing absolute vorticity patterns. Because 
of the large scale character of the systems the 
discrepancies are on the whole not great. 
Some are connected with the slight uncertainty 
in the actual drawing of the absolute vorticity 
lines. However, it is seen that there is a large 
discrepancy in the two patterns around the 
low off the west coast of Portugal: on the 
one map absolute vorticities of greater than 
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7 units were measured, whereas on the other 
the greatest measured absolute vorticity was 
but 4 units. 

According to the basic assumption under- 
lying the computations one should expect to 
find that the absolute vorticities “move with 
the fluid”. Merely by inspecting the map for 
January s it becomes clear that it is impossible 
to account for some of the major changes in 
the pattern of absolute vorticity around the 
low over the United States from the advection 
of absolute vorticity. Truncation errors have 
here to a large extent influenced the results and 
are probably the main reason for the strange 
computed deformation of this low. 


Forecast of February 13, 1949, 0300 GMT. (See Fig. 2 for explanation of diagrams.) 


The assumption of the conservation of 
absolute vorticity requires in particular that 
the extreme values of the absolute vorticity 
remain unchanged. An inspection of the 
individual minima brings out the fact that 
nearly all of them intensified during the 
computation (a similar computational error 
did not occur for the maxima). These dis- 
crepancies account for some of the main errors 
in the forecast. Thus on January 31 the maxi- 
mum height increase, which is found over 
Newfoundland, was about 600 ft. too great 
but corresponding to this there was a rather 
large decrease in absolute vorticity that cannot 
possibly be accounted for by advection. A 
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similar fictitious decrease of vorticity on Feb- 
ruary 13 off the east coast of the United States 
may very likely have been responsible for the 
exaggerated height increase in the Atlantic. 

A comparison between computed and ob- 
- served displacements of well-defined troughs 
and ridges gave the result that out of 20 
troughs 16 moved too slowly and the rest with 
about the right speed, and of 10 ridges 5 
moved too slowly, 3 too fast, and 2 about right. 
Thus on the whole the disturbances were 
computed to move too slowly. The dis- 
crepancy is apparently not due to truncation 
error since the finite difference approximation. 
does not systematically underestimate the 
wind speed. Nor is the explanation to be 
found in the geostrophic assumption. The 
geostrophic wind is a poor approximation to 
the actual wind in regions where the height 
contours are strongly curved; but the winds 
are overestimated in cyclones and under- 
estimated in anticyclones, and this is just the 
contrary of what must be true if the computed 
motions are to be explained. A possible 
explanation is given by the fact already 
mentioned in Section II that the level of 
nondivergence is probably higher than soo mb. 
If this level were appreciably higher, the speed 
of propagation would increase significantly 
because of the greater speed of the wind. 

2. Errors due to non-fulfillment of the vertical 
wind variation assumption. — The approxi- 
mation (1) that the winds are parallel at all 
heights is admittedly a crude one. An indi- 
cation of how the motion is modified when 
this assumption is dropped is given by the 
following consideration. 

We replace (1) by the synoptically more 
tenable assumption that the isotherms are 
parallel at all heights. In particular, we assume 
that v may be written 

V = V + Vy, 


(30) 


where v is the mean wind, supposed equal 
to the actual wind at p = p, and vr is the 
thermal wind referred to this level. 

With slight approximation the vertically 
averaged vorticity equation may be written® 


d > A 
am 7 In (¢ +f) & 0, 


6 CHARNEY (1948). 
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or, since f> ¢, except near the center of an 
intense cyclone, 


(31) 


Inserting the expression (30) and noting that 
by definition vy = 0, we obtain 


a 


dt =—v-VEtfl—w VG 


(32) 


where £, is the vorticity of the thermal wind. 


If v, were parallel to v we should again 
obtain equation (3), but the interesting case is 
where the isothermal pattern is out of phase 
with respect to the streamline pattern. In this 
case the last term in (32) produces vorticity 
changes which lead to an increase or decrease 
in the total kinetic energy. In the special case 
of a symmetric trough with a north-south 
axis, this term gives an increase of cyclonic 
vorticity along the axis when the isotherm 
trough is displaced towards the west. The 
increase of absolute vorticity in the major 
trough near the east coast of Canada on January 
30 may thus be attributed to a phase lag of the 
isotherm pattern which was observed at 
0300 GMT. As this development was the 
principal discrepancy in what was otherwise 
a good forecast, an attempt is made at the close 
of this article to give a quantitative explana- 
tion in terms of a simplified baroclinic model. 

3. The effects of baroclinicity. — The most 
typical conditions for the breakdown of the 
barotropic model will probably occur when 
potential energy is converted into kinetic 
energy. This usually is reflected in the inten- 
sification of the extrema of absolute vorticity 
which, as remarked earlier, is inconsistent 
with the assumption of conservation of 
absolute vorticity in a horizontal non-divergent 
flow. Therefore, an inspection of the observed 
changes at the extreme points in the field of 
n can be expected to give some information 
as to what extent the model has applied in 
the present situations. As a further test the 
initial and verification maps may be consulted 
to ascertain whether areas enclosed within 
isolines of absolute vorticity were actually con- 
served, and whether such regions if initially 
simply connected remained so in the observed 
developments. The February 13 forecast offers 
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two typical examples of the failure of the model 
in this sense. In the trough over the east coast 
of the Hudson Bay the maximum absolute 
vorticity was initially between s and 6 units 
but became, after 24 hours, closer to 7 units, 
corresponding to the observed, but not 
predicted, intensification of the trough. The 
other example is the trough over the west 
coast of the United States. Here the region 
with the isolines of absolute vorticity labeled 
by 4 was split into two separate regions in the 
observed development, and this splitting was 
confirmed even more markedly when the 
absolute vorticities were recomputed for a 
finer grid. The corresponding destruction of 
vorticity was responsible for the cutting off of 
a closed low which the computations failed 
to give. 

4. Barotropy versus baroclinicity. — It should 
be clear that the fact that vorticity may be 
created or destroyed cannot annihilate the 
physical effects from the horizontal transport 
of vorticity. The problem is to determine the 
relative importance of these effects in compari- 
son with those resulting from the action of 
forces creating or destroying circulation. The 
following points of view may be put forth in 
the light of the results of the present forecasts. 
Considering first the barotropic effects, it may 
be said that these were seldom if ever negligible 
in magnitude compared with the observed 
changes and were often the predominating 
ones. Moreover the transport of vorticity 
did not contribute merely to the translational 
propagation of systems but also to their de- 
velopment by the dispersion of kinetic energy. 
The forecasts for January 30 and January 31 
are illustrations of this point. 

As to the effect of the non-conservation of 
absolute vorticity, the evidence, both from 
the observed developments and from the 
failure of the forecasts in a number of instances 
where the error could not be attributed to 
the use of finite differences, supports the com- 
monly accepted view that baroclinicity cannot 
be ignored even in day to day changes in the 
upper flow patterns. Even when small, the 
non-barotropic effects may be responsible for 
important structural changes in the velocity 
field, as for example, in the cutoff low on 
January 13. However, this importance is not 
everywhere and at all times the same. Baro- 
clinic effects were apparently decisive in 
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determining the development of the major 
trough over eastern Canada and the United 
States on January 30 as the trough moved 
into the Atlantic, but on the following day 
the greatest part of the change was accounted 
for by the transport of the absolute vorticity. 

In addition to the independent effects of 
barotropy and baroclinicity there is also the 
problem of their mutual interaction. An 
intrusion of vorticity in one locality, by chang- 
ing the velocity field, automatically changes 
the conditions for the advection of vorticity 
in the surroundings. In this sense the addition 
of even small vorticities may be important, 
especially in regions of intense extrema of 
absolute vorticity. These regions generally 
possess large amounts of kinetic energy, and 
small changes in the vorticity field may 
produce large changes in the vorticity advec- 
tion. The motion of the intense low over the 
United States on January 5 was possibly an 
example of this effect. From an inspection of 
the initial map it is seen that the height contours 
and isolines of absolute vorticity coincided to 
a high degree, whereas on the observed map 
24 hours later the center of absolute vorticity 
was displaced slightly to the south of the center 
of the low pressure, thus giving more favorable 
conditions for the displacement of the whole 
system. Actually it was found from the 
observed map for 1500 GMT that the center 
speeded up in the second half of the 24-hour 
period. 

We may add that while baroclinic effects 
may be important in changing the condi- 
tions under which advection of absolute 
vorticity takes place, barotropic effects may 
distort the mass field and thereby influence 
the conditions under which potential energy 
is converted into kinetic energy. It is not, 
however, within the scope of this paper to 
discuss such effects, as little light can be thrown 
upon them by a study of purely barotropic 
processes. 

5. Suggestions for the improvement of the 
barotropic forecasts. — Strictly speaking the 
quasi-barotropic equation that has been used 
applies to the vertically averaged motion of 
the atmosphere. If, therefore, the actual 
motion at the level whose motion most closely 
approximates the mean has superimposed upon 
it other motions that do not appear in the 
vertical average, these motions will not be 
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governed by the equation of mean motion 
and will consequently exert a distorting 
effect on the forecasts. Hence, if it is possible 
in some approximate sense to ascribe to these 
motions a quasi-independent behavior, it 
would seem preferable to forecast the mean 
motions themselves and then to identify the 
mean motions with those at a particular level, 
rather than to operate with the motions for 
this level from the beginning. Thus it is 
suggested that errors due to motions that do 
not satisfy (1), as well as analysis errors, may 
be reduced by defining the z* in 


ZX 08 
Auen (8 Ae +h =") (33) 


dt f 
by 
AX =. ,A® 7,80 
Pa — Ze zd 34 
A A po J P ( ) 


instead of by z at the level p*. 

It is, of course, obvious that the length of 
the space interval should be reduced in future 
forecasts, but care must then be exercised in 
selecting the method of interpolation of the 
grid values. If the interpolation were performed 
subjectively, decreasing the size of the grid 
interval would lead to increasingly large 
errors in the difference quotients because of 
the increasing difficulty in determining small 
differences by subjective estimation. One 
might argue that a very small interval should 
not in any case be used because it would 
exaggerate the noise motions. But the way 
to smooth these motions is not to use a grid 
interval that distorts the large-scale motions; 
the smoothing is much more efficiently and 
exactly accomplished by fitting the data by 
polynomials or other mathematical functions 
which can be chosen in advance to give any 
desired degree of smoothing. A reduction in 
the size of the grid interval would then give 
more rather than less accuracy in approxi- 
mating a derivative by .a finite difference 
quotient. > 

6. À simplified baroclinic model. — In an 
effort to assess the relative importance of pure 
advection of absolute vorticity as against the 
circulation producing forces of the atmosphere, 
the following baroclinic model has been 
adopted because of its simplicity and adapta- 
bility to numerical analysis. It is also put 
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forward as the next in a hierarchy of models 
whose study is expected to lead to a better 
understanding of the atmospheric motions. 
The essential simplification has been brought 
about by the geostrophic assumption and 
the assumption of horizontal advection of 
potential temperature. As it turns out, the 
reasoning is much in accord with that of 
SUTCLIFFE (1947, 1950). 

If pressure is adopted as a vertical co- 
ordinate, the hydrostatic equation may be 
written 

OZ I 


RT: (35) 


from which one derives 
Of dea I 7 In # 
Op (5 ott. wi Outta 


where # is the potential temperature. Integra- 
tion from p,, the pressure at the ground, to 
an arbitrary level p and application of the 
surface spherical Laplacian operator to both 
sides of the resulting equation then gives 


Oz (Oz 
A, (=) as! (5) +4,A (37 


where the subscript o denotes a surface value, 
and 


(36) 


à 
I dln? 
A=— | — ——— dp. 
ge ot “FP 
Po 


Taking the vertical pressure average we get 


dz dz 
s — | = à — ele 
À (=) A. i: | 


/0z dz 
The terms A, (5) and A, (=) may now be 
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dt ot 


eliminated between (37), (38), 
the form 


to give 


NUMERICAL INTEGRATION OF THE 


(eS —_ NS 


Fig. 6. Barotropically and baroclinically computed 

soo mb height tendencies at 0300 GMT, January 30, 

1949. The barotropic tendencies are represented by 

continuous lines and the baroclinic tendencies by broken 
lines. The unit is 100 ft/24 hours. 


The advective hypothesis is finally introduced 
by ignoring w in the approximate adiabatic 
equation 


and substituting the resultant expression for 
à In #/0t in A. 

Equation (40) is then seen to be a generaliza- 
tion of the quasi-barotropic equation (3), for 
if the assumption (1) is made, streamlines an 
isolines are everywhere parallel, and A = A 


at the level p. It also can be shown to reduce 
to (32) under the assumption (30). Its solution 
is obtained as before by solving a Poisson’s 
equation in two dimensions. 

Equation (40) has been applied to the com- 
putation of the initial height tendency for the 
January 30 situation as a means of explaining 
the major discrepancy in the barotropic 
forecast for that date. The vertical integrations 


needed for evaluating J (7, z), and A were 
based on data obtained from the 1000, 850, 
700, 500, 300, 200, and 100 mb charts. 

It was to be expected that the advective 
hypothesis would be found untenable in the 
stratosphere because of the great statical sta- 
bility there. Indeed, it turned out that the 
tendencies were greatly improved by ignoring 
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Fig. 7. The broken lines represent the 24-hour height 

change computed by translating the baroclinically com- 

puted tendency field shown in Figure 6 in the direction 

of the mean current and with the speed of the trough. 

The solid lines represent the observed 24-hour height 
change. 


entirely the local potential temperature changes 
in the stratosphere, i.e., by assuming that the 
change in potential temperature due to hori- 
zontal advection is exactly compensated by 
its vertical transport. On this basis the initial 
height tendencies were calculated for all levels 
up to 300 mb for a region containing the 
trough at the east coast of Canada and United 
States. The 500 mb tendencies are shown in 
fig. 6 together with the computed initial 
tendencies obtained from the barotropic 
model. The units are in hundreds of feet per 
24 hours. It is immediately apparent that the 
baroclinically calculated height falls on the 
east side of the trough are more favorable, 
both with respect to intensity and position, 
for the observed development. In order to 
obtain a more direct comparison with obser- 
vation, the computed baroclinic tendency 
field was translated for 24 hours with the past 
speed of the trough in the direction of the 
mean current. The resulting 24-hour height 
change, together with the observed height 
change, is shown in fig. 7. It is seen from 
the diagram that the correspondence is close 
and by comparison with fig. 3 is much 
better than the barotropically computed 
change. 

It may also be mentioned that the computed 
baroclinic tendencies for the 1000 mb level 
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corresponded reasonably well with the ob- 
served tendencies on the sea-level map. The 
pressure falls to the northeast of a well de- 
veloped surface cyclone were, if anything, 
somewhat too great. An effort is now being 
made to see whether the effect of vertical 
motions is to reduce the falls. Some preliminary 
calculations indicate that this is the case. 
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On the Solar Work of the Stockholm Observatory 
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Abstract 


A brief report is given on the regular solar work of the Stockholm Observatory started in 
1947. Besides scientific research in special fields the purpose of this work is an attempt to furnish 
Swedish geophysical institutions and the Royal Swedish Board of Telegraph with solar data 
in addition to the information given in the Ursigrams. In consequence of its high geographic 
latitude Sweden is exposed to aurora zone disturbances, and therefore a local and more 
detailed solar service has been held in high estimation in our country. The observations 
may be of a general interest, however, as they show that a fairly complete record of many of 
the principal features of solar activity can be obtained with a modest equipment. Some special 
methods of observation have been developed which are well adapted for this kind of work. 


A. Sunspots and faculae. Every clear day at 
least one solar photograph is taken in the H and 
K region of the spectrum. For these pictures a 
Schott and Gen. U. G. 2 filter of 1 mm thick- 
ness is used selecting a wavelength range of 
200 A. With Ilford Lantern contrasting -plates 
such photographs reveal many more details 
than images taken in say the yellow region of 
the spectrum. The spots are very dark and 
outstanding objects in the ultraviolet in con- 
sequence of their low temperature compared 
with that of the photosphere. Moreover, the 
bright faculae can be traced not only close to 
the limb but all over the surface of the sun. 
In fact, the ultraviolet pictures resemble spectro- 
heliograms taken in the outer wings of the K 
line and in metallic lines and it seems that the 
general depression of the radiation curve in the 
H and K region by accumulated line absorp- 
tion produces the characteristics of the low 
chromosphere though with less contrast than 
on a real spectroheliogram. Among details 
present on the plates dark markings sometimes 
of a filamentary structure resembling the dark 
flocculi on Kr spectroheliograms are of a 
particular interest. Fig. 1 shows such dark 


flocculi crossing bright ones near a small spot 
close to the limb. It seems to the writer that 
this filamentary dark structure should if 
possible be included among the solar phe- 
nomena studied for the purpose of judging the 
general solar activity. 

The use of the ultraviolet spectrum has 
another advantage when taking solar photo- 
graphs and that is the improved resolution of 
the image. Though the Stockholm Observa- 
tory pictures are made with an objective hav- 
ing an aperture of only 30 mm and a focal 
length of only 800 mm they show a consid- 
erable number of details. With this small 


‘aperture the effect of seeing is not so im- 


portant as with large instruments. By using 
enlargements by eyepiece solar images with 
a diameter of 56 mm are obtained. 

B. Chromospheric faculae, flares and filaments 
in the Ho line. Since July 1948 a small spectro- 
heliograph built in the instrument shop by 
S. LINDGREN has been in operation at the 
Stockholm Observatory. The instrument con- 
sists of a grating spectrograph (f = 1 700 mm) 
of Littrow type with two symmetrically placed 
slits, one on each side of the tube. The whole 
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Fig. 1. Ultraviolet photographs of the sun’s limb obtained 

in the H. and K region on April 15, 1950. Though the 

spectral region transmitted by the filter has a width of 

about 200 A dark and bright flocculi are well visible. The 
times of exposure are given in U.T. 


spectrograph is moved perpendicular to the 
incident beam of light and the speed is con- 
trolled by driving weights and an oil-damper. 
The grating was ruled by Professor Sieg- 
bahn’s laboratory and has 720 grooves/mm 
and a very strong concentration of light in 
the first order. 

The solar image used for taking spectro- 
heliograms is as small as 22 mm. Also an image 
of only 16 mm has sometimes been tried. In 
spite of such a small size the spectroheliograms 
show a good number of details and it is easy to 
classify the sunspot groups as active ones or 
not active ones. Flares can in general be well 
detected and Hx filaments can be easily seen 
as well as prominences outside the limb. A 
spectroheliogram of the whole disc is generally 
obtained in about 20 seconds. 

Fig. 2 shows a very typical flare of in- 
tensity 3 observed with this instrument on 
May 3, 1950. The observation at 9.56 U.T. 
was started because of a call from the Royal 
Swedish Board of Telegraph informing us 
of a fade out in the short wave com- 
munications which started at 9.40 U.T. The 
conditions were nearly normal again at 10.16 
U.T. The striation in the images is mainly 
due to some dust on the slit. 

C. Prominences. For a detailed study of 
prominences a polarizing H « monochromator 
is used. It is an instrument of Lyot’s type 
independently built by the writer in 1937. The 
monochromator is attached to a 190 mm 
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refracting telescope with a focal length of 
2,500 mm. 

A new monochromator is planned to be 
used on a coronagraph now being constructed. 
The coronagraph will be used also for spec- 
trographic recordings of the solar corona. 

D. Reports and forecasts distributed to Swedish 
institutions. Twice weekly reports as to sunspots 
and general ‘solar activity are mailed to a 
number of Swedish institutions. The areas of 
the sunspots are measured and the groups are 
given a + sign when being surrounded by 
H« faculae (flocculi) and a — sign when no 
noticeable bright H« faculae are seen with 
our instrument. The measures are presented 
in a table giving ‘‘day-sectors’” before and 
after the central solar meridian and north and 
south of the equator (compare fig. 3). In the 
table one and the same region of the sun is 
represented by one and the same column and 
in this way the change in the size of the spots 
can easily be traced. At the same time days 
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Fig. 2. A typical flare of intensity 3 observed at the Stock- 
holm Observatory on May 3, 1950. 
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giving no observations in consequence of a 
cloudy sky are the less disturbing as an inter- 
polation between the values on clear days can 
generally be made without difficulty. 

In table I and table II an extract of the tables 
for the Stockholm Observatory sun spot 
areas are given for 1948 and 1949. The areas 
are measured in 10°? of the half sphere of the 
sun as unit and the figures refer to the day- 
sector after the sun’s central meridian observed 
at about 8:30" a.m. U. T. (This sector is to the 
right of number 1 in fig. 3). The upper number 
in the tables refers to the north hemisphere and 
the lower number to the south hemisphere. 
When a figure is given within parentheses no 
observation was possible that day but the value 
has been obtained by extrapolation. If this 
is the case for several days in succession the 
extrapolation has been extended over several 
days. Zero means no spot in the day-sector 
after the meridian and — that a long period of 
cloudy weather has made the determination 
impossible or very uncertain. A figure without 
+ or — sign means that no spectroheliogram 
has been obtained for the spot in question. 

Table I and table II are fairly complete, at 
least for the spring and the summer, and may 
present useful material to geophysicists. Re- 
ference is here made to previous reports from 
the Research Laboratory of Electronics, Chal- 
mers University of Technology, Gothenburg, 
where some of these figures have been used. 
Sunspot numbers or areas are generally not 
published from other observatories for one 
day-sector only and the present material may 
become handy if one wants to study the 
possible influence on terrestrial phenomena of 
some special spot. 

The solar work reported here has been 
carried out in collaboration with Mr L. 
DAHLMARK (1948—1949) and Mr H. KRISTEN- 


Fig. 3. Ultraviolet photograph of the sun obtained on 
April 25, 1950. The different sectors correspond to the 
rotation of the sun in 24 hours. 


SON (1949). The latter was in charge of the so- 
lar section during the writer’s stay at the High 
Altitude Observatory in Boulder, Colorado, 
in 1949. Both these scientists have taken a very 
active part in the observations and reductions 
and contributed with several improvements 
in the methods of observation. The work has 
been supported by grants given by the Swedish 
Council of Scientific Research. 


Stockholm Observatory, Saltsjöbaden. 
June, 1950. 


YNGVE OHMAN 
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Annual and Diurnal Temperature Variations in the 


Upper Atmosphere 


By H. WEXLER, U. S. Weather Bureau, Washington, D. C. 


(Manuscript received 20 March 1950) 


Abstract 


An attempt is made to determine the annual and diurnal variations of temperature to a height 
of 80 km by knitting together available observations and theoretical results. The observations 
include radiosonde data, meteor-determined densities, V-2 rocket observations, anomalous 
propagation of sound, and infrared spectroscopy. The theory most heavily relied upon is the 
radiative computation made by Gowan. For the annual temperature range there is reasonably 
good agreement between theory and observations, such that in the upper ozonosphere (30 to 
55 km), summer temperatures are higher than winter temperatures by amounts as large or 
larger than than these at the surface. As in lower levels, the annual temperature range aloft 
increases with latitude. For the diurnal temperature variations there is a large disagreement 
between the computations of Gowan and Penndorf. However, available evidence appears to 
favor Gowan’s computations, which, together with other considerations, indicate a diurnal 
temperature range of the order of 10 to 15° C. 


The magnitude of the annual and diurnal 
variations of air temperature below about 20 
km are well-known from years of observa- 
tions gathered from a widespread aerological 
network. For altitudes above 20 km there 
exist only sporadic observations, both direct 
(sounding balloons, rockets) and indirect (me- 
teors, sound, spectroscopy, etc.) taken at a 
few places. Despite the scarcity of data there 
have been attempts (Martyn and Puttey, 
1936), (WARFIELD, 1947) based on indirect 
observations to construct the average tempera- 
ture-height curve to over 100 km, which 
compared quite well with later direct observa- 
tions taken by the V-2 rocket soundings in 
New Mexico (DURAND, 1949). This verifica- 
tion of a result derived from speculative 
reasoning based on a few observations, en- 
courages a similar attempt to arrive at some 
idea of the annual and diurnal temperature 
variations of the upper atmosphere to a 
height of 80 km. This height is chosen because 


it represents the rough upper limit of the most 
abundant observations, both direct and in- 
direct. 

The problem of determining the annual 
and diurnal temperature variations of the 
upper atmosphere is more than just an attempt 
to extend aeroclimatological statistics to greater 
heights than attained by routine soundings. 
For, just as many important dynamical 
consequences follow from the annual and 
diurnal march of temperature in the surface 
layer of air, so it follows that significant 
dynamical influences may also result from 
annual and diurnal heating and cooling of the 
upper portion of the ozonosphere (30 to $5 
km) if it can be shown that the magnitudes 
of the temperature variations are appreciable. 
For example, in a search for possible non- 
cyclic solar controls of sea-level weather it has 
been suggested that the ozone layer, through 
its strong absorption of certain wave-lengths 
of ultraviolet radiation, may serve as the 
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Fig. 1. January and July vertical cross-sections of the average temperature field, Northern Hemisphere, in ° C. 


connecting link between sporadic bursts of 
intensified ultraviolet radiation from the sun 
and the changes in the circulation pattern of 
the terrestrial atmosphere. If it can be shown 
that the absorption of normal intensity ultra- 
violet solar radiation by the top layers of the 
ozonosphere yields only insignificantly small 
annual and diurnal temperature ranges, then 
it would be difficult indeed to expect that a 
short burst of ultraviolet radiation of double, 
or even triple, normal intensity could cause a 
significant warming up of these same layers 
and resulting important dynamical conse- 
quences on the general circulation. On the 
other hand, if annual and diurnal temperature 
variations turn out to be large, then the pros- 
ect for demonstrating an important sporadic 
ultraviolet effect on sea level weather-is more 
promising. 
Annual Temperature Variations in the 


Lower Atmosphere 


First, let us examine the annual course of 
average temperature in the lowest 19 km of 


the Northern Hemisphere from latitudes 10° 
to 70° N, as illustrated in average vertical 
cross sections for January and July (fig. 1). 
These data represent the final product of 
analysis of all available pre-World War II 
aerological data for the Northern Hemisphere 
(1945). The surface data were taken from 
HANN-SÜRING (1939). The work was done 
during the war at the U.S. Weather Bureau 
under the leadership of R. C. Gentry with 
the assistance of Air Weather Service weather 
officers. A description of these, as well as 
other war-time meteorological chart prod- 
ucts, is given elsewhere (WEXLER and TEP- 
PER, 1947). 

In fig. 1 the isotherms for each 5° C are 
shown as thin lines with thicker lines repre- 
senting the tropopause. The latter lines rep- 
resent the latitudinal variation in mean 
height of the tropopause taken from FLOHN 
(1947). Their presence is not to be inter- 
preted as implying that the tropopause is a 
continuous surface extending from tropical 
to polar regions. Indeed, it is well-known 


JULY MINUS JANUARY TEMPERATURE (°C) 


Fig. 2. Range of average temperatures from winter 
(January) to summer (July), Northern Hemisphere, in ° C. 


that there may be more than one tropo- 
pause over a given point during certain 
meteorological conditions and during other 
types of situations it may be missing al- 
together. 

The annual range of average temperature as 
approximated by the difference of the average 
temperatures, July—January, is shown in fig. 2. 
In general, the temperature range decreases at 
all altitudes from north to south. There appears 
to be distinct evidence of two layers where 
minima in the temperature range occur with 
increasing altitude. The first is at 4 to 5 km, 
and is especially distinct at higher latitudes. 
The second is found in a sloping layer, 
roughly corresponding to the average height 
of the tropopause. An attempted explanation 
of these two interesting minima is outside 
the scope of this paper; for our purpose, it 
is significant to note that at great heights 
(above 15 km) there is some evidence of 
an increasing temperature range with height 
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Fig. 3. Points: densities computed from meteor observa- 


tions (Whipple, Jacchia, Kopal). T is the normal surface 
temperature for the observation day at Boston, Mass. 
Curve: NACA density-height curve (Warfield). 


at all latitudes. The problem is to see if 
this is so, and to find the magnitude of the 
temperature range at heights above the level 
of direct observations. 


Computed Annual Temperature Range for 
an Isothermal Atmosphere 

In seeking to extend the annual range to 
greater heights, use is made of the seasonal 
variation of densities as determined by me- 
teors observed over Massachusetts (WHIPPLE, 
Jaccuta, Korat, 1949), and the observations 
are reproduced in fig. 3. 

Taking the observed annual range of the 
densities at a mean height of 78 km over 
Massachusetts as A logo © = 0.5 (loc. cit. pp. 
154—5), let us convert this value to that of 
the absolute temperature of the atmosphere 
from o to 78 km under the assumption that the 
atmosphere is isothermal at temperature To. 
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Fig. 4. Annual temperature ranges, observed and com- 
puted, expressed as ratios of the absolute temperature, 
summer to winter. 


Using the well-known hydrostatic formula 


me 
| aes Poe re () 
nee 
ER (2) 


where p is barometric pressure at height z, 
Po Is pressure at z = 0, 
R is the universal gas constant, 
is the acceleration of gravity, 

m is the molecular weight of air taken equal 
to 28.97, the sea-level value, since the 
composition of air at 70 km is the 
same as that at sea level (CHACKETT, 
PANETH, and WILSON, 1949), 

o is air density. 


Now since average sea level pressure at 
Massachusetts is exactly the same (1016 mb) in 
summer as in winter, differentiate (2) loga- 
rithmically holding py and z (78 km) constant, 
thus obtaining 


ae ) 
HR NE, (3) 
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where Ho = RT )/mg is the height of the 
homogeneous atmosphere corresponding to 
surface temperature, T). 

do 


From (3) it is seen that < 0 for z < H, 


do 
ge ovtor 22H 


Since H,~8 km and since Whipple’s 
densities are given for an average height of 
z = 78 km, it is thus seen that do/dT, > o. 
This means, if we accept Whipple’s result 
that the summer density at 78 km is higher 
than that for the winter, that the temperature 
of the lowest 78 km of the summer atmosphere 
is higher than that of the winter atmosphere. 
Let us find this difference in temperature. 
Rewriting (3) as 


CL dino 2.3 *dlog:00 (4) 
To CD Te m 
Ay Ay 


and taking dlog,,e = 0.5 at z = 78 km from 
Whipple's data, Ty =1264° Kho (Hy = 7272 
km), we find that: dT,/T, = 0.1265, dT y= 
334° andı Tal In ew. Pay where Par gy 
are the summer and winter values, respectively, 
of the temperature of an assumed isothermal 
atmosphere extending from o to 78 km over 
Massachusetts. 

In fig. 4 the ratio of 1.13 is drawn as a 
straight line from o to 78 km. Also in this 
figure are drawn two curved lines: (a) a dotted 
curve extending from o to 19 km over Mas- 
sachusetts and showing the ratio T,/T,, as 
computed from aerological observations (Ano- 
nymous, 1945); and (b) a solid curve showing 
the ratio T,/T,, calculated from Gowan’s 
computations (Gowan, 1947 a) of the vertical 
temperature distribution satisfying radiative 
equilibrium for an atmosphere at 50° N 
characterized by 0.03 % by volume of CO:, 
a constant water vapor mixing ratio corre- 
sponding to 10 % humidity at the tropopause, 
and the average distribution of ozone at that 
latitude, with winter having a greater ozone 


1 264° K is the harmonic mean temperature of the 


. NACA tentative standard atmosphere (WARFIELD, 1947 


from o to 78 km. If the arithmetic mean (273° K) is used, 
then dT = 36° and Tos/ Tow= 1.14. 
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content, 0.28 cm (NTP), compared to summer, 
0.20 cm. In both seasons the atmosphere is 
considered to be irradiated by a 6,000° K 
black-body sun. 

Gowan’s computations were based upon 
division of the atmosphere into 5 km layers 
from 55 to 15 km, and a bottom layer from 15 
km to 11 km (the latter figure being the 
average height of the tropopause at 50° N). 
He gives as his final result the mean tempera- 
ture of each layer separately for winter and 
summer, under conditions of radiative equilib- 
rium. It was a simple matter, then, for the 
present writer to draw a smooth curve through 
the ratios of these various mean temperatures 
to achieve the result shown in fig. 4. 

The ratio, T,/T,,, as found from Gowan’s 
results, increases to unity at 14 km, a height 
s km lower than the increase found in the 
observed curve over Massachusetts. This dis- 
crepancy which may be due to a difference in 
the heights of the tropopause is not important 
for the present discussion. What is important, 
however, is that above a certain altitude, both 
curves, the computed and observed, increase 
to unity, an indication that at the higher alti- 
tudes summer temperatures are probably 
greater than winter temperatures. 

Since some observational evidence has been 
found indicating that the ultraviolet intensity 
is lower than that expected from a 6,000° K 
sun, Gowan also computed the radiative 
equilibrium temperatures for a summer at- 
mosphere irradiated by a 4,000° K sun. In this 
latter case the temperature at s50—s5 km is 
347° K, or 100° less than for the 6,000° K sun; 
while at 11—15 km the temperatures are almost 
equal, 215° K vs. 217° K (see Table II). Un- 
fortunately, no computations are at hand for 
a winter atmosphere irradiated by a 4,000° K 
sun, so it is not possible to prepare a curve 
showing T,/T, against height for this case. 
It is likely, however, that the ratio curve 
would be quite similar to that for a 6,000° K 
sun. 


Computed Annual Temperature Range for 
the NACA Atmosphere 


Let us now perform new calculations in- 
volving use of the seasonal density variations 


observed by Whipple, and operating not in an 


isothermal atmosphere assumed previously but 


H. WEXLER = 


in an atmosphere whose temperature distribu- 
tion is that of the NACA tentative standard 
atmosphere (WARFIELD, 1947). (See fig. 7 in 
this paper.) 

Let T,, be the NACA temperature at a given 
height, z, and let T, and T,, be the summer 
and winter temperatures, respectively, at that 
same height in the atmosphere over Mas- 
sachusetts where it is assumed that 
AT; Ls = die Er 


“a a 


1; = 1% ne 


A T is the annual range of temperature at 2. 
The hydrostatic relationships for the sum- 
mer and winter atmospheres, respectively, are 


78 km 78 km 
RN. PP 
R | Ir R A Ty 


Ps = Pos € ° > Pw = Pow € © (6) 


where p,, pw are summer and winter pressures 
at 78 km, and pos, Pow are summer and winter 
pressures at o km. Dividing equations (6) 
and eliminating p, and p, by the respective 
Equations of State we have 


where 0,, T; are summer values of atmos- 
pheric density and temperature, 
respectively, at 78 km, 

Ow, Ty, are corresponding values for 
winter. 


Now from (s) 


A 
yar x (8) 


m 


since A T < Ta 


Substituting (8) in (7) and solving for the 


“mean value” of A T,A T, for the layer 0 to 78 
km we have 


Th; E 
In T. In" im 
Eee s ow Ow 
Al =——~ a 
Ef ale 
R 16 
Le) 
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Evaluating the numerator of this expression 
for the atmosphere over Massachusetts, we 
have seen that In pys/pow = 0 and In ¢,/0,, = 
= 0.5 ~2:3 = 1.15. For lack of definite in- 
formation of the value of In T,,/T, we shall 
assume it to be near o since it is likely that at 
78 km, where a minimum point in the tem- 
perature height curve occurs, there exists only 
a very small annual range of temperature 
because of lack of significant solar radiation 
absorption. However, even if T, were 20% 
larger than T,, then In T,/T,, =.0.18 which 
is much smaller than In 0,/0,, = 1.15. 

The integral in the denominator of (9) is 
computed from the NACA tentative standard 


atmosphere. The final result is A T'=30°C and 


The value of T,/T,, given above is drawn 
as a straight line in fig. 4. It is seen that this 
line lies slightly to the left of that computed 
for an isothermal atmosphere and is thus closer 
to the line showing the ratio of the harmonic 
mean temperatures from II to ss km of 
Gowan’s summer and winter stratospheres. 
However, the fact that the line does lie appre- 
ciably to the right of Gowan’s mean indicates 
that the annual variation of temperature in the 
lower 78 km of an NACA atmosphere, 
using Whipple’s observations, may be even 
larger than Gowan’s computations show". 


Comparison with Temperatures Determined 
by Sound Measurements 


The annual variation of temperature in 
the ozonosphere described here is in sharp 
disagreement with results found by sound 
measurements over N. Germany (GUTENBERG, 


ı From additional meteor observations the seasonal 
density variation turns out to be /\logro @ = 0.39 
instead of 0.5 (see Harvard College Observ. and Cen- 
ter of Analysis, MIT Tech. Rep. No. 2, Ballistics of 
the upper Atmosphere, L. G. Jacchia, p. 17, Cambridge, 
Mass. 1948). This will decrease AT to 23° and 
i Te et: L100; which is thus very close to Gowan’s 
value. 


267 


1946). In figure 2 of Gutenberg’s paper, 
curves are presented showing that in the layer 
20 to 60 km, the summer temperatures are 
lower than those for winter, with the maximum 
difference of 90°C occuring at 40 km. It is 
likely that this is a false result, arising from 
failure to take into adequate account the wind 


’ difference between summer and winter in the 


European stratosphere. 


Wind observations of smoke from artillery 
shell bursts at 100,000 ft. (30 km) taken over 
southeast England in 1944—45 (JOHNSON, 
1946; MURGATROYD and CLEWS, 1949) show 
that there is a monsoonal reversal from 
westerly winds in winter ( + 37 m/sec.) to 
easterly in summer (— 13 m/sec.), a difference 
of so m/sec. Assuming that this seasonal 
difference in wind is characteristic of those 
winter and summer days on which sound 
measurements were made in Germany and 
which form the basis of Gutenberg’s curves, 
it will be shown that the wind difference can 
more than account for the apparent excess 
temperature of winter over summer at 30 
km, which in Gutenberg’s figure 2 amounts 
to nearly 50° C. 


In figure 3 of Gutenberg’s paper the observed 
sound velocities are given for the layer 20 to 
60 km; it is not clear whether these sound 
velocities were corrected for wind, and if so, 
how the winds were measured. At 30 km the 
winter sound speed is 37 m/sec. greater than 
the summer speed, a value, however, which 
is less than the so m/sec. wind speed seasonal 
difference. This may signify that the greater 
winter sound speed observed at 30 km need not 
necessarily be attributed to higher temperature, 
but may very readily be caused by wind, 
provided that Gutenberg’s sound speed meas- 
urements were made to the east of the sound 
source during both winter and summer. In- 
formation on this latter point could not be 
found. 


Let us use the above data in finding T,/T,, 
at 30 km over western Europe, T, and T,, 
being the absolute temperature in summer 
and winter, respectively. The velocity of 
sound, according to Gutenberg, is given by 
20.1 VT m/sec. Now according to our assump- 
tions, the observed velocity of sound is given 
by a summation of the above temperature 
term and a wind speed. The winter-summer 
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difference in observed sound velocity is then 
given by 


W,, —W,+ 20.1 (via age VT,) = 37 meters/sec. 


where W,,— W,= 50 meters/sec., the ob- 
served winter-summer wind difference. 

Substituting this value in the preceding 
equation, we have 


VT;—VTy = 0.65, 


showing that the summer temperature at 30 
km is greater than the winter temperature. 

If we take T,, = 218° K (the NACA value), 
then we find T,/T, = 1.09, which when 
plotted in fig. 4 falls slightly to the right of 
Gowan’s curve. Again, because of various 
uncertainties, the close agreement should not 
be taken too seriously, but at least it shows 
that the annual temperature range as de- 
termined by sound measurements over Ger- 
many is not irreconciliable with our earlier 
conclusions if account is taken of the wind in 
the manner described above. 

After the original manuscript of this paper 
was sent to the Editor, additional observa- 
tions of temperature at great heights came 
to the writer, and have been plotted on 
the original figure 4: (1) “Winds and Tem- 
peratures in the Lower Strathosphere” by 
C. J. Brasefield, Journal of Meteorology, Feb- 
ruary 1950, pp. 66—69; (2) ,,Investigation 
of Stratosphere Winds and Temperatures 
from Acoustical Propagation Studies” by 
Tr Ar Peoples, Jr and A. P./ Crary which 
will be published shortly in the Journal of 
Meteorology; (3) ‘Reflection of Sound Waves 
from the Stratosphere over India in Differ- 
ent Seasons of the year” by L. S. Mathur, 
Indian Journal of Meteorology and Geophysics, 
Vol. 1, No. 1, January 1950, pp. 24—34. 
The Brasefield data included two typical 
soundings, one made in July 1948 (his figure 
5). These two soundings were taken as 
representatve of summer and, winter condi- 
tions, respectively, for New Jersey. “Tie 
second set of observations was taken at 64° 
latitude over Alaska in the summer of 1949 
and winter of 1948. The Indian observa- 
tions were made at latitude 20° N during 
the period 30 May—2 June 1946, 16—20 
October 1946, and 29 January—ı February 
1947. The temperatures for the first period 
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were considered to be representative of sum- 
mer while those of the third period were 
taken as winter temperatures. The October 
temperatures were slightly (~ 3° C) warmer 
than the winter- temperatures. An attempt 
was made to eliminate the wind effect on 
temperatures found from the Alaskan and 
Indian sound observations by having the re- 
cording stations arrayed on perpendicular lines. 

These additional data are plotted on fig. 4, 
and form some interesting comparisons with 
the material originally there. The new Jer- 
sey point falls directly on the Gowan curve. 
The Alaskan data show a steady increase of 
the temperature ratio, T,/T,, with height 
to value of 1.28 at 50 km. The Indian 
curve starts from a much smaller value of 
T,/T, and increases very slowly with height 
to a value of 1.06 at 5so—s5 km. 

Comparison of the data in fig. 4 for low, 
middle, and high latitudes reveals much of 
interest. Taking the observations as represen- 
tative of their latitudes, it is seen that above 
20 km, at all latitudes, the temperature ratio 
increases with height, but much slower at 
20° latitude than at 64° latitude. At so km 
the temperature ratio at 64° latitude is lar- 
ger (1.28) than at 20° latitude (1.06). A 
similar increase northward is noted for the 
surface temperature ratio, 1.02 at latitude 
20° and 1.14 at 64° (Hann-Süring, 1939, 
p. 180). The cause of these parallel varia- 
tions is the same. namely, a much greater 
increase in insolation from winter to sum- 
mer at higher latitudes than at lower lati- 
tudes. For example, at the bottom of an 
atmosphere having a transmission coefficient 
of 0.8, the ratio of the direct solar radiation 
received for the warm half of the year to 
the cold half of the year is 1.4 at 20° la- 
titude and 10 at 64° latitude (Milankovitch, 
1930, p. A64). 

Although we have been dealing with 
comparisons of the temperature ratio T;/T,, 
the same conclusions would apply qualita- 
tively to the annual temperature range, 
T,—T,,. For example, at 50 km the annual 
range at 64° N latitude in Alaska is 65° C 
while at 20° N latitude in India, it is 20° C. 
At the surface the corresponding ranges are 
41° C and 6° C, respectively. 

The increasing positive temperature dif- 
ference, T,—T,,, with increased latitude is 
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consistent with the change in winds ob- 
served above 17 km in middle and high 
latitudes from westerly in the winter to 
easterly in the summer. 

Although the summer-winter temperature 
ratios are larger in the ozonosphere than 
they are at the surface, this is not the case 
for the heat energies because of the much 
smaller densities aloft. For example, the ratio 
of the heat energy gained from winter to 
summer in the ozonosphere (pressure at base 
taken as 50 mb) to that gained by the whole 
atmosphere below is 0.10 at latitude 64° 
and 0.22 at latitude 20°. 


Diurnal Variation of Temperature 


The evidence regarding the diurnal varia- 
tion of temperature in the upper atmosphere 
is less consistent than that concerning the 
annual variation. Ballard found from serial 
sounding balloons that at 25 km over Nebraska 
in July the air temperature drops 10° C in the 
8 hours from 6 p.m. to 2 a.m. and less than 2° C 
at 16 km (BALLARD, 1941). It is not possible 
to assess the significance of these findings since 
only one sounding a day was taken and not 
much was known about radiational effects on 
the early temperature measurements in the 
stratosphere. There is however some indirect 
evidence based on infrared spectroscopy 
regarding the diurnal temperature variations in 
the ozonosphere. Adel by measuring the 
absorption in the ozone layer of 9.6 u radiation 
coming from the sun at day and from the 
moon at night, and making a separate measure- 
ment of the 9.6 u emission of the ozone layer 
alone was able to determine “ozonosphere” 
temperatures for several months in the spring 
of 1948 at Alamogordo, New Mexico (ADEL, 
1949). Only one of his series of measurements 
went from night to the following day, thus 
revealing partial information on the diurnal 
variation of temperature in the ozone layer. 
The observations are shown in Table 1. 

Although the temperature varies from 
— 37° to — 41° C there is no systematic trend 
indicating a diurnal tendency. In the same 
table for comparison is shown Gutenberg’s 
summary (GUTENBERG, 1949) of the White 
Sands, New Mexico, radiosonde observations 
at the highest elevation, 27.5 km, capable of 
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Table 1. 
Adel’s Ozonosphere temperature over 
Alamogordo, New Mexico 


(Observed by radiation measurements at 9.6 u) 


Date Time Temperature 

AS AS 23: 43 Mst — 39°C 
237,95 — 41 
A12S ASE serbia poe wee 00: 10 — 4I 
UPB ey 39 
T2535 37 
12: 53 — AI 


Gutenberg’s Summary of White Sands, New 
Mexico, Raob data at 27.5 km 


(April 1947—May 1948) 


Month | Jan./Febr. March | April 
Temperature —48 42°C | —so + 4 | —44 + 6 
NIGH Ol s0.6 56 6 7 24 

Month | May | July September 
Temperature —38 7 37 + 6 42 + 4 
Nor (‘Obsiene. Gi 12 5 


summarization. White Sands is 20 miles south- 
west of Alamogordo. The average monthly 
temperatures and standard errors are shown, 
and it is seen that Adel’s late April data fit in 
very nicely between the April and May means. 
This good agreement brings up the question 
as to what atmospheric layer Adel’s tempera- 
tures really refer to. Since most of the ozone 
is found between 15 and 30 km, one might 
think that Adel’s temperatures should be 
indicative of the mean temperature of the 
layer 15 to 30 km. However, the White Sands 
mean seasonal temperature-height curves 
plotted in fig. 1 of Gutenberg’s paper show 
that the mean temperature of this layer should 
be — 50° C or about 10° colder than Adels 
values.1 On the other hand the close agreement 
between Adel’s temperatures and the White 
Sands mean temperatures at 27.5 km may be 


ı Average value of Adel’s temperatures is —40° C, 
with extremes of —33° and —48°C. An independent 
comparison by the writer of the White Sands soundings 
taken within a few hours of Adel’s measurements verified 
the relative warmth of Adel’s temperatures. 


270 


interpreted as an indication that Adel’s values 
refer to the layer 25 to 30 km. 

The White Sands temperatures at 27.5 km 
indicate a definite annual variation, March 
being coldest and July warmest of the months 
shown. The ratio of the absolute temperatures! 
of these two months is 1.06, which when 
plotted on fig. 4 falls almost exactly on the 
curve drawn from Gowan’s computations 
described previously. In view of the large 
standard error of the temperatures at these 
extreme altitudes where data are sparse, too 
much significance should not be given to this 
close agreement. 

Returning to the diurnal temperature varia- 
tion it is indicated by the single case from 
Adel that the variation, if any, is very small 
in the layer from 25 to 30 km. This result 
is in agreement with Gowan’s computations 
of the nocturnal cooling of the ozonosphere 
shown in fig. 5. 

Two cooling curves are given in this 
figure: A and B. They both apply to the 8-hour 
cooling (summer night) of a stratosphere of 
constant mixing ratio corresponding to 10 % 
relative humidity at the tropopause, and 
containing 0.280 cm of ozone. The initial 
temperature distribution of each stratosphere 
varies, however. Stratosphere A has the initial 
temperature distribution corresponding to 
radiative equilibrium under a 6,000°K sun 
while Stratosphere B had the initial temper- 
ature distribution corresponding to radiative 
equilibrium under a 4,000°K sun. These 
initial temperature distributions are shown 
in Table I. 


Table II. Initial Vertical Temperature Distribu- 
tions, Gowan, Penndorf (in °K) 
7 Gowan ; 
ae si Penndorf 
A en 

5035 | 448 | 347 50 | 300 
45—50 424 333 47 280 
40—45 397 329 41.5 261 
35—40 332 301 36 253 
30—35 295 282 30.5 250 
25—30 272 266 DIES 242 
20—25 249 245 | 
15—20 2132 229 17 233 
II—IS 217 215 | 


ı Solar radiational errors were eliminated to within 
1° C by use on many of the White Sands radiosondes of 
the new unshielded white-coated thermistor (Brasefield, 
1948), so that the annual variation indicated is be- 
lieved to be close to that of the true air temperature. 
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Fig. 5. Full curves: Ozonosphere cooling after an 8 hour 
summer night (Gowan). Curve A applies to a warmer 
initial atmosphere than curve B. Upper scale in abscissa 
applies to these curves. 
Dashed curves: Daily cooling (curve C) and heating 
(curve D) in the ozonosphere (Penndorf). Lower scale 
in the abscissa applies to these curves. 


The cooling curve for Stratosphere A was 
taken from Table 1, Column 1 of Gowan’s 
paper (Gowan 1947b) while the cooling 
curve for Stratosphere B was obtained in a 
private communication from Gowan (7 March 
1950). They both show cooling of 1°C at 
25—30 km, a result in good agreement with 
Adel’s observations. 

The lack of appreciable diurnal tempera- 
ture variations at 2$—30 km may also be 
considered in good agreement with Penndorf’s 
cooling rates (PENNDORF, 1936) computed in a 
different manner from Gowan and also 
shown in fig. 5 as Curve C. This cooling 
curve represents the cooling in °C per day 
for an atmosphere containing 0.26 cm of 
ozone of vertical distribution corresponding 
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to the summer atmosphere over Arosa (lat. 
45° N). The vertical temperature distribution 
on which Penndorf based his cooling calcula- 
tion shown in Curve C, increases from 233° K 
at 17 km to 300° K at so km. 

The two cooling curves of Gowan and Penn- 
dorf, while agreeing in shape and showing 
relatively little cooling in lower levels, depart 
widely aloft; in the layer 40 to so km where 
both computations show the maximum values, 
the cooling rates differ by a factor of more than 
1,000. Many reasons may be advanced to 
account for this major discrepancy; for 
example, both computations are quite sensitive 
to the vertical distribution of ozone, water 
vapor and temperature. The two distributions 
of ozone adopted are apparently not much 
different. Although Gowan took into account 
the loss of heat in the infrared bands of water 
vapor and carbon dioxide, Penndorf did not; 
however, Gowan states that the effect of the 
latter gas is negligible. Most important appears 
to be the effect of the vertical temperature distri- 
bution assumed at the beginning of the cooling 
period; the results of two differing temperature 
distributions (Gowan A and B, Table II) are 
illustrated by comparing curves A and B in 
fig. 5. As may be seen from Table II, Penndorf’s 
temperatures are appreciably lower than 
Gowan’s throughout the layer 30 to so km. 
A cooler atmosphere of course should be 
expected to cool less rapidly than a warmer 
atmosphere, other things being equal. But it 
is not likely that this difference in initial tem- 
perature distributions can account for the 
large difference in cooling rates. Even if 
Penndorf had assumed an initial temperature 
of 500° K at so km, his cooling would still be 
100 times smaller than Gowan’s value at that 
height as explained below. 

In fig. 5 is plotted Curve D which is the 
heating curve in °C per day for the same Penn- 
dorf atmosphere described previously and 
irradiated by a 5,910°K sun at lat. 45° on 
June 22. The heating values shown by Curve 
D are much larger in magnitude than the 
cooling values shown in Curve C, a result 
which indicates lack of radiative equilibrium. 
However, for radiative equilibrium to exist, 
Penndorf states that a temperature of at least 
500° K would be required at 50 km (loc. cit. p. 
267). This would then mean that at 500° K 
Penndorf’s cooling rate would be equal to 
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the heating rate, but would still be too small 
compared to Gowan’s results by a factor of 
at least 100. 

In absence of a more positive explanation 
of the discrepancy, we may try to see which 
of these two theoretical results agrees better 
with observations. The observations will be 
of two types: Those based on what is believed 
to be known of the annual temperature 
variations at the top of the ozone layer (40 to 
55 km) and those based on the temperature 
lapse-rate between ss and 80 km. 

Referring to Gowan’s nocturnal cooling 
values in fig. 5, it may first be said that for his 
results to be physically reasonable, the daytime 
heating values must be of the same magnitude, 
but of opposite sign, as the cooling values; for 
otherwise the air at these levels would be 
rapidly cooling or heating, and would soon 
reach impossibly low or high values. 

The net warming up from winter to sum- 
mer occurs as the days grow longer and the 
nights shorter, thus allowing a small positive 
temperature gain to accumulate each day. 
Adopting the Gowan computations — both 
seasonal and diurnal — one would expect at 
40 to 45 km an annual temperature cycle of 
maximum range of 75°C on which is su- 
perimposed daily oscillations of range at least 
10 to 12°C. In contrast to Gowan’s values, 
Penndorf’s results if plotted on a similar 
graph would show practically no diurnal 
temperature variation at any place in the 
ozonosphere. Let us estimate from the full 
value of Penndorf’s daily heating rate the 
maximum annual range of temperature to 
be expected in the layer 20 to so km. Penn- 
dorf’s average daily heating for this layer at 
the spring equinox is about 0.12°C per day 
(loc. cit. p. 266). In 180 days (i.e., from mid- 
winter to midsummer) this layer would then 
warm about 22° C. From the earlier discussion. 
on annual temperature variations, as illustrated 
in fig. 4, it appears likely that the average 
annual range for the layer 20 to so km is 
larger than 22°. Thus even if non-radiative 
heat dissipation, such as advection and con- 
vection, are neglected and the full amount of 
Penndorf’s daily radiative heating value is 
used, the expected annual temperature varia- 
tion cannot be accounted for. It appears from 
this consideration that Penndorf’s heating 
rates are too low. 
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Fig. 6. Comparison of measured upper atmosphere tem- 
peratures (Helgoland blast), V-2 rocket results, and 
NACA tentative standard (after Cox). 


But what proof have we that Gowan’s 
large daily temperature oscillations of magni- 
tude at least 10° C exist at the top of the 
ozonosphere, from 40 to 55 km? In other 
words, is it not possible for the winter to 
summer warming of the top layers of the 
ozonosphere to occur as a day-to-day accumu- 
lation of small differences of two much smaller 
quantities, say 1/1oth of Gowan’s values as 
shown in fig. 5. Arithmetically, of course, this 
solution is possible, but there is additional 
evidence which lends some support to the 
magnitude of Gowan’s cooling values. 

If a diurnal temperature variation of at 
least 10° C is observed at 55 km and zero tem- 
perature variation exists at 25 km, this would 
mean that the temperature inversion found 
between 30 and 55 km, and illustrated in fig. 
6 (Cox, 1948), would be strengthened by 
day and weakened by night, but remain as 
an inversion day and night. What can one 
say of the layer from 55 to 80 km? First, it 
appears that in the upper portion of this layer, 
say, from 70 to 80 km, there is relatively little 
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known absorption of solar radiation. Below 
70 km absorption by ozone heats the air with 
the maximum heating between 40 and 55 
km, as described previously. Thus the layer 
55 to 80 km is analogous to the layer of air 
above the ground which is heated by insola- 
tion at day, forming a steep, or even a su- 
peradiabatic lapse rate in the surface layer, 
and cooled at night, forming an inversion. 
The absorbing layer at the top of the ozone 
layer is not nearly so sharply bounded as 
the ground in its absorptive properties, but 
the analogy is fairly good. Thus, if during 
the daytime heating period the temperature 
at 55 km should increase by at least 10° while 
the temperature at 70 to 80 km remains un- 
changed, one might expect that a steep, 
possibly even a superadiabatic lapse rate, might 
be established in the lower portion of this 
layer. 
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Fig. 7. Relation of temperature to altitude above White 
Sands, New Mexico (after Newell). 
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To support the contention of the existence 
of a steep or superadiabatic lapse rate in the 
layer 55 to 80 km, we can cite tempera- 
ture measurements made by the V-2 rocket 
soundings released at 1123 MST on 7 March 
1947, at 1313, 22 January 1948, and 1837, 
5 August 1948 at White Sands, New Mexico; 
the mean sounding of which is shown in fig. 7 
(NEWELL, 1950). These observations show a 
temperature fall of 75° C between 60 and 66 km, 
or a superadiabatic lapse rate amounting to 
125% of the dry adiabatic rate. For this lapse 
rate to become adiabatic, the temperature at 
60 km would have to be decreased by 15° C. 
No night rocket soundings at which air tem- 
peratures have been measured are available 
for comparison with those made at day. Cox, 
in his determination of stratospheric tempera- 
tures from the Helgoland blast of 1,100 GCT, 
18 April 1947, was forced to assume the exist- 
ence of a superadiabatic layer from 57 to 62 
km to render his data interpretable (fig. 6). 

It should be pointed out that there exists 
considerable uncertainty in the determination 
of the air temperatures from the V-2 rockets. 
Because of the high speed of the rocket, tem- 
peratures are not measured by a temperature 
element, but by two methods, not independ- 
ent, which are based on air pressure. The first 
method uses the slope of the log p vs. altitude 
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curve to compute the mean temperature of 
5 or 10 km thick layers. The second method 
is to determine the density from the speed 
of the rocket and the ram pressure (pressure 
measured at the nose of the rocket), and then 
by use of atmospheric pressure to find the air 
temperature by the gas law. For a further 
discussion of the fairly large probable errors 
(~ 20° C) involved in the temperature values, 
see a paper soon to be published in the Journal 
of Geophysical Research (HAvENS, Koti, La Gow, 
1950) and kindly loaned to the writer by the 
authors. 

In concluding this section on diurnal temper- 
ature variations, it may be stated that while 
evidence points to practically zero diurnal 
temperature change in the lower part of the 
ozonosphere (up to 30 km), evidence is not 
so convincing as to the magnitude of the change 
in the upper portion of the ozonosphere 
(from 30 to 55 km). It is likely that Penndorf’s 
values are too low and that the much higher 
values of 10 to 15° C, as deduced from 
Gowan’s computations and considerations of 
lapse rate and annual temperature range, 
are nearer the truth. It would be most 
helpful to the study of the diurnal tempe- 
rature variations in the ozonosphere if serial 
day and night rocket soundings could be 


made. 
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II. The Climatology of Blocking Action 
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Abstract 


The results of a semi-statistical study of 112 cases of blocking action in the upper 
westerlies are presented; areas of most frequent occurrence, characteristic movement 
and persistence, and seasonal and yearly trends of blocking activity are determined. The 
regional precipitation and mean surface-temperature anomalies associated with blocking 
over Europe are illustrated and compared with those produced by strong zonal flow 
aloft. The effects produced by variations in blocking activity upon European climatic 
trends are discussed. The influence of blocking action upon Scandinavian glaciation 
and certain indications as to the röle played by blocking in an ice-age climatic cycle 


are discussed. 


Introduction 


In Part I, previously published, the results 
of a detailed analysis of five selected cases of 
strong blocking action were presented. The 
former study provides a description of at- 
mospheric blocks at mid-tropospheric level 
and some understanding of the initiation, 
development and dissipation of this phenom- 
enon. The reader is referred to Part I, Rex 
(1950), for a description of the concept of 

locking action as used herein. 

In order to assess the climatological signifi- 
cance of blocking action, which is the primary 
aim of this paper (Part II), it becomes at once 
desirable to investigate the characteristics of 
the phenomenon more generally, i.c. statisti- 
cally in so far as data limitations permit. Further- 
more, before any analysis of the effect of 
blocking upon regional climate can be under- 
taken, it is necessary to determine the specific 
climatic anomalies associated with blocked 
periods. Accordingly, the immediate objective 


of this paper is two-fold; first, to compile a 
catalog of blocked cases from which the areas 
of most frequent occurrence, normal move- 
ment, seasonal and yearly trends, etc. may be 
statistically computed and second, to determine 
the precipitation and mean surface-temperature 
anomalies observed over the region dominated 
by blocking action in typical cases. Utilizing 
these results an attempt will be made to relate 
blocking variability to recently observed 
yearly climatic fluctuations in Europe and to 
discuss its rôle in connection with present day 
glacier regimes and with past ice-age climatic 
cycles. For the latter purpose, it is necessary 
to investigate the climatic influences of a 
zonal flow-pattern aloft over Europe. 


Statistical Characteristics of Blocking Action 


The concept of blocking used in the compila- 
tion of the blocking catalog is based upon the 
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development of a particular contour pattern 
(circulation pattern) in the troposphere. In 
view of the fact that large-scale features of 
the atmospheric motion are more easily 
identified at upper levels, away from complex 
surface structures, it will be most desirable to 
focus principle attention upon the lower- 
middle troposphere. Any attempt to identify 
blocked cases with certainty from an inspection 
of the longer available record of surface anal- 
yses would require a prohibitive expenditure 
of time. In order to obviate the necessity for 
time consuming subjective analysis, and to 
enjoy the security of an objective method 
one may introduce a definition of blocking 
based upon the occurrence of persistent high 
pressure within certain “key” areas (ELLIOTT 
and SMITH 1949). In so doing, the selection 
of the “key” areas limits the generality of the 
results obtained for the purposes of this study. 
Perhaps more important, it seems likely that 
certain pressure developments will satisfy 
a definition based upon persistent high pressure, 
which are not examples of blocking in the 
sense used here. Unfortunately, useable middle 
tropospheric analyses are available only for 
recent years. Nevertheless it seemed best to 
accept this limitation rather than those imposed 
by the selection of one of the alternate tech- 
niques mentioned. 


If one employs both 3,000 dynamic-meter 
constant level and soo mb constant pressure 
charts, published daily analyses over the north- 
ern hemisphere are available for about ten 
different years. Similar charts covering the 
longtitudinal sector from 160° W through o° 
to 90° E are available for an additional three 
years. It must be emphasized, however, that 
the analysis is uncertain or altogether missing 
in certain areas on the older charts due to the 
sparsity of data; this is particularly true over 
the western Pacific and eastern Asiatic con- 
tinent. In performing the search each daily 
chart was inspected by the author in order to 
identify pressure or contour patterns satisfying 
the following definition of blocking action. 
A blocked case must exhibit the following 
characteristics at upper levels: 


a. The basic westerly current must split into 
two branches; 


b. each branch current must transport an 
appreciable mass; 
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c. the double-jet system must extend over 
at least 45 degrees of longitude; 


d. a sharp transition in the westerlies from 
a zonal type flow upstream to a meridio- 
nal type downstream must be observed 
across the current split; and, 


e. the pattern must persist with recognizable 
continuity for at least ten days. 


In a situation satisfying these requirements 
blocking is said to be initiated when (a) occurs 
and is said to have dissipated whenever (a), 
(b), (c) or (d) are no longer satisfied. It is 
perhaps worthwhile to mention that the 
definition used can be rather easily expressed 
in quantitative form in terms of zonal and 
meridional circulation components and their 
derivatives. The less exact subjective form was 
used in view of the rather extensive computa- 
tional procedure involved in the application 
of an objective definition. 


The results of the search just described are 
shown in Table I. For each blocked case 
identified the initial, final and approximate 
mid-period date and longitudinal position are 
given. For the purposes of this paper the 
longitudinal position of a block is defined as 
the longitude of the westerly current split. 
In addition the area, type and pattern as well 
as the data source and certain remarks are 
given in the table. Because of the frequency 
distribution of the location of block initiation 
obtained (figure 2), it is possible to classify 
blocking by area as cither Atlantic (ATL) or 
Pacific (PAC). In both areas blocked cases of 
two types are distinguished in the table — those 
having a warm anticyclonic vortex (WA) 
immediately downstream from the jet split, 
and those having a cold cyclonic vortex (CC). 
The designations strong (S), average (M) and 
weak (W) used with reference to the blocked 
pattern indicate crudely the clarity of the 
pattern as compared with the “strong” cases 
discussed in Part I (June—July 1949 and 
January 1947). It should be understood that 
these pattern designations are not intended to 
measure the intensity of the blocking action. 
As has been shown in Part I blocked cases often 
appear to pass through several cycles of 
rejuvenation during their life, also Pacific and 
Atlantic cases may exist simultaneously; such 
features have been indicated under Remarks 
in the table. 
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Table I. Blocking case catalog 


Initial Mid-period Final 
NE À A es item Remarks Source 
Day/ Day/ Day/ 
Mo. Long. Mo. Long. Mo. Long. 


1932 ZEIT BOSSE ET DE | OV | 78/7 LOEW) PATE. TEE M JMC 
6/12 | 30 W| ıı/ı2| 30 W| 15/12] ro W » » W ir » 
8/12 |145 W] 14/12 |140 W| 19/12 |120 W| PAC | WA NM iif Sonne Block » 
Tossa LS). ıo W| 24/1 35 W| 28/1 40 W| ATL | WA M JMC 
Tee 35 WI 20/2 45 W| 25/2 so W » » » \ » 
23/2. 855 RB | 28/2 Iv7s R | 52/3) 1x80 PAC » NP nt N 
8/3 o 14/3 o 18/3 TOME || MAMET: » W » 
21/3 Se 323 15 W| 10/4 20 W » » M several cycles » 
13/4 |135 WI 21/4 |140 W| 28/4 |135 W| PAC] » » pute ATL . 
25/4 10 W| 2}5 22 Wi 20/5 | ro W| ATE, » S ee 7 » 
18/5 | 10 W| 3/6 | 18 W| 16/6 o » » » several cycles » 
5/7 o 10/7 ES PAPA 25 E » » M » 
26/7 |150 W| 30/7 |140 W| 4/8 |ızo W| PAC 2 » » 
18/10 | 28 W| 23/10] 25 W| 28/10] so W| ATL | CC S ae vote » 
ı6/ıı | 12 W| 27/11 | 15 W| r0o/12| 35 W| » WA Se Ball Nios ee RS » 
9221375 E | 15/12 |180 21/12 |170 E | PAC » M J en » 
1934 | 20/1 10 W| 26/1 10 W| 30/1 10 W| ATL | WA M low latitude JMC 
26/2 15 W| 2/3 By NWA © Gill 10 E » » » » 
20/3 8 W| 28/3 174 WAN A so W » CE » \ » 
20/3 |150 W| 24/3 |165 W| 30/3 |150 W| PAC » S Map block » 
27/4 5 WI 3/5 10 W| 9/5 SEE peewee » M » 
0m ES SWITCH fi-22) SW}, 126 |) 25. W| © WA W ) 
31/7 15 E 4/8 10 W| 9/8 o » CE M » 
5/12 |140 W| 11/12 |140 W| 15/12 |125 W| PAC| WA » \ » 
78/122 070, BE 24/12 | LOW IE T2 12 EB | ATE CE » ATL after PAC » 
1935 6/1 35 WI 14/1 20 W]| 22/1 30 W| ATL| WA M JMC 
1/3 10 W| 8/3 15 W| 18/3 sE » » S \ » 
nae 140 W ae TSW Leesa EL eA cl) Ce Ww PAC after ATL » 
31/3 |1ss W| 7/4 |170 W] 16/4 |130 W » » M two cycles » 
23/4 | 20 W| 11/5 | 20 W| 1/6 | ro W| ATL | WA S Nees hon » 
16/5 |120 WI 23/s |130 W| 31/5 |130 W| PAC | CC WwW iif » 
19/6 5 WI 25/6 15 WI 29/6 10 W| ATL| WA M » 
1936 1/2 |175 W 2 |160 WI 18/2 |150 W| PAC | CC W JMC 
5 As + Ww Si 5 W| 16/2 5E| ATL | WA > N aoubte nn » 
5/3 o 223 5 W| 19/3 30 W » Ge M » 
1/4 | 45 W| 6/4 40 WI 14/4 | 40 W » » » k » 
19/4 | 10 W| 6/s | 20 W| 26/5 o » WA Sn Be » 
18/4 |160 W| 24/4 |175 W| 28/4 |170 E | PAC | CC » » 
27/5 |155 W| 2/6 |140 W| 7/6 |1s0 W » ? Ww > 
18/6 | 30 W| 22/6 | 40 W| 27/6 | 30 W| ATL | CC M ) 
29/9 | 20 W| 4/10] 20 WI 12/10] 40 WI » WA » » 
6/12 o 11/12 o 227 12) Tey de » CC WwW » 
I I WI 18/1 10 BY 27) x Zone, MAD) WEA S JMC 
Ae ies Ww a 180 12/2 |145 W| PAC » » PAC after ATL N 
26/3 15 W| 10/4 10 E | 24/4 o ATL » M » 
2/5 20 W| 10/5 ıs W| ıs/s s W] » » » » 
30/5 |140 W| 9/6 |130 W|:17/6 |115 W| PAC » » » 
16/8 6 W| 23/8 o 28/8 o DE || KCXC; » » 
2/10| ro W| 9/10} 28 W| 17/10] 35 W| » » WwW » 
31/10| 20 W| 8/11] 30 W SE 5 OM » » M » 
1938 | 14/2 | 30 W| 22/2 | ro W| 28/2 | 20 E | ATL | WA S JMC 
5/4 |130 W| 15/4 |115 W| 25/4 |130 W PAC » M » 
Ww 35 WI s/s 30 W| ATL » S » 
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Table I (cont.). 
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Blocking case catalog 


ea a ce 
Final 


Year 


1939 


1940 


1945 


1946 


1947 


1948 


1949 


Initial 


Day/ Day/ | Day! | L 
Mo: Long. Mo. Long. Mo. Long 


2/8 | 25 W 
15/10 |135 W 
TO TON SEN, 

4/12] 10 E 
25/1 27 W 
20/2 TORE: 
25/3 | 25 W 
28/4 | 25 W 
14/6 SE 
17/8 10 W 
24/10| 25 W 

7/1 0 
17/1 |130 W 

8/2 10 W 
25/9 | 40 W 
24/11| 20 W 

9/12 |145 W 
17/12 | 15 W 
9/1 40 W 
15/2 5 W 

6/4 | 30 W 

3/7 10° W 

1/8 | 20 W 

3/9 | 30 W 
29/10| 18 W 

8/12} so W 
11/12 |157 W 
14/1 | 20 W 

2/3 10 W 

2/5 |160 W 
4/6 |170 W 
6/10} 10 W 
10/12] o 

1/1 0 
2/1 |155 W 
19/1 10 W 

5/2 |175 W 
10/2 30 W 
15/3 |165 W 
4/5 s E 
17/6 o 
EB MSN 

8/8 15 W 
20/10 | 10 W 
16/2 | 25 W 
26/3 25 W 
27/3 |175 W 
7/s |140 W 
9/6 |1so W 
16/12 | 20 E 
23/1 |1so W 
2310 


Mid-period | 


10/8 | 20 
19/10 |125 
24/10 | 15 
DS 121 Te 
31/1 5 
28/2 5 

2/4 5 
18/5 | 15 
20/6 0 
25/8 10 
30/10 | 25 
12/1 25 
25 TE |ITS 
722201025 
2/10} 20 

1/2130 
15/12 |125 
24/12 | 25 
13/1 50 
19/2 5 
16/4 | 25 
9/7 o 
9/8 20 
10/9 0 
din XS} 
14/12] 25 
18/12 |ıss 
22/1 10 
11/3 | 15 
12/5 |140 
18/6 |135 
20/10 | 25 
201122 RO 
6/1 ) 
10/1 |175 
27/1 30 
172 1250: 
16/2 20 
223125 
14/5 5 
28/6 0 
20/7 fe) 
24/8 | 15 
28/10 | 20 
24/2 | 40 
3/4 | 20 
6/4 |180 
19/5 |140 
15/6 |150 
2412) 25 
2/2 |165 


25 W| 31/1 40 


Ww 
Ww 
Ww 


44444 444 4 288 Seeeeee <q TPE 


EEE 


ig zz "a 4e 


18/8 

24/10 
28/10 
26/12 


6/2 
5/3 
12/4 
7/6 
25/6 
31/8 
4/11 


16/1 
30/1 
1712 
4/10 
5/12 
21/12 
30/12 


18/1 
2512 
24/4 
15/7 
16/8 
15/9 
23/11 
17/12 
23/12 


31/1 
19/3 
25/5 
25/6 
DT 
SAT TS 


12/1 
18/1 
2/2 
2/3 
222 
29/3 
26/5 
6/7 
25/7 
7/9 
SI TER 


1/3 
11/4 
16/4 

2/6 
19/6 
29/12 


14/2 
9/2 


Area 
30 BY) >ATE 
130 W| PAC 

o ANNE 
Io E » 
10 W| ATL 
TS » 

5 W » 
20 W » 
30 W » 
10 W » 

5 E » 
20 W| ATL 

115 W| PAC 
20 W| ATL 
10 E » 
Is W » 

115 W| PAC 
to W| ATE 
55 W| ATL 
20 E » 
40 W » 
15 W » 
Io W » 
20 E » 
15 W » 
20 W » 

155 W| PAC 
Ay AES || Wa IRL 
10 W » 

120 W| PAC 

125 W » 
10 W| ATL 

Oo » 

TorB || ASG 
155 W| PAC 
o AT 
165 W| PAC 
so W| ATL 

125 W| PAC 

o AIRE 
25 E » 
10 E » 

oO » 

o » 
20 W| ATL 
5 W » 
165 W| PAC 
135 W » 
135 W » 

2.50) ALE 
170 W| PAC 


TosE | ALE 


Type | Pattern Remarks Source 
@E S JMC 
» M » 

WA ‘ double block 5 
CE S » 
WA S low latitude JMC 
EC » » 
WA W » 

» S several cycles » 
CS M » 

2 WwW » 
WA S » 
WA W Il JMC 

> S f PAC after ATL > 

» W » 

» » » 

» M » 

» » \ » 

: S | ATL after PAC N 
WA W BB 

» » » 

» M two cycles » 

» » » 

» » » 
GG » » 

» S two cycles AAF 
WA » \ ree. » 

5 Ww | double block N 
WA S AAF 
CC W several cycles » 
WA M » 

» W » 

» S two cycles BB 
Ce » cont. Jan. 1947 N 
Ce S \ began Dec. 1946 BB 
WA M | double block WxME 

» S BB 

» two cycles 

ar \ double block WxMF 

» » | preceding zonal flow BB 

» » several cycles WxMF 

» » » » BB 

» » » » » 

» W » 

» S » 

» M » 
WA S BB 
» » \ » 

$ ‘ N double block WxME 
CC W three cycles » 
WA S » 

3 M BB 

? WxB 
WA ÿ double block s 
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Table I (cont.). Blocking case catalog 


| Initial Mid-period Final 
Year RE men Vee ul le Area | Type | Pattern Remarks Source 
ay ay 
Me. Long. re Long. The. | Long. 
Esel 30 WiEr8/30 | 40: Wi 4 20 Bal) ATI | WA M || several cycles > 
3/3 |160 W| 12/3 |160 W| 18/3 1165 W| PAC » S ae ee » 
16/6 17 WI 27/6 10 W| 7/7 TOMBA TL » » TWB 
8/9 o 20/9 30 WI 29/9 | 25 W » » » » 
16/11 s WI 24/ıı | 30 W| 30/ıı | 15 E » » M high latitude » 
un EEORW ly) ta |e rok. | ar/ro [Mrs E » CE W » 
19$0 | 18/1 25 W| 25/1 25 W| 2/2 TORE, | CADIS WA: S TWB 
18/3 TS EU 26/3 30 WI 1/4 50 W » » » » 
ATL — Atlantic S — strong 
PAC — Pacific M — average 
WA — Anticyclonic vortex W — weak 


CC 


— Cyclonic vortex 


See text for source abbreviations. 
Longitudes given in whole degrees. 


The following published source material 
has been used in the compilation of Table I: 
Historical Daily Weather Maps, Northern Hemi- 
sphere, 3,000 dyn. meters, published by the 
Joint Meteorological Committee, Wash., D. C. 
(JMC), Northern Hemisphere 500 MB Maps, 
published by Headquarters, Air Weather Ser- 
vice, Wash., D. C. (AAF), Northern Hemi- 
sphere Daily 500 MB Maps, published by the U.S. 
Weather Bureau, Wash., D. C. (WxB), Upper 
Air Section, Daily Weather Report, 500 mb, 
published by the Meteorological Office, 
London (BB), Taglicher Wetterbericht, 500 mb, 
published by Zentralamt Bad Kissingen (TWB), 
and upper-air analyses, soo mb, published in 
micro-film form by the U. S. Weather Bureau, 
Wash., D.C. (WxMF). 


From the data compiled in Table I certain 


characteristics of the blocking phenomenon 
at mid-tropospheric level may be demonstrated 
which have at least quasi-statistical significance 
in spite of the limited number of cases exam- 
ined. In figure ı is shown a histogram of the 
frequency distribution of the position of block 
initiation in longitude. To construct the dia- 
gram, the total number of cases initiated in 
successive longitudinal bands 10 degrees in 
width were determined and plotted as per- 
centages of the total number of cases in eac 

of the two groups; cases initiated on boundary 
longitudes were assigned one-half to each 
adjacent longitudinal band. A band width of 
. 10 degrees was chosen in view of the fact that 


the position of initiation can not be determined 
more exactly than + 5 degrees of longitude. 
The 112 cases identified are grouped into two 
rather definite preferred regions; one centered 
between 5° and 15° W longitude (containing 
82 cases) and one centered between 145° and 
155° W longitude (containing 30 cases). 
Examples of blocking action contained in the 
first group will hereafter be referred to as 
Atlantic cases and those contained in the 
second group as Pacific cases. Both the Atlantic 
and Pacific distributions are relatively sym- 
metrical with only a slight skewness toward 
the west in each case; accordingly, the median 
and mode lie close together near the center 
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Fig. 1. Geographical distribution of block occurrence 

in the northern hemisphere; the histograms show the 

frequency distribution in longitude of 112 cases of block 

initiation at mid-tropospheric level. Modes shaded black. 

Cases included occurred during the years 1933—1940 
and 1945—1949 inclusive. 
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Fig. 2. Seasonal trend in the mean longitude of Atlantic 

block initiation (82 cases) compared with the seasonal 

longitudinal movement of the mean Atlantic (Icelandic) 

low-pressure center. Positions were determined in each 
instance at mid-tropospheric level. 


of cach group. In the two preferred zones 
thus determined, 68 % of all Atlantic cases 
occur between 5° E and 25° W longitude 
and 60% of all Pacific cases occur between 
135° and 165° W longitude; a similar result 
has previously been obtained by ELLIOTT and 
SMITH ois) in their study of blocking 
development at the surface. No cases of 
block initiation were found between 55° W 
and 115° W longitude or between 25° E and 
145° E longitude. It must be emphasized, 
however, that because of data limitations some 
cases may have escaped detection. The in- 
clusion of such cases is not likely to appreciably 
alter the form of figure 1, except possibly in 
the western Pacific—eastern Asiatic area where 
the source material is uniformly inferior. In 
the author’s opinion, however, such Asiatic 
cases are infrequent, if not altogether non- 
existent. 

The appearance of two preferred zones for 
block initiation is not altogether surprising if 
one accepts the concept of block initiation 
described in Part I and considers the location 
aloft of the Atlantic (Icelandic) and Pacific 
(Aleutian) semi-permanent cyclonic centers. 
These semi-permanent cyclones tend in the 
mean to produce in the lower troposphere a 
southward movement of the cold air-masses 
lying over the North American and eastern 
Asiatic continental basins which, in turn, 
tends to produce locally a stronger and more 
concentrated jet stream aloft. Furthermore, 
such cyclonic centers represent the locus of 
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most frequent and intense “cyclonic im- 
pulses” (see Part I, p. 206). We should therefore 
expect the initiation of block development to 
be favored in the region lying downstream 
from the semi-permanent cyclonic centers. 
To illustrate the association just described 
figure 2 has been prepared showing by month 
the longitude of the mean Atlantic low center 
(at the soo mb level) and that of most frequent 
Atlantic block initiation; positions for the 
cyclonic center were taken from Normal 
Weather Maps, Northern Hemisphere Upper 
Level published by the Joint Meteorological 
Committee, Washington, D. C., and those 
for the Atlantic block were computed from 
the data given in Table I. Approximations 
for both functions were obtained by the 
method of semi-averages and are shown as 
the solid curves in the figure. The seasonal 
displacement in longitude of both the Atlantic 
cyclone and block are related although very 
slightly out of phase. Drawing rather heavily 
upon the hydraulic analogue previously dis- 
cussed we may expect, in any blocking de- 
velopment, the distance from the point of the 
initiating impulse downstream to the point of 
block initiation to increase as the mean stream 
velocity increases. In figure 2 we notice that 
the distance from the cyclonic center to the 
block decreases from a maximum in January 
or February to a minimum in August or 
September paralleling changes observed in 
the strength of the mean westerlies aloft. We 
suggest that the phase difference mentioned be- 
fore may result from the effect just described. 
Comparisons made for Pacific cases produced 
similar results but are not shown herein. It 
has been suggested that block development is 
primarily the result of the dynamic influence 
of continental land masses or of the thermo- 
dynamic effect of continental-oceanic bound- 
aries; such explanations are more difficult to 
apply to Pacific cases which are most frequently 
initiated some 30 degrees in longitude seaward 
from the North American west coast and overan 
ocean surface free from appreciable latitudinal 
thermal gradients. It must be emphasized how- 
ever that the dynamic and thermodynamic ef- 
fects just mentioned tend to support blocking 
development, and that they may be a prime 
cause for the observed stability (position) of 
blocks but not, in the author’s opinion, for 
the initiation of blocking development. 
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The 82 Atlantic cases occurred during a 
time interval equal to 13% years or at an 
average frequency of about one case every 
two months; the 30 Pacific cases occurred 
during 10 years and 8 months or at an average 
frequency of about one case every four 
months. More exactly the number of Atlantic 
cases is equal to 68 %, or slightly more than 
two-thirds, of all cases occurring around the 
hemisphere. Although this predominance of 
Atlantic cases may be in part explained by 
less adequate analyses over the Pacific area, 
a real excess of Atlantic cases no doubt remains. 
If one accepts the concept of block initiation 
given in Part I, such an excess would be 
expected in view of the relatively stronger 
and narrower westerly current system aloft 
which is observed over the North American 
continent and western Atlantic as compared 
with conditions found over the western 
Pacific. In both Atlantic and Pacific blocks 
the warm-anticyclone type accounts for 72% 
of the cases in each group. The contour 
(circulation) pattern associated with this type 
(Part I, figure 1) may then be taken as the 
most “typical” blocked pattern. Although 
the latitudinal position of the block is not 
given in Table I examples of extreme high 
and low latitude cases are so indicated under 
Remarks. Latitudinal variations in the location 
of block initiation and development are 
relatively small being governed by the lati- 
tudinal position of the mean jet aloft. Typical 
block developments are centered between 
45° and 55° N latitude. 

A seasonal trend of blocking action has been 
computed from the data of Table I and is 
shown for both Atlantic and Pacific cases in 
figure 3. The percentage of days (each month) 
during which blocking action was observed 
is shown by months. An approximation for 
each function was obtained by the method of 
semi-averages and is shown as a solid curve 
in the figure. The two trends obtained are 
approximately similar and in phase, rising 
from a flat minimum in the late summer and 
early autumn to a sharper maximum in the 
late spring. Pacific blocking activity varies 
from zero (August through September) to a 
maximum of about 23 % (April and May). 

Essentially different from a climatological 
point of view, Atlantic blocking activity 
never falls below a minimum of about 15 % 
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Fig. 3. Seasonal variation in blocking activity; Atlantic 

(82 cases) and Pacific (30 cases) trends in the percentage 

of days dominated by blocking action are shown by 
month. 


(July through September) and rises to a 
surprisingly strong maximum of about 40% 
(April and May). The occurrence of maximum 
blocking activity as late as April or May 
emphasizes the importance of stream width 
in the process. In accordance with the pre- 
viously discussed concept of block initiation, 
we should expect increasing activity during 
the autumn and winter as mean westerly 
velocities and the intensity of initiating im- 
pulses increase. A continuation of this rise into 
the late spring, long after impulse and stream 
intensity have begun to weaken, must be 
explained by a rapid decrease in width of the 
mean westerly flow aloft which is, in fact, 
observed during the spring and early summer. 

The persistence of blocked situations is 
analyzed in figure 4 which shows a frequency 
distribution of blocked cases by duration. The 
numbers of Atlantic cases, having durations 
in whole days from 10 to 41 (the longest 
found), were determined from the data of 
Table I and are shown in the figure as per- 
centages of the total number of cases. Per- 
centage values to the nearest one percent are 
used, and 4 % of the cases, those having 
periods longer than 34 days, are not shown 
in the figure. The distribution obtained is 
irregular but shows a mode at 14 days. A 
secondary, and probably superficial, maximum 
occurs at the boundary limit of 10 days and is 
most likely the result of psychological in- 
fluences operating in the application of the 
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Fig. 4. Persistence of blocking action in the Atlantic 

area; the straight-line frequency graph shows the distri- 

bution by duration of 82 cases of Atlantic block devel- 

opment. Cases (4%) with periods longer than 34 days 
are not shown in the figure. 


definition of blocking. The median coincides 
with the mode (14 days) and the average 
duration was calculated to be 16.6 days. An 
analysis of Pacific case duration shows similar 
results with a most probable duration of 12 
days. We may conclude that although blocking 
cases of shorter duration than to days most 
certainly occur, the most probable duration 
lies between 12 and 16 days. The application of 
the definition used, in so far as it pertains to 
persistence, has not excluded a significantly 
large percentage of actual cases; certainly all 
cases of climatological significance have been. 
considered. : 

The longitudinal movement of Atlantic 
blocks is illustrated in figure 5 by means of 
a series of histograms showing the frequency 
distribution of block position initially, at 
mid-period and just prior to dissipation. These 
diagrams were constructed in the same man- 
ner as that described previously in connection 
with figure 1, the uppermost (initial) being 
identical with that shown in figure 1. At mid- 
period (about 7 days after initiation) both the 
mode and median have shifted toward the 
west. The distribution remains quite regular 
but is skewed slightly toward the east. The 
final distribution (about 14 days after initiation) 
shows a tendency toward the development of 
two maxima, one located quite far to the west 
at about 50° W longitude and the other, more 
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well developed, toward the east at 10° E 
longitude. We conclude that Atlantic blocks 
typically move toward the west some 10 
degrees during the first half of their develop- 
ment. Subsequently, during their final half- 
life most blocks move eastward some 30 
degrees, resulting in a net displacement toward 
the east of some 20 degrees; a very small 
percentage continue their westerly motion 
reaching a final position some 40 degrees 
west of their initial position. The rather re- 
markable stability (position) of blocking action 
is most striking; more than 65 % of all cases 
observed remain within a longitudinal span 
40 degrees in width and centered about 5° W 
longitude. In a recent paper Bourn (1950) has 
discussed the stable perturbations produced in 
the upper westerlies by orographic barriers; 
notably the North American Rocky Moun- 
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Fig. 5. Movement characteristics of Atlantic blocks; the 

series of histograms shows the frequency distribution in 

longitude of 82 cases of Atlantic block development 
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tains. He describes the characteristics of the 
quasi-stationary long-wave pattern produced. 
It is suggested that the ever present tendency 
for secondary ridge formation, which is 
found over the North Atlantic, may signi- 
ficantly influence the movement of Atlantic 
blocks. Blocks, once formed in the area, 
would tend to move into or remain in such a 
position that the blocking and orographically 
formed ridges are in superposition. 

In order to disclose any pronounced seasonal 
variation in the movement characteristics of 
blocks, Table II was prepared. The percentages 
of all Atlantic cases each month with an 
observed net displacement toward the west, 
toward the east, and stable were computed 
from the data of Table I. Net displacement is 
defined as the longitudinal difference between 
the final and initial positions of the block. No 
definite trends are evident from an inspection 
of the table, due in part perhaps to the limited 
number of cases available. In every season net 
displacements toward the east dominate those 
toward the west. There are also indications of 
an increase in the percentage of stable cases, and 
a decrease in the percentage of westward 
moving cases, during the summer months; 
the percentage of eastward moving cases 
remains essentially constant throughout the 
year. 

In concluding the statistical study of blocking 
action a yearly trend of Atlantic blocking 
activity was computed from the data of Table 
I and is shown in graphical form in figure 6. 
A similar computation for Pacific cases was 
not attempted in view of the fact that sufficient 
data is not available. The number of days 
during which blocking was observed during 
each year, each winter (Jan.—Feb.—Mar.) and 
each summer (July—Aug.—Sept.) was com- 
puted and expressed as a percentage of the 
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Fig. 6. Yearly variation in blocking activity in the 

Atlantic area; yearly, winter and summer trends in the 

percentage of days dominated by Atlantic blocking 
action are shown by year. 


total number of days in each instance. The 
yearly and seasonal yearly trends thus obtained 
are shown in the figure for the periods 1933 
through 1940 and 1945 through 1949." It 
should be emphasized that, although these 
trends account for blocking variability with 
regard to frequency of occurrence and per- 
sistence, they do not include any factor re- 
lated to the “intensity” of the phenomenon. 
Their use in the sense of an “index of blocking” 
must therefore be made with caution. The trends 
of figure 6 show one salient difference between 
the 1930 and 1940 decades, namely, that larger 
fluctuations in amplitude characterize the more 
recent period; this is most clearly seen in the 
summer yearly trend. 


1 After completion of this manuscript an arith- 
metical error was found in the 1947 values shown 
in figure 6; the correct yearly and winter values are 
38 % and 44 % respectively. 


Table II. Block movement by months — Atlantic cases 


Percentage of Total Cases Each Month 


Net Displace- 


ward an | Feb. | Mar. | Apr. | May | June | July | Aug. | Sept. | Octzl Nova ||, Dec: 
ANESESE sr 25 22 45 25 o © o 17 33 43 34 20 
Stable". 33 33 22 37 71 57 50 so 0 43 33 40 
Se en sels 42 45 23 38 29 43 50 33 67 I4 33 40 


t Stable — net displacement less than 5° E or W. 
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Fig. 7. Mean surface-pressure distribution over the Atlantic-European area during a wintertime case of strong Atlantic 

blocking action (16—31 January 1947). Surface isobars are shown as thin solid lines; pressures given in millibars. 

The mean contour pattern at the 500 mb level is shown by the thin long-dashed lines; contour heights given in 
dynamic decameters. Cyclone tracks at the surface are indicated by the heavy, solid and dashed, arrows. 


The Effect of Blocking Action upon 
European Climate 


As has been shown in Part I, the circulation 
pattern typically associated with a blocked 
zone is markedly cellular in character and 
represents a strong departure from the nor- 
mal. In view of the stability and persistence of 
such cases (figures 4 and 5), we should expect 
blocking activity to be associated with well 
defined climatic anomalies over the underlying 
area. It is the purpose of this section to describe 
the isanomal patterns associated with Atlantic 
blocks; no attempt has been made to analyze 
the corresponding Pacific cases in view of the 


difficulty of obtaining climatic data over the 
area concerned. 


In order to obtain results representative of 
the seasonal extremes, one winter (January 
1947) and one summer (June—July. 1949) case 
have been analyzed in detail. The winter case 
selected was initiated on 19 January at 10° W 
longitude and dissipated 15 days later on 2 
February over the Greenwich meridian. Be- 
cause of the availability of certain data, the 
period discussed in the following study extends 
from 16 to 31 January inclusive. The slight 
discrepancy involved will not affect the general 
form of the results obtained. In figure 7 the 
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Fig. 8. Isanomals of precipitation over the Atlantic-European area during a wintertime case of strong Atlantic 
blocking action (16—31 January 1947). Isanomal lines are given for 10, 50, 100 and 150 percent of normal and are 
shown as thin solid lines, except in areas of uncertain analysis where they are shown as thin short-dashed lines. Areas 
recording greater than normal amounts of precipitation are indicated by the darker shading; areas recording pre- 
cipitation amounts greater than so % of, but less than, normal are indicated by the lighter shading. The mean 
contour pattern at the 500 mb level is shown by the thin long-dashed lines; contour heights given in dynamic 

decameters. Stations used in constructing the isanomal pattern are indicated by the solid dots. 


mean soo mb contour pattern for this period 
(thin dashed lines) is shown superimposed over 
the corresponding mean surface pressure 
pattern (thin solid lines). The tracks followed 
by surface cyclonic systems during the period 
are shown as heavy arrows in the figure. 
Four principal storm tracks may easily be 
identified in figure 7. The first, situated over 
the north-central Atlantic, extends from the 
seaward area south and east of Greenland over 
Iceland and on toward the east-northeast to 
Spitzbergen and the North Cape region. 
This track is associated with the strong south 


19—005140 


to southwesterly current aloft around the 
western side of the blocking-ridge (the north- 
ern current-branch of the block). Storms 
moving in this path, favored by a strong north- 
west—southeast thermal gradient and hav- 
ing’ access to an adequate supply of relatively 
warm moist air to the southeast, developed 
rapidly becoming deep slow-moving occlu- 
sions before reaching Iceland and there, often 
generating secondary systems which moved 
on toward the Spitzbergen area. From the 
standpoint of Continental climate the principal 
influence of these storm centers was felt as a 


286 


series of diffuse cold frontal passages across the 
Scandinavian peninsula and northwestern Eu- 
rope producing scattered showery precipita- 
tion (snow) over upland and coastal areas. 
This storm track is shown in figure 7 by 
heavy solid arrows. 

The second track, lying over the Medi- 
terranean area, extends from the central 
Atlantic over Gibraltar and then easterly 
across the Italian peninsula and on to the 
Balkan and Black Sea regions. It is associated 
with the strong westerly jet observed over the 
Mediterranean at 500 mb (the southern current- 
branch of the block). Four storm systems 


moved along this path during the period under: 


discussion and produced light to moderate 
precipitation (snow in upland areas) over the 
southern slopes of the Italian, Swiss, Austrian 
and Balkan Alps. This storm track is shown 
by heavy solid arrows in figure 7. 

The third track, lying principally over the 
Baltic, extends from Swedish Lappland and 
extreme northwestern Russia toward the south- 
west across northwestern Europe to the Low 
Countries. This cyclone track is associated with 
the northeasterly jet aloft which flows down 
the western side of the extended White Sea 
trough and is shown in figure 7 by heavy short- 
dashed arrows. Disturbances moving along this 
path were inverted, having a relatively warm 
maritime air mass to the northwest and a rela- 
tively cold dry continental air mass to the south- 
east. Three such low-pressure centers, moving 
rapidly as open frontal waves toward the south- 
west, produced light snowfall over Finland, 
south-central Sweden and Denmark. Any 
intense development of such disturbances, 
with resulting heavier precipitation, is pre- 
cluded by the fact that much of the moisture 
originally in the warm air supply is orographi- 
cally removed during its flow over the Scandi- 
navian mountain divide. The cyclone shown 
over the English Channel, although associated 
with this third storm track, illustrates a special 
type of development not infrequently observed 
over western Europe, but one not necessarily 
associated with blocking action. This disturb- 
ance developed following the cutting-off and 
subsequent intensification of a cyclonic vortex 
aloft over western Europe. The resulting 
stagnant surface center produced sustained 
heavy snowfall and squally weather over the 
Channel and over the southern British Isles 
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and parts of the French north coast during 
the last days of January. 

The fourth track extends from the Black 
Sea area toward the north and west over 
southwestern Russia. It is associated with the 
northeasterly flow aloft around the south- 
eastern side of the White Sea trough and is 
shown in figure 7 by heavy long-dashed 
arrows. Two surface cyclones of moderate 
intensity were observed during the period 
and produced some precipitation in the form 
of rain over the Ukrainian plains. 

The precipitation isanomal pattern com- 
puted for the European area for the winter 
period is shown in figure 8. The mean soo mb 
contour pattern for the period (thin dashed 
lines) is shown superimposed over the pre- 
cipitation isanomals (thin solid lines) in order 
to facilitate orientation. The anomalies were 
computed with respect to the fifty-year 
normal, 1881—1930, and are indicated in the 
figure by a system of isanomal lines and 
shading. Normal values were in most cases 
taken from published data given by KÔPPEN 
and GEIGER (1932 (a) and (b), 1935 and 1939), 
CLAYTON (1927 and 1934) and KENDREW 
(1922). In some instances 1881—1930 normal 
values were not available and in such cases the 
best available substitutes were used. The 
pattern obtained reflects in general variations 
from the 1881—1930 normal. In view of the 
deflection of the normal westerly storm 
tracks to the north and south, away from con- 
tinental Europe, we obtain precipitation values 
generally much below normal over the entire 
area. Only in certain upland and windward- 
coastal regions are normal or slightly above 
normal values observed. It should perhaps be 
noted that the areas of less than 10 %, of normal 
precipitation shown in western Russia and 
over the North Atlantic are valid only in the 
region immediately adjacent to the analyzed 
area. Most probably precipitation amounts 
over east Greenland, the North Atlantic and 
Spitzbergen exceeded normal. 

Figure 9 shows the anomalies of mean 
surface-temperature computed for the Euro- 
pean area during the winter period. The same 
reference normals and an essentially similar 
method of representation were used as those 
described in connection with figure 8. Mean 
24-hr. surface temperatures were used in com- 
puting the anomaly pattern. The European 
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Fig. 9. Isanomals of mean surface-temperature over the Atlantic-European area during a wintertime case of strong 
Atlantic blocking action (16—31 January 1947). Isanomal lines are given for each 2° C above and below normal 
and are shown as thin solid lines, except in areas of uncertain analysis where they are shown as thin short-dashed 
lines; the isonormal line (zero anomaly) is shown as a heavy solid (or dashed) line. Areas recording mean tempera- 
tures more than 4° C above normal are indicated by the darker shading; areas recording mean temperatures more 
than 4° C below normal are indicated by the lighter shading. The mean contour pattern at the soo mb level is 
shown by the thin long-dashed lines; contour heights given in dynamic decameters. Stations used in constructing 
the isanomal pattern are indicated by the solid dots. 


continent, British Isles and Mediterranean 
area were subjected to below-normal tempera- 
tures. These negative anomalies are associated 
with flow at the surface of cold continental 
polar air into a well-developed White Sea 
trough and in fact reach their greatest values 
(8° to 10° C below normal) below and near 
the position of this trough line aloft. A tongue 
of positive anomalies extends up the western 


side and around the top of the blocking-ridge, : 


reaching as far south as the Black Sea area and 
including most of Norway, Sweden, Finland 
and western Russia. These above-normal 


surface temperatures are located in the path of 
warm air advection around the blocking-ridge 
and reach maximum values ranging from 
10° to 12° C above normal over east Green- 
land and Spitzbergen with a secondary 
maximum of from 6° to 8° C above normal 
over central Sweden. 

We now shift our attention to the second 
case of blocking action analyzed; the summer 
case. This blocked period was initiated on 
16 June 1949 at 17° W longitude and ended 
22 days later at 10° E longitude on 7 July 1949. 
A rather detailed aerological description of 
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Fig. 10. Mean surface-pressure distribution over the Atlantic-European area during a summertime case of strong 

Atlantic blocking action (16 June—7 July 1949). Surface isobars are shown as thin solid lines; pressures given in milli- 

bars. The mean contour pattern at the soo mb level is shown by the thin long-dashed lines; contour heights given 
in dynamic decameters. Cyclone tracks at the surface are indicated by the heavy, solid and dashed, arrows. 


this period, together with illustrations of the 
development of the contour (circulation) 
pattern, has been given in Part I previously. 
Computations of mean 500 mb contour and 
surface pressure patterns and of the precipita- 
tion and mean surface-temperature anomalies 
were made and are shown, together with the 
surface storm tracks, in figures 10, 11 and 12. 
The method of computation and representa- 
tion used in the construction of these figures 
is similar to that used in figures 7, 8 and 9. 
Comparing figures 7 and 10, we note that 
both the surface and upper-air pressure 
patterns show a greater development of the 
cyclonic center over western Russia and a 
less pronounced development of the blocking 


high in the summer case. As a result the storm 
track previously referred to as the fourth is 
eee somewhat toward the north and 
carried five well-developed cyclonic systems 
as far north as central Finland, producing 
moderate to heavy rains over northwestern 
Russia. The two patterns occupy essentially the 
same relative position with, as would be 
expected, the north and south jet branches not 
so well developed in the summertime case. 
The southern jet branch in this case is displaced 
far to the south passing over North Africa 
at about 30° N latitude and carrying the Medi- 
terranean storm track (the second track) 
beyond the limits of the base chart used. 
An inspection of figure 11 shows again in 
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Fig. 11. Isanomals of precipitation over the Atlantic-European area during a summertime case of strong Atlantic 
blocking action (16 June—7 July 1949). Isanomal lines are given for 10, 50, 100, 150, 200 and 300 percent of 
normal and are shown as thin solid lines, except in areas of uncertain analysis where they are shown as thin short- 
dashed lines. Areas recording greater than normal amounts of precipitation are indicated by the darker shading; 
areas recording precipitation amounts greater than 50 % of, but less than, normal are indicated by the lighter 
shading. The mean contour pattern at the soo mb level is shown by the thin long-dashed lines; contour heights 
given in dynamic decameters. Stations used in constructing the isanomal pattern are indicated by the solid dots. 


this summer case a rather remarkable reduction 
of rainfall over the European continent gener- 
ally. We note some increase in the amounts 
of precipitation recorded in northwestern 
Russia, Poland and the Ukraine in association 
with the intensified Russian low-pressure 
center mentioned above. Along the North 
African coast rather heavy rainfall was recorded 
in contrast with the extremely low values 
associated with the wintertime case. As shown 
by figure 11, the more southerly position of the 
Mediterranean storm track produced moist 
on-shore winds along this coastline in contrast 
to the dry off-shore air flow shown in figure 8. 


Normal or greater than normal rainfall was 
recorded only over small and widely scattered 
areas on the Continent. 

A comparison of figures 9 and 12 reveals no 
significant differences between the temperature 
anomaly patterns associated with the two 
cases except a relative displacement of the 
area of negative anomalies some s or 10 
degrees eastward in the summer case. As a 
result of this displacement extreme western 
Europe and the British Isles are exposed to 
small positive temperature anomalies. As 
would be expected in view of the less intense 
summertime circulation, the anomaly extremes 
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Fig. 12. Isanomals of mean surface-temperature over the Atlantic-European area during a summertime case of 
strong Atlantic blocking action (16 June—7 July 1949). Isanomal lines are given for each 2° C above and below nor- 
mal and are shown as thin solid lines, except in areas of uncertain analysis where they are shown as thin short-dashed 
lines; the isonormal line (zero anomaly) is shown as a heavy solid (or dashed) line. Areas recording mean tempera- 
tures more than 4° C above normal are indicated by the darker shading; areas recording mean temperatures 
more than 4° C below normal are indicated by the lighter shading. The mean contour pattern at the soo mb 
level is shown by the thin long-dashed lines; contour heights given in dynamic decameters. Stations used in 
constructing the isanomal pattern are indicated by the solid dots. 


shown in figure 12 are less intense than in 
the former case. A thermopleion of 4° C 
above normal in Ireland and in the White and 
Caspian Sea areas and a thermomeion of 6° C 
below normal in Yugoslavia are observed. 
The validity of any general conclusions as to 
the climatic regime associated with European 
blocking which may be drawn from the two 
studies just discussed will, of course, depend 
directly upon how representative the cases 
analyzed are. In this connection, it should be 
noted that the winter case selected was initiated 
more toward the east than normal, whereas 


the summer case selected was initiated to the 
west of the normal longitude (see figure 2). 
This rather unfortunate happenstance is not 
likely to make any extensive reconsideration 
necessary in so far as the wintertime anomaly 
patterns over Europe are concerned, but would 
lead one to expect considerably higher tempera- 
tures over western and central Europe, the 
British Isles and Scandinavia, in conjunction 
with summertime blocking, than those indi- 
cated in figure 12. No satisfactory answer to the 
question of representativity can be given until 
sufficient blocking cases are analyzed to permit 
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the construction of a “mean” blocked precipita- 
tion and mean temperature anomaly pattern 
for the European area. The computation of 
such a “mean” pattern is at present underway; 
meanwhile, from the preliminary results so 
far obtained it seems safe to assume that the 
climatic anomalies described in the preceding 
paragraphs are, in main outline, representative 
of blocked situations in general. 

We may then briefly summarize the regional 
climatic characteristics of blocked periods 
over Europe as follows: 

1) Precipitation amounts over the Continent, 
Scandinavia and the British Isles are gener- 
ally much below normal. Normal or 
slightly above-normal values are recorded 
only in scattered small areas along wind- 
ward coasts and on windward mountain 
slopes. Precipitation amounts over the 
Scandinavian and Alpine mountain com- 
plexes are generally much below normal. 

2) Wintertime surface temperatures are 2° to 
6° C above normal over central and north- 
ern Scandinavia but are below normal over 
the entire Continent and the British Isles, 
reaching a minimum of 8° to 10° C below 
normal over the Balkan area. Summertime 
blocking may be expected to produce even 
higher surface temperatures over Scandi- 
navia and positive anomalies over the 
British Isles and westcentral Europe. 

3) Orographic and convective (within the 
polar air masses brought down over the 
Continent by the extremely well-developed 
White Sea trough) rain and snow showers 
account for a major portion of the total 
precipitation recorded over the European 
area. 

4) Moderate, mostly north-easterly, surface 
winds prevail over the Continent generally. 
Storminess, as associated with the passage 
of intense cyclonic systems, is more or less 
absent over the entire area. 

In view of the strong and rather clearly 
defined precipitation and mean surface-tem- 
perature anomalies associated with blocking 
action over Europe and because of the rela- 
tively large percentage of the time dominated 
by blocking (figure 6), one might expect 
yearly variations in blocking activity to be 
reflected in European climatic trends. Such 
would most certainly be the case unless some 
process systematically acts so as to produce a 
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reversal of the anomaly patterns following 
each blocked period. Since no reason for such 
a tendency is known it was thought desirable 
to make a comparison between the blocking 
trends shown in figure 6 (hereafter referred 
to as blocking indices) and yearly European 
mean surface-temperature and precipitation 
trends computed for the same period. As has 
been previously mentioned, the use of the 
trends of figure 6 as indices is a somewhat 
unjustified procedure in view of the fact that 
no account is thereby taken of variations in 
the “intensity” of blocking. Intensity varia- 
tions may be expected to produce changes of 
magnitude, although not of isanomal geo- 
metry, in the associated climatic anomaly 
patterns. 

Climatic trends were computed for 35 se- 
lected stations and areas in Europe and have 
been compared with the indices of blocking 
activity. The climatic trends were computed 
from data given by Crayton (1947) and that 
published in the meteorological and climato- 
logical bulletins, summaries and yearbooks of 
the various meteorological services in the 

rea. In figure 13 are shown the results of the 
mean surface-temperature trend comparisons 
for certain stations; in figure 14 the corre- 
sponding precipitation trend comparisons are 
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Fig. 13. Comparison between the yearly variations in 
Atlantic blocking activity and fluctuations in the mean 
surface-temperature recorded at various European sta- 
tions; winter season and summer season trends are shown. 
Daily mean 24-hr. surface temperatures were used in 
computing the seasonal means for each station. “K”’ 
values indicate the correlation coefficient obtained be- 
tween the blocking index and the particular temperature 
trend in question. 
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Fig. 14. Comparison between the yearly variations in 

Atlantic blocking activity and fluctuations in the amount 

of precipitation recorded at various European stations; 

winter season and summer season trends are shown. “‘K”’ 

values indicate the correlation coefficient obtained be- 

tween the blocking index and the particular precipitation 
trend in question. 


shown!. In these figures the seasonal indices 
of blocking (winter and summer) have been 
used rather than the yearly index. Those sta- 
tions selected for illustration showed the best 
correlation with the corresponding blocking 
indices but reflect, in a general way, the results 
obtained at all stations sampled within the 
adjacent area. Correlation coefficients between 
index and climatic trends have been computed 
in some cases, using a standard procedure as 
described by Conrap and POLLAK (1950) and 
have been corrected for “sample size”. 
Returning to figure 13, Berlin and Frankfurt 
are representative of the thermomeion, while 
Jan Mayen is representative of the thermo- 
pleion. The computed mean temperature 
trends correlate rather strongly, negatively and 
positively respectively, with the winter index 
trend. In the summer season, the Gibraltar and 
Vienna trends, representative of the ther- 
momeion, and the Oslo and Karesuando trends, 
representative of the thermopleion, correlate 
with the summer index trend in the same sense 
as that described for the winter case. It is 
interesting to note that in the summer season 
the positive correlations are relatively stronger 
and the negative correlations relatively weaker 


1 As before mentioned the correct 1947 blocking- 
index winter value is 44 %. The Berlin and Norr- 
land negative correlation coefficients shown in figures 13 
and 14 are very slightly improved by this correction. 
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suggesting from another point of view that a 
“typical” summertime block produces a tem- 
perature anomaly pattern at the surface which 
is shifted somewhat to the east with respect to 
that shown in figure 12. Both the winter and 
summer precipitation trends computed for all 
European stations selected show a negative 
correlation with the corresponding index 
trends; note especially the negative correlation 
found between average wintertime precipita- 
tion over Norrland and the winter index. 

The modest correlation coefficients obtained 
and the short length of the trend record 
available obviously limit the interpretation 
which may be made of the results. Never- 
theless it seems rather certain that yearly 
variations in blocking activity are directly 
associated with the climatic fluctuations over 
Europe which have been observed during 
recent years. We may speculate further in 
suggesting that a similar relationship between 
variations in Pacific blocking activity and the 
observed climatic fluctuations over western 
North America may be found. 

The correlations obtained, being uniform 
with respect to sign and systematically or- 
dered in the sense of the typical blocked cli- 
matic patterns just described, support our 
previous assumption that blocked situations 
are not followed with any regularity by 
periods of “opposite” climatic extreme. The 
occasional. occurrence however of such a se- 
quence of events is well exemplified by the 
European winter of 1947. In early December 
1946 strong blocking developed over Europe 
which persisted almost without interruption 
until the third week in February 1947. As 
we should expect, the weather over northern 
and central Europe and most of the British 
Isles during this period was extremely cold 
and relatively dry. This circulation regime 
was followed in very early March by a 
strong zonal circulation type which, as will 
be discussed in the following section, may 
be considered to be a ‘climatic opposite” 
to the blocked type. This period, ending in 
mid-April, was characterized by milder tem- 
peratures and an almost unprecedented 
amount of rain- and snowfall over the area 
in question. The combination of heavy snow- 
fall following persistent cold produced an 
extremely “severe” winter. It is of course 
clear that a season of this type (Jan. —Mar. 
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Fig. 15. Isanomals of precipitation over the Atlantic-European area during a wintertime case of strong zonal flow 
aloft (1—31 January 1949). See figure 11 for an explanation of the method of representation used. 


1947) can not be adequately described in 
terms of the typical blocked pattern and 
the blocking index since a high index value 
results from the first half-period (see figure 6), 
whereas large total precipitation amounts 
result from the second half-period. 


The Effect of a Zonal Circulation Type 
upon European Climate 

In the sense advocated by RossBy and 
“WiILLETT (1948) and by WiLLETT (1949) the 


general circulation of the atmosphere may be 
said to periodically fluctuate between two 
extreme states; one characterized by strong 
meridional motions (in quasi-horizontal planes) 
with an expanded circumpolar vortex and 
one characterized by a strong zonal motion 
with a contracted circumpolar vortex. Blocked 
situations, as discussed in the preceding section, 
produce a hemispheric circulation of the 
former type. In order to permit a comparison 
of the effects upon European climate produced 
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by the two opposed types, blocked and zonal, 
an analysis of the climatic anomaly pattern 
produced over the European area by a zonal 
circulation type aloft has been made. For 
the purposes of this comparison a zonal 
type circulation is considered to be one in 
which a single westerly jet stream, without 
strong north—south flow components, is 
observed at the soo mb level over the sector 
lying between 30° W and 60° E longitude. 
Two periods of persistent strong zonal flow 
over Europe, one occurring in January and 
the other in June, were selected for analysis 
in order to minimize the effects of seasonal 
variation in comparing them with the blocked 
cases. The winter period includes the entire 
month of January 1949 and produced the 
precipitation and mean surface-temperature 
distributions shown in figures ıs and 16. 
These figures were constructed from prepared 
data and charts for the month of January 1949 
published in Die Grosswetterlagen Mitteleuropas, 
Zentralamt Bad Kissingen, II, 8. The method 
of representation used is the same as that 
described in connection with figures 8 and 9. 
The mean 500 mb contour pattern produced 
during the winter period was essentially zonal 
over the area in question (except in the Medi- 
terranean area) with a well-developed westerly 
jet passing across the British Isles and into 
northern Europe at about 55° N latitude at 
the soo mb level. One principal storm track 
was observed during the month; it extended 
from the west across Ireland over southern 
Sweden and on toward the northeast into the 
Baltic region. Frequent deep disturbances 
following this track produced heavy precipita- 
tion (snowfall) over the entire region lying 
adjacent to and north of the track; a positive 
precipitation anomaly reaching values over 
250 % of normal was recorded in the 
Scandinavian mountains where orographic 
effects augmented frontal action. Elsewhere 
in Europe below-normal amounts were re- 
corded except on the windward (north- 
westerly) slopes of the Pyrenees, Alps and 
Caucasus where normal to twice normal 
values were observed. Further south along 
the North African coast and over the southern 
Mediterranean light to moderate frontal rain, 
produced in association with a quasi-stationary 
thermal front extending along the south 
Mediterranean coast, built up positive ano- 
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malies of 100 to 150 % of normal. Mean 
surface-temperatures over most of Great 
Britain and Europe were near normal or 
slightly above (1° to 2° C above normal) with 
an area of relatively strong positive anomalies 
(8° C above normal) situated over north- 
western Russia. This anomaly may be inter- 
preted as the result of an abnormal easterly 
penetration of the relatively warm Atlantic 
maritime air mass, which was carried inland 
by the strong westerlies at least as far as central 
Russia and there replaced the relatively cold 
dry polar air mass which normally dominates 
that region. 

In figures 17 and 18 are shown the precipita- 
tion and mean surface-temperature anomalies 


: observed over Europe in connection with a 


period of strong summertime zonal flow. The 
most suitable period found was interrupted by 
the development of a stagnant cold-vortex 
aloft over western Europe. Accordingly, the 
mean contour and anomaly patterns were 
computed over a split interval, 1—10 and 
26—30 June 1950 inclusive. The method of 
computation employed and the technique of 
representation used in the figures is identical 
with that described previously in connection 
with figures 8 and 9. 

The mean contour pattern aloft for the 
summer period shows a weak ridge at about 
10° E longitude which extends as far north as 
the North Cape region; as a result, the air- 
flow over the British Isles and Scandinavia has 
a significant southerly component. The pattern 
is also displaced slightly northward in compari- 
son with the winter case just discussed. As 
before one principal storm track and a pre- 
cipitation anomaly pattern essentially similar 
to the winter case in its general features, but 
displaced slightly northward, are observed. 
Only in the Mediterranean area, where the 
effects of the ridge just mentioned are more 
pronounced, is the isanomal pattern signi- 
ficantly changed; these changes are not dis- 
cussed as they lie outside the area of principal 
interest of this paper. Mean surface tempera= 
tures over the European area were generally 
high, reaching values of 6° C above normal 
over parts of Great Britain, the Swiss Alps 
and eastern Spain. Only in the mountains and 
along the western coast of Norway, and over 
the Ural mountains of central Russia were 
negative anomalies of 2° to 4° C below normal 
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Fig. 16. Isanomals of mean surface-temperature over the Atlantic-European area during a wintertime case of 
strong zonal flow aloft (1—31 January 1949). See figure 12 for an explanation of the method of representation used. 


recorded. It is believed that the strong positive 
anomalies observed over the British Isles 
and western Europe are largely the result of 
southerly flow components introduced over 
this area by the ridge before mentioned. In a 
similar way the negative anomalies appearing 
over Russia appear to be associated with a 
weak trough observed over the Urals. In the 
writer's opinion the isanomal pattern found 
over the Scandinavian peninsula is more typical 
of conditions produced by the “ideal” sum- 


mertime zonal flow pattern; i.e. negative 
anomalies which extend inland from the west- 
ern coast lines and progressively weaken 
toward the east with foehn effects producing 
moderate to strong positive anomalies in the 
lee of mountain complexes (note the Baltic 
region in figure 18). 

We may then summarize the regional 
climatic characteristics of a period dominated 
by strong zonal flow aloft over Europe as 
follows: 
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Fig. 17. Isanomals of precipitation over the Atlantic-European area during a summertime case of strong zonal 
flow aloft (1—10 and 26—30 June 1950). Sce figure 11 for an explanation of the method of representation used. 


1) Precipitation amounts over Scotland, Scan- 


dinavia and adjacent regions reach values 


several times normal (maximum 250 to 


300 % of normal). A secondary area of 


weak positive anomalies is produced in the 
mountains of southern Europe and along 
the North African coast. 

Abnormally high mean temperatures (2° to 
4° C above normal) are observed over the 
Continent, Scandinavia and the British 
Isles during wintertime zonal periods. The 


isanomal patterns associated with zonal 
flow during the summer are more complex 
but indicate slightly below normal tempera- 
tures (1° to 2°C below normal) over the 
coastal slopes of Norway and western 
Europe with areas of positive anomalies 
(2° to 4°C above normal) in the lee of 
mountain complexes. 

The major portion of all precipitation 
recorded is produced by frontal action 
(cyclonic activity) in combination with 
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Fig. 18. Isanomals of mean surface-temperature over the Atlantic-European area during a summertime case of strong 
zonal flow aloft (1—10 and 26—30 June 1950). See figure 12 for an explanation of the method of representation used. 


orography and falls over the Scandinavian 
area which lies in the major cyclone track. 
4) Strong westerly surface winds and a stormy 
weather regime, associated with the passage 
of intense cyclonic systems, prevail over 
northern Europe and Scandinavia with 
light winds and settled weather generally 
observed over southern Europe. 
We immediately recognize pronounced 
differences between the European climatic in- 
fluences of the blocked and zonal flow patterns 


just discussed. Considering first precipitation 
characteristics, we note that the zonal flow 
pattern aloft is associated generally with an 
increased total fall and an increased area of 
positive precipitation anomalies as compared 
with the blocked case (compare figures 15 
and 17 with 8 and 11). Furthermore, the posi- 
tion of the primary positive anomaly in the 
case of zonal flow is principally determined by 
the latitude of the westerly jet stream aloft, 
whereas in the blocked case, so long as pro- 
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nounced displacements of the contour pattern 
do not occur, the precipitation isanomal 
pattern is essentially fixed in position by geo- 
graphical features. This results from the fact 
that in the former case frontal precipitation 
produced by storm centers moving along the 
polar front are of prime importance, whereas 
in the second case orographic- and/or coastal- 
type precipitation dominates. The area of 
greatest contrast between the zonal and blocked 
types is found over the Scandinavian peninsula 
and northern British Isles and over the adjacent 
seaward areas to the northwest with a second 
area, of lesser importance from our point of 
view, situated along the North African coast. 

Comparing the temperature isanomal pat- 
terns we note that in blocked-type cases its 
position is fixed primarily by the location of 
the blocking-ridge with relatively strong 
positive anomalies to the west and north of 
the ridge and with relatively strong negative 
anomalies to the south and east. In contrast 
zonal flow aloft is generally associated with 
much weaker anomalies over the entire area, 
with the possible exception of central Russia. 
The position of the temperature isanomal 
pattern for the zonal case is primarily fixed by 
geographical features and is therefore quite 
stationary except that during the summer the 
pattern is quite sensitive to small changes in 
the direction of the westerly current aloft; 
i.e. variations in the magnitude of north— 
south flow components produce large distor- 
tions in the “idealized” temperature anomaly 
pattern. Two arcas of maximum contrast 
between zonal and blocked periods are found. 
One area is located over Scandinavia, Scot- 
land and the North Atlantic, and the other 
is located over southern Europe and the Bal- 
kans; in each area a reversal of the sign of the 
anomalies produced, occurs in shifting from 
one circulation type to the other. 

Particularly in northern Europe and Scan- 
dinavia a clearly marked difference between 
the surface weather regimes associated with 
these two circulation types is observed. On 
the one hand (blocked-type) generally moder- 
ate northeasterly winds, extensive convective 
cloudiness and shower activity, and a lack of 
storminess are observed; on the other (zonal- 
type) generally strong westerly winds and 
intense “cyclonic” (stormy) weather are ob- 
served. 


DANTE ERB 


= 


« In the foregoing paragraphs it has not been ° 
intended to imply that characteristic and 
strong climatic anomalies are exclusively asso- 
ciated with the two circulation types discussed. 
Variations of the blocked and zonal patterns, 
as well as other forms of an essentially distinct 
character, occur and collectively determine 
the observed circulation regime over any 
period of moriths or years. With respect to 
blocked patterns however, the zonal type 
may be said to be representative, at least in a 
dynamical sense, of the opposite extreme. 
Accordingly the description presented in this 
section has been included primarily to permit 
a comparison of the climatic influences of 
extreme circulation types. 


The Effect of Blocking Action upon 
Glaciation 


It may be of some interest to briefly discuss 
certain implications of the foregoing results 
with regard to present-day glacial regimes. 
The remarks to follow refer to glaciation in 
the Scandinavian area, although they apply 
with some modification to the Scottish and 
Alpine regions. If we may assume that the 
results reported by AHLMANN (1941) and those 
obtained by WALLÉN (1949) on the Kärsa 
Glacier are applicable to the Scandinavian area 
in general, mean summertime air temperature 
becomes the most important single factor 
influencing ablation, while the total winter 
snowfall more obviously controls accumula- 
tion. The results of the third section therefore 
suggest that blocking action is, in the mean, 
unfavorable to glacial growth both summer and 
winter. During the summer, above normal 
surface temperatures and a tendency toward 
clear or only partly cloudy skies are associated 
with blocking over the Scandinavian area (see 
figures 12 and 13). These conditions are ampli- 
fied by the tendency for summertime blocks 
to be initiated east of the yearly normal (see 
figure 2) and for westerly movements during 
summertime development to disappear (see 
Table Il). During the winter much below 
normal amounts of precipitation, occurring 
together with much above normal tempera- 
tures (see figures 8, 13 and 14), are detrimental 
to glacial growth. 

In contrast, the meteorological conditions 
over Scandinavia associated both winter and 
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summer with a zonal-type flow aloft appear, 
in the mean, to favor glacial growth. In winter, 
particularly over the western slopes of the 
Scandinavian mountain divide, heavy snow- 
fall, with surface temperatures only slightly 
above normal (but well below freezing), favors 
large accumulation. In summer, the presence 
of a cool, moist air mass, with mean air 
temperatures below normal (over western 
Scandinavia) and overcast skies, favors low 
values of ablation generally. Reasonably ac- 
curate estimates of the total accumulation and 
ablation on Scandinavian glaciation are avail- 
able only for the Kärsa (WALLÉN 1949) and 
Stor (AHLMANN 1949) Glaciers and for a four- 
year period. A comparison between these 
values and observed blocking activity supports 
the hypothesis suggested in the preceding 
paragraphs. 

In concluding this section we may speculate 
upon the röle played by blocked and zonal 
circulation types in the development of a 
glacial period. Evidence at hand (FLINT 1947 
and AHLMANN 1948) suggests that the great 
ice sheet, which partially covered the con- 
tinental land mass of Europe during the most 
recent (Würm or Quartenary) ice-age, was 
initiated in the mountains of the Scandinavian 
peninsula. It seems reasonable to suppose that, 
barring any appreciable change in land-mass 
profiles, a future European glacial epoch would 
find its genesis in the source region just de- 
scribed. We should therefore expect conditions 
existing in the Scandinavian mountains to be 
indicative of a developing ice-age. 

As has been suggested, a circulation regime 
predominately zonal in type, with the westerly 
jet aloft located at about 55° N latitude, appears 
to be a circulation pattern favourable both 
winter and summer for glacier growth in 
Scandinavia. The relatively low ablation and 
high accumulation values thereby favored 
would result in local glacier growth in the 
Scandinavian source region which if unchecked 
would eventually produce an ice cap. As the 
permanent ice cover begins its slow move- 
ment toward the south (gravitational flow), the 
size and intensity of the surface polar anti- 
cyclone will increase and, in turn, support an 
expanded circumpolar cyclonic vortex aloft. 
Such an expanded vortex in the upper tropo- 
sphere must, as it develops, cause the westerly 
jet stream aloft, together with its associated 
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surface cyclone track, to take up mean positions 
further and further toward the south. This 
will carry the area of maximum precipitation 
away from the original source region and 
nearer to the advancing edge of the ice sheet, 
thereby permitting glaciation to advance more 
rapidly and to reach more southerly latitudes 
than would be possible under the action. of 
gravity flow alone. This selfconditioning of 
the climatological environment of developing 
glaciation has been discussed in several in- 
stances by AHLMANN and THORARINSSON 
(1938). 

The southward displacement of the jet 
stream aloft, in association with the advancing 
ice-margin, must lead progressively to an 
intensification and concentration of the upper 
westerlies which, in accordance with the theory 
of block development advanced in Part I, will, 
in the mean, lead to more frequent and per- 
sistent blocking action. This progressive in- 
crease in blocking activity must eventually 
arrest and finally reverse the regime of the ice- 
mass. It is significant to point out that for a 
long period after total ablation exceeds total 
accumulation, an advance of the ice-mass may 
be expected under the action of gravitational 
flow alone. Accordingly, the establishment 
of a predominantly blocked-type circulation 
pattern over Europe, although in itself un- 
favorable to glaciation, would precede the 
period of maximum extension of the ice- 
margin by several millenniums. 

In a recent paper VIETE (1950) has de- 
scribed the general atmospheric circulation 
during glacial periods as reconstructed from 
geomorphological, geological, botanical and 
meteorological evidence. If one may interpret 
the phrase “increased general circulation” as 
principally denoting an increased exchange 
of air between high and low latitudes, the 
circulation regimes which he has associated 
with initial and developing glaciation are 
in general agreement with the foregoing 
outline. 

‘The importance of time-lag relationships in 
any discussion of the behavior of an ice-mass 
of such colossal proportions as those of the 
Wiirm Inland Ice Sheet cannot be over- 
emphasized. In particular and as shown 
above, we may not imply that the meteoro- 
logical regime associated with a period of 
maximum extension is, in itself, necessarily 
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favourable to the initiation and growth of an 
ice-mass. Although this fact has been previ- 
ously pointed out by several writers, it has 
frequently been overlooked in current litera- 
ture dealing with the meteorology of glacial 
epochs. A certain amount of confusion inevit- 
ably results. It is also interesting to compare the 
sequence of changes in general atmospheric 
circulation, which were suggested in the pre- 
ceding two paragraphs as being associated with 
the development and eventual decay of a 
glacial epoch, with those accompanying the 
“index cycle” (see Rosssy and WILLETT 1948). 
These two cycles, one of extremely-long and 
the other of relatively-short period, closely 
resemble each other in main outline as has 
been suggested by WILLETT (1949) 


ne 


Conclusion 


It may be desirable to summarize the 
principal results of this study in the form of 
the following conclusions: 

1) Blocking action is most frequently initiated 
in two relatively narrow longitudinal zones 
in the northern hemisphere, one (Atlantic) 
centered aloft at 10° W longitude and the 
other (Pacific) at 150° W longitude. 

2) Both Atlantic and Pacific blocks, once 
formed, normally persist for a period of 
12 to 16 days and are relatively stable in 
position. 

3) Regular seasonal and irregular yearly varia- 
tions in both Atlantic and Pacific blocking 
activity occur. These variations are strong, 
in excess of 50 %, and in the case of 
the seasonal trends are essentially in phase 
in both Atlantic and Pacific. 

4) Atlantic blocks exceed Pacific blocks 
in frequency of occurrence, apparently by 
a factor of two to one. 

5) European blocking action produces a char- 
acteristic and pronounced distortion of the 
normal surface-temperature and precipita- 
tion distributions over the underlying area. 
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6) Yearly variations in blocking activity over 
the European area produce recognizable 
fluctuations in the associated climatic trends 
(precipitation and mean surface-tempera- 
ture). 

7) Both winter and summer, European block- 
ing action produces climatic conditions 
which, in the mean, appear to be un- 
favourable to Scandinavian glaciation. The 
climate associated with a zonal flow pattern 
aloft appears, in the mean, to favor the 
growth of Scandinavian glaciers. 

8) It is suggested that the circulation regime 
associated with a developing ice-age should 
display the following salient characteristics. 
During the initial development of local 
mountain glaciation and the subsequent 
active development of an ice cap a pre- 
dominantly zonal-type circulation aloft 
prevails over the area in question. During 
a later stage in the development of the ice 
mass, but near the time of most rapid 
advance of the inland ice margins, blocking 
activity in the upper westerlies progres- 
sively increases until active growth of the 
ice mass ceases. Forward movement of the 
ice margin continues under the action of 
gravitational flow until, at the time of 
maximum extension, a predominantly 
blocked-type circulation aloft prevails. 
After this time a rapid retreat of the ice 
margin commences in association with a 
progressive, but lagging, decrease in block- 
ing activity in the upper westerlies. 
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Abstract 


The results of three high-altitude ascents of an experimental radiosonde incorporating an elec- 
tronic dew point hygrometer and carried by a large plastic balloon are presented. Significant 
features of the moisture distribution in the stratosphere revealed by these flights include at least 
one saturated layer in the vicinity of the tropopause on each occasion, and marked increases of 
the mixing ratio with height in the layer between 16 and 30 km in the two summer soundings. 
Charts and tabulations of the data from the three soundings are presented. 


I. Introduction 


In a short previous publication (BARRETT et 
al. 1949) the authors presented the results 
obtained from one flight of an experimental 
radiosonde to a height of 30 km, the out- 
standing feature of this sounding being a 
detailed picture of the distribution of water 
vapor as measured by an automatic, continu- 
ously operating electronic dew point hygro- 
meter. Since the date of writing of that paper, 
two additional flights to comparable altitudes 
have been carried out successfully, yielding 
new and valuable data on the distribution of 
water vapor at altitudes beyond the reach of 
aircraft and far above the levels at which any 
standard radiosonde humidity element yields 
meaningful results. 


1 Work performed under Contract N6ori-20, Task 
Order I, with the Office of Naval Research, U. S. Navy. 

2 Presently associated with Bendix Aviation Corp., 
Detroit, Michigan. 


2. Instrumentation 


Because of the fact that a complete discussion 
of the principles and operation of the sounding 
instrument requires a rather lengthy presenta- 
tion, only a brief outline will be given here. 
The complete discussion will appear in a paper 
by Suomi and Barrett to be published at a 
later date. Although the discussion of the 
instrument should logically precede the pre- 
sentation of the results of its use, it is felt that 
the data from the flights is of such interest to 
workers in meteorology and other fields that 
it should receive priority in publication. 

The complete sounding instrument consists 
of (1) an electronic dew point hygrometer, (2) 
a thermojunction for free-air temperature 
measurement, (3) a hypsometer for pressure 
measurement, (4) an ice-water bath for the 
reference thermojunction used in all three 
measurements, (5) a stable DC amplifier to 
raise the voltages developed by the thermo- 


THE DISTRIBUTION OF WATER VAPOR IN THE STRATOSPHERE 


couples to a level suitable for telemetering, 
(6) a standard cell for internal standardiza- 
tion of the amplifier and telemetering system, 
(7) a radio telemetering unit which establishes 
a unique correspondence between the DC 
output of the thermocouple voltage amplifier 
and the repetition frequency of pulses modu- 
lated upon a radio-frequency carrier wave, 
(8) a sequencing switch which causes the 
various telemetered quantities to be transmitted 
in turn, and (9) the batteries required for power 
supply. The batteries and all components 
except the sensing elements and the radio 
telemetering unit are contained in a pressurized 
and thermaily insulated chamber. 

The choice of transducers for the three 
meteorological variables is made from the 
standpoint of maximum accuracy at very high 
altitudes. The dew point method of hygro- 
metry is the only procedure in existence at the 
present time which is capable of reasonably 
uniform accuracy over the entire range of 
atmospheric temperature, pressure and mois- 
ture content. The thermojunction for air 
temperature measurement is so constructed as 
to have negligible thermal capacity and high 
albedo, thus minimizing errors due to lag and 
solar radiation. The hypsometer, with an 
approximately logarithmic relation between 
pressure and boiling temperature of the fluid 
used, tends to have a nearly constant per- 
centage error in pressure regardless of the 
absolute value of the pressure. 

From the practical viewpoint, this selection 
of transducers makes possible the conversion of 
the moisture and pressure measurements to 
equivalent temperature measurements so that 
thermojunctions can be used to make all three 
determinations. This greatly simplifies the 

- problem of calibration, standardization and 
telemetering. The temperature of the dew or 
frost condensed on the metal surface of the 
hygrometer mirror and the temperature of the 
boiling hypsometric fluid are measured by 
thermojunctions of the same type used for air 
temperature measurement, and the same 
amplifier and telemetering system is utilizable 
on a time-sharing basis for all three measure- 
ments. 

An additional simplification results from the 
choice of the hypsometer as a pressure trans- 
ducer. Since the dew point hygrometer 
requires a refrigerant to cool the metal surface 
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to the condensation point, it is possible to select 
a pure chemical substance with a low boiling 
temperature and a known vapor pressure curve 
to fulfill simultaneously the functions of re- 
frigerant and hypsometric fluid. 

The complete sounding instrument is admit- 
tedly quite large and complex as compared to 
the ordinary radiosonde, although it is prob- 
able that modern techniques of miniaturiza- 
tion could effect a considerable saving in 
size and weight. The device as flown has a 
weight of 46 pounds. This requires that a 
large balloon or cluster of smaller balloons be 
used for a flight. The 70-foot plastic sounding 
balloon manufactured by General Mills, Inc., 
Aeronautical Research Laboratories for the 
U.S. Navy was used as a vehicle for the 
flights reported herein. This balloon is capable 
of reaching altitudes of 103,000 feet or over 
with the 46-pound instrument plus a parachute 
and radiolocation transmitter as a load. Un- 
fortunately, the large size of this balloon makes 
ground handling and launching impossible un- 
der windy conditions, so that flights must be 
launched in the early morning hours or when 
the gradient wind is extremely weak. 

Ground equipment for reception and re- 
cording of the telemetered data consists of the 
SCR-658 rawinsonde receiver, electronic fre- 
quency meter and strip-chart recorder, as used 
for radiosonde reception, without modifica- 
tion of any kind. 

Because of the high stability of calibration 
of the thermojunctions and the completely 
negligible variation of gain of the thermo- 
couple voltage amplifier, all of the significant 
errors of measurement are concentrated in 
the telemetering and recording apparatus. A 
more complete discussion of these errors will 
appear in the later paper. Because of the internal 
standardization feature, these errors involve 
only small residual nonlinearities in the tele- 
metering and recording systems. Although no 
extensive checks have been made, it appears 
that the relative accuracy of the temperature 
measurements is represented by a standard 
error of about 0.3° C or less. This is only 
serious in the case of the hypsometer; it 
corresponds to an error of about 15 mb at sea 
level. Comparison of the pressure-height data 
from the three flights with synoptic radio- 
sonde data obtained from surrounding stations 
shows, however, that no great systematic 
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errors are introduced by this uncertainty in 
the pressure. 


3. Flight Procedure 


Because of the fact that the handling and 
launching of the large balloons requires an 
experienced crew and specialized equipment, 
all flights were performed in cooperation with 
Project SkyHook, the U.S. Navy project 
concerned with the testing of the performance 
of the General Mills balloons. Two of the 
flights to be described were made from the 
regular testing center at Camp Ripley, Min- 
nesota; the third was a special demonstration 
flight made at the St. Louis, Missouri, Naval 
Air Station. 

The restriction imposed by surface wind 
conditions described above, and an additional 
restriction by the Civil Aeronautics Authority 
on the maximum permissible cloudiness 
during a flight, made it necessary to confine 
flights to essentially anticyclonic situations. All 
three flights to be discussed were made in such 
situations from one to three days after passage 
of a cold front..For this reason the synoptic 
pattern was rather similar on all three flight 
occasions. It is to be regretted that at least one 
sounding in a deep cyclone could not have 
been made so that the differences between 
stratospheric conditions in typical cyclones and 
anticyclones could have been studied. 

After arrival at the launching base, the in- 
strument was calibrated, checked and placed 
in a stand-by condition. Synoptic charts and 
Weather Bureau forecasts were studied each 
day to determine the feasibility of launching 
on the following morning. If a favorable 
situation appeared imminent, all personnel 
were alerted in the evening and assembled by 
0500 hours on the following morning. If 
examination of teletype data indicated con- 
tinued absence of cloudiness and the local 
surface wind was low, the instrument was 
taken to the launching site, given a final test, 
loaded with the ice-bath and refrigerant and 
placed in operation. The load line was laid 
out and the balloon inflated during the same 
period. Launching was usually completed by 
0630. The balloons ascended at a rate of 700 to 
1000 feet per minute and remained at peak 
altitude for two to four hours, after which 
time the instrument and parachute were 
released by a clockwork mechanism. At the 
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time of release of the load an aircraft was in 
the air under the balloon position and the pilot 
kept the parachute under observation during 
the descent. A truck was dispatched to the 
region of expected landing and maintained . 
radio contact with the aircraft. In this way the 
load was usually recovered immediately after 


landing. 


4. Results of Flights 


The first successful sounding was launched 
at Camp Ripley at 0633 CST on ı July 1949. 
This flight has been discussed in the previous 
paper (BARRETT et al. 1949); however, the 
adiabatic charts are reprinted here in figs. ı 
and 2 for convenient comparison with the 
later flights. 
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centered near Bismarck, North Dakota. Cen- 
tral Minnesota was located in a region of 
very weak southwesterly gradient wind. The 
flow pattern at the 700-mb level was from 
west-southwest with little or no curvature of 
the streamlines. The main zonal current (jet 
stream) was located to the north of the Cana- 
dian border; thus the sounding was made in 
an air mass associated with the strong sub- 
tropical anticyclone which dominated most 
of the country at the time of the flight. 

The lower portion of the sounding (fig. 3) 
is characterized by a very steep lapse rate of 
temperature and a relatively smooth dew point 
trend. Relative humidities are rather low up 
to the 600-mb level and extremely low in the 
500- to 400-mb region. 

The dashed temperature trace below the 
820-mb level is an extrapolation required by 
the fact that the telemetering system did 
not provide for temperatures warmer than 
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It should be mentioned here that the term 
“dew point” as well as all numerical values of 
this quantity shown on the charts and in the 
tables, refers to the dew point for all values 
above 0°C and to the frost point for all values 
below 0°C. Any uncertainty with respect to 
the phase actually present on the hygrometer 
mirror is confined to the region between 0° C 
and -—8° C and is therefore of no great concern 
when attention is concentrated on the upper 
portion of the atmosphere. Tests made in the 
laboratory indicate that (at least in city-pol- 
luted air) the transition from water to ice 
takes place at —3° C on the average. 

The second flight, the results of which are 
shown in figs. 3 and 4, was sent aloft at 0617 
CST on 26 August 1949 from the Camp 
Ripley location. The surface synoptic pattern 
was characterized by a weak wave cyclone 


marae: 

Above the 400-mb level (fig. 4) the relative 
humidity begins to increase, reaching a value 
of so percent at the base of the complex 
tropopause at the 200-mb level. A sharp rise 
of the dew point at 150 mb in combination 
with a fall of the temperature produces a 
saturated layer (with respect to ice) at the base 
of what may be regarded as the real tropopause 
at 135—125 mb. Above this level the dew 
point falls gradually while the temperature 
first becomes isothermal at —60° C and then 
rises to a value of —39.5° C at the top of the 
sounding. Above the 45-mb level the dew 
point remains approximately constant near 


° 
23192 Ce 

No accurate wind data are available for this 
or the other flights since theodolite readings 
were taken only at five minute intervals; 
however an examination of the plotted balloon 
trajectories for both summer flights shows 
that the zonal wind component decreased 
with height above the 200-mb level and be- 
came easterly above the 70—80 mb region on 
both occasions. 

The third sounding, plotted in figs. 5 and 6, 
was made at NAS St. Louis at 1029 CST on 
7 January 1950. This ascent took place in the 
center of an intense polar anticyclone (central 
pressure 1040 mb, at sea level) which followed 
behind one of the most severe cold fronts of 
the season at St. Louis on 3 January. The. 
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700-mb pattern had been characterized by a ne mer TETE 
trough over the western U.S., giving south- if | SIN SESE Se ile 
westerly steering aloft and mild, foggy surface N ee N 
weather in Illinois and Missouri. On 2 January  * SI ha) FR 
this trough began to move eastward; this Ne ii a SJ 
development was accompanied by the growth en I den 
of a surface cyclone with a well-marked cold BEN! DE 


front which passed St. Louis on the afternoon 
of 3 January causing very heavy rainfall, high 
winds and tornadoes. The temperature fell 
from about 65° F to below freezing, and by 
the following evening had dropped to 8° F. 
Precipitation continued in the form of snow 
and freezing rain for most of the day following 
the frontal passage. 

By 7 January St. Louis was located to the 
west of the trough line and the flow at 700 mb 
was from west-northwest. The main zonal 
current was fairly broad and located mostly 
to the north of the station. Temperatures at 
the surface and aloft had moderated, skies had 
cleared completely and the surface wind had 
become light. The flight was launched at a 
later hour than usual in order to permit the 
audience (Naval reserve aerological officers) 
to assemble. Launching was successfully ac- 
complished in spite of slightly higher surface 
winds than were considered safe and in spite 
of the ice-covered ground which made running 
with the equipment difficult and hazardous. 

The principal point of interest in the lower 
portion of the flight is the extreme stratifica- 
tion of the moisture in the “cold dome’ below 
the weak frontal zone at 750—760 mb as 
contrasted with the region above. It is also 
interesting to note the very close agreement 
among the three flights on the dew point in 
the 400-mb region in spite of the difference 
in season. The persistence of values near — 50° C 
leads one to believe that this value might 
almost be regarded as a definite characteristic 
temperature in anticyclonic situations. 

The upper portion of the soundings (fig. 6) 
has many points of similarity to that of the 
26 August sounding in spite of the seasonal 
difference. As before, the relative humidity 
rises, on the average, above the 400-mb level 
and becomes high in the region of the base of 
the very complex tropopause region (220 mb). 
This high relative humidity is maintained, in 
spite of a considerable drop in temperature, 
until a sharp rise of the dew point produces 
the saturated layer at 115—130 mb. 
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Fig. 6. Flight of 7 January 1950: Upper levels. 


Above the saturated layer the sounding 
differs greatly from the two summer ascents. 
The temperature remains almost constant 
while the dew point falls to the lowest value 
(— 79.5° C) recorded on any of the flights. 

It is to be regretted that weak and fading 
radio signals caused the loss of accurate data 
between 90 and 30 mb and above the ıs-mb 
level. The weak signals resulted from the 
fact that the winds were westerly at all levels; 
the elevation angle of the balloon was approxi- 
mately 6 degrees at the conclusion of the 
flight. 


5. Analysis and Conclusions 


The distributions of temperature and pressure 
for the three soundings are in reasonable 
accord with expectations. The summer stra- 
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tospheric temperatures show an increase with 
height, while the winter temperature remains 
isothermal. As mentioned above, the zonal 
wind in the stratosphere was easterly on the 
two summer occasions and westerly for the 
winter case. The observed summer tempera- 
ture distributions obtained at Camp Ripley 
are in very good accord with the stratospheric 


Table 1. Flight Data for 1 July 1949 
eh lie Tem- Dew Vapor Mixing 
ea = perature| Point | Pressure | Ratio 
des. © des G mb gm/kg 
0.35 | 987 2 FTO) 17328 8.83 
0.58 961 27.52 ——- 5.6 | 3.82 PISE 
0.95 920: 16.0 —3-65| 45 3.05 
1.60 852 13.3 —4.3 | 4.2 372 
22S 788 10.7 — 6.9 | 3.38 2.67 
2.70 746 7.0 | — 12.5 | 2.10 AS 
3-17 | 704 9. | ISA 1.57 1.38 
3.46 680 3.0 | — 13.2 1.95 1.78 
3.83 650 1.8 | —13.5 | 1.88 1.80 
3.97 638 1.2 | — 16.2 | 1.48 1.44 
4-32 610 — 2.6 | — 16.3 | 1.48 RS 
4.94 564 — 6.7 | — 24.6 | 0.653 0.722 
5.58 519 — 10.4 | — 25.8 | 0.585 0.703 
GE 479 — 13.2 | — 39.2 | 0.142 0.184 
7.00 432 — 19.1 | — 40.9 | 0.116 0.167 
7.42 | 407 — 21.3 | — 50.2 | 0.036 0.059 
7.89 381 — 23.5 | — 46.8 | 0.054 0.088 
8.32 | 359 a fed 55:4: 10.0:089) 0.033 
8.46 353 — 27.4 | — 46.0 | 0.064 0.113 
8.SI 349 — 28.3 | — 52.0 | 0.031 0.054 
8.90 331 — 31.4 | — $1.2 | 0.033 0.062 
9.00 297 — 31.8 | — 46.4 | 0.061 O.116 
9.58 301 — 36.8 | — 48.2 | 0.048 0.098 
TOUTE 276 — 42.0 | — 55.8 | 0.018 0.042 
11.10 238 — 45.5 | — 64.4 | 0.0058 | 0.015 
12.07 206 — 50.5 | — 67.4 | 0.0038 | 0.0II 
12.70 187 — $$.1 | —.54.8 | 0.021 0.069 
13.82 157 — 62.4 | — 74.2 | 0.0013 | 0.0052 
14.00 152 — 60.5 | — 77.0 | 0.00084| 0.0034 
14.90 132 — 55.3 | — 64.0 | 0.0062 | 0.029 
16.04 TE — 53.7 | — 70.4 | 0.0024 | 0.014 
16.28 107 — 49.2 | — 48.0 | 0.050 0.293 
16.51 103 — 49.5 | — 53.2 | 0.027 0.164 
16.60 IOI — 49.1 | — 70.0 | 0.0026 | 0.016 
17.00 96.0 | — 48:5 | — 75-3 | 0.00II | 0.0073 
17.78 85.0 — 47.5 | — 76.4 | 0.00096] 0.0070 
18.30 78.0 — 46.9 | — 73.2 | 0.0015 | 0.012 
18.84 72.4 | —44.5 | — 73.0 | 0.0016 | 0.014 
19.75 65.0 | — 45.0 | — 73.3 | 0.0015 | 0.015 
20.45 57.2 | — 44.4 | — 73.7 | 0.0013 | 0.014 
21.04 52.0 | — 41.0 | — 67.9 | 0.0035 | 0.042 
22.00 45.4 | — 40.0 | — 70.0 | 0.0026 | 0.035 
23.07 38.6 | — 38.8 | — 65.2 | 0.0055 | 0.088 
24.21 32.8 — 38.5 | — 64.5 | 0.0058 | 0.109 
26.28 24.4 — 38.3 | — 68.7 | 0.0031 | 0.078 
30.07 14.0 — 38.2 | — 68.0 | 0.0035 | 0.154 
1 Est 
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data obtained at Bermuda by Crary (1950) 
through sound-propagation studies; no quanti- 
tative wind data are available for comparison 
with the winds computed from the acoustical 
data. The winter sounding at St. Louis appears 
very similar to the Alaskan winter data of 
Crary if one considers the difference in lati- 
tudes of the two locations. 

The moisture data obtained from these three 
flights represent wholly new information for 
the region above 12 km altitude, and a great 
improvement in resolution and precision at 
all levels as compared to usual radiosonde in- 
formation. The only previous direct informa- 
tion concerning the vertical distribution of 
water vapor at high elevations which is known 
to the authors is that obtained by Dosson 
et al. (1946) using a manually-controlled con- 
densation hygrometer installed in a high- 
altitude aircraft. The data reported in the 
present paper are in reasonable agreement 
with Dobson’s results up to 12 km if one con- 
siders the geographical separation of the loca- 
tions at which each set of soundings was made. 

Since moisture data for the region above 12 
km have not been previously available, the 
results for the three soundings are presented 
in numerical form in Tables I, I and II. Water 
vapor concentration has been expressed both 
as vapor pressure and mixing ratio. Such in- 
formation should be of considerable value to 
those concerned with problems of radiation 
equilibrium in the upper atmosphere. 

One of the most significant pieces of in- 
formation obtainable from the tabulated data 
is the fact that, for the two summer soundings, 
the mixing ratio increases markedly with 
height in the region from 16 to 30 km. In the 
case of the flight of 1 July the mixing ratio at 
30 km is more than twenty times as great as 
that at 17 km and the vapor pressure of water 
at the top of the flight is actually more than 
three times as great as that at 17 km. For the 
flight of 26 August the mixing ratio increase 
is more moderate; the value at 30 km is 
slightly more than seven times as great as 
that at 16.4 km. 

In the winter case the trend of the mixing 
ratio with altitude is quite the opposite, 
although the decrease is not as striking as the 
summer increase. The value at 28 km is 
slightly less than half that at 17 km. The fall 
in the mixing ratio appears to be concentrated 
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Table 2. Flight Data for 26 August 1949 
| h Tem- Dew | Vapor | Mixing 
Height | Pressure perature] Point |Pressure| Ratio 
2 zu des, @|des., © mb gm/kg 
OS 17.977 | 20.0 75-2, 1117226 11.00 
1.06 900 2208 11.8 |13.84 9.58 
1.40 865 20.0! LAE || OLS 6.72 
1.84 822 15.8 5-4 | 8.97 6.78 
PBs || GP 10.5 — $.0 | 4.02 3.24 
3208, 7x0 6.0 — 6.4 | 3.56 3.12 
3.52210.079 1.4 | — 11.6 | 2.25 2.09 
3 740 O54 — 0.8 | — 12.0 | 2.18 2.07 
4.51 $92 —7.§ | — 25.2 | 0.609 0.641 
4.73 | 576 — 9.0 | — 36.5 | 0.192 0.207 
4.93 SOI — 9.5 | — 29.5 | 0:405 0.449 
Swigh ||) 528 — 15.0 | — 37.0 | 0.180 | 0.214 
5.70 506 — 15.4 | — 39.3 | 0.147 0.180 
6.31 468 — 19.$ | —41.4 | O.III 0.147 
6.90 | 431 — 24.2 | — 47.0 | 0.058 0.085 
7.58 392 — 28.5 | — 52.8 | 0.028 0.044 
Sas ONO — 34.7 | — 53.0 | 0.027 0.048 
8.68 334 — 38.0 | — 55.5 | 0.020 0.037 
9.20 | 310 — 42.3 | — 56.7 | 0.017 0.035 
10.38 260 — 48.9 | — 60.0 | 0.0II 0.026 
TOYA. || BAG — $1.5 | —57.5 | 0.015 0.038 
11.58 215 — 57.6 | — 65.2 | 0.0041 | 0.012 
12.06 200 — 58.6 | — 64.2 | 0.0062 | 0.019 
12.44 188 — 58.5 | — 62.5 | 0.0076 | 0.025 
13.09 170 — 56.4 | — 64.0 | 0.0062 | 0.023 
13.46 160 — 56.3. | — 63.4 | 0.0067 | 0.026 
13.76 TS?) — 57.6 | — 64.2 | 0.0063 | 0.025 
14.40 138 — 59.0 | — 59.0 | 0.012 0.056 
15.02 126 — 60.6 | — 60.6 | 0.009$ | 0.047 
16.25 102 — 60.2 | — 68.0 | 0.0035 | 0.021 
17.80 80.0 | — 60.5 | — 65.4 | 0.0050 | 0.039 
19.58 60.0 | — 58.8 | — 67.0 | 0.004I | 0.042 
20.48 52.0 — 56.5 | — 66.5 | 0.0044 | 0.052 
22.22 39.5 | — 52.3 | — 69.4 | 0.0028 | 0.044 
23.05 34.8 — §1.§ | — 71.7 | 0.0020 | 0.036 
24.86 26.6 — 50.3 | — 71.8 | 0.0019 ‚| 0.044 
25.22 25.2 | — 49.0 | — 72.8 | 0.0016 | 0.041 
25.90 22.9 — 46.8 | — 71.4 | 0.0024 | 0.065 
27.15 18.7 — 46.7 | — 72.3 | 0.0018 | 0.062 
29.40 13.4 | — 40.4 | — 72.1 | 0.0019 | 0.088 
29.88 T2rS — 39.6 | — 69.1 | 0.003I | 0.150 
32.09 | 9.0 | — 39.4 | —72.4 | 0.0018 | 0.123 
1 Est. 


in the layers above 24 km; the loss of data 
between 17 and 24 km may, however, have 
caused the loss of some significant variations. 

The presence of at least one saturated layer 
on every ascent, when considered together 
with the similar observations of Dosson and 
co-workers (1946), leads to the conclusion that 
such a layer is a semi-permanent feature of 
the atmosphere, at least in middle latitudes. 
This in turn suggests that the temperatures in 
the lower stratosphere may be regulated to a 
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Table 3. Flight Data for 7 January 1950 


5 Tem- Dew | Vapor | Mixing 
Height | Pressure perature} Point |Pressure| Ratio 
nn des. @ | des, € mb gm/kg 
0.17 | 1,015 0.0 | — Io.5 | 2.5I 1.54 
0.42 984 2.0 — 6.7 | 3.50 2.22 
0.60 962 2, — 7.5 | 2.06 138 
0.74 946 2.4 — 8.5 | 2.94 1.94 
0.90 928 || BE, I.09 
1.45 865 0.3 | — 13.0 | I.99 1.43 
1.58 849 — 0.9 | — 29.3 | 0.405 0.296 
2.00 808 — 3.6 | — 13.3 | 1.96 1.50 
2.28 780 —4.5 | — 24.8 | 0.647 0.516 
VAT, 762 — $.1 | — 24.0 | 0.701 0.572 
2.56 WSS — 4.2 | — 23.1 | 0.758 0.628 
3.49 668 — 8.6 | — 27.8 | 0.478 0.445 
4.01 627 — 10.3 | — 33.2 | 0.202 0.200 
4.35 599 — 12.5 | — 32.8 | 0.216 0.224 
Sa7 537 — 16.6 | — 43.6 | 0.086 0.099 
5.61 $05 — 18.7 | — 44.3 | 0.080 0.098 
6.31 463 — 19.8 | — 40.0 | 0.129 0.174 
6.38 457 — 20.8 | — 44.8 | 0.073 0.099 
6.62 442 — 23.4 | — 48.5 | 0.053 0.075 
7.00 419 — 25.0 | — 43.4 | 0.086 0.128 
702 400 — 28.1 | — 50.1 | 0.039 0.061 
8.16 355 — 34.0 | — 47.4 | 0.054 0.094 
8.42 341 — 35.7 | — 52.0 | 0.031 0.056 
8.69 329 — 36.1 | — 53.1 | 0.027 0.05I 
9.40 297 — 41.1 | — 47.3 | 0.057 0.120 
10.03 269 — 44.2 | — 52.8 | 0.027 0.063 
10.28 259 — 47.9 | — 50.1 | 0.039 0.095 
10.56 249 — 48.7 | — 55.9:| 0.018 0.046 
10.95 235 — 52.3 | — 55.5 | 0.023 0.062 
11.66 | 212 | — 54.2 | — 57.2 | 0.016 0.048 
12.08 197 — 54.I | — 55.2 | 0.021 0.067 
13,76 166 — 53.2 | — 54.5 | 0.022 0.088 
14.33 139 — 57.0 | — 59.I | 0.012 0.056 
14.70 131 — 57.6 | — 57.6 | 0.015 0.072 
15.49 116 — $7.0 | — 57.0 | 0.016 0.087 
17.20 88.5 | — 60.0 | — 63.6 | 0.0068 0.048 
24.07 29.3 | — 59.0 | — 69.5 | 0.0028 0.060 
24.88 25.8 | — 59.2 | — 75.7 | 0.0011 0.026 
26.48 20.0 | — 59.2 | — 76.3 | 0.000906] 0.030 
26.78 19.0 | — 59.2 | — 78.4 | 0.00068] 0.022 
27.48 16.9 | — 59.1 | — 77.5 | 0.00080] 0.029 
28.21 15.8 | —59.I | — 79-5 | 0.00059| 0:023 


great extent by the absorption of long-wave 
terrestrial radiation by the water vapor (and 
possibly carbon dioxide) in that region. This 
appears likely in view of the fact that the other 
principal absorber, ozone, is present only in 
small concentrations in the lower stratosphere 
and should not contribute heavily to the 
heating of the stratosphere below an altitude 
of about 20 kin. Even at the higher altitudes 
(20 to 30 km) it appears that water vapor 
contributes greatly to the heat balance; the 
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close parallelism between the lapse rates of 
temperature and mixing ratio for the two 
summer soundings gives strong evidence for 
this. 

Some further evidence in support of the 
frequent occurrence of saturated layers in the 
stratosphere is provided in a private com- 
munication received by the senior author from 
Fr. E. Gherzi, S. J., Director for Meteorology 
and Seismology of the Zi-Ka-Wei Observatory 
at Shanghai, who reports the presence of 
radar echoes returned from altitudes between 
to and 30 km at frequencies of 60 and 214 
mc/s on many occasions. The radar equipment 
employed peak power of too kilowatts and 
Yagi arrays aimed vertically. Echoes were 
stronger at the lower frequency. + 

À chart showing the variation of mixing 
ratio (on a logarithmic scale) with height for 
the three ascents is shown in fig. 7 as a supple- 


309 


ment to the tabulated data. From this plot it is 
seen that the mixing ratio in the lower 6 km 
of the atmosphere is higher in summer (as 
would be expected), that the values for both 
summer and winter are nearly identical near 
the 7-km level, that the summer values are 
lower in the 9- to 22-km layer and that this 
situation reverses above 22 km so that higher 
values prevail at 30 km in summer. It should 
be remembered that all three soundings were 
made under anticyclonic conditions; these 
distributions are probably not representative 
of other types of synoptic patterns. 

The data forcefully point out that the con- 
centration of water vapor in the middle 
stratosphere is not small enough to be ignored 
in studies of the atmospheric heat balance, as 
has been assumed in the past. It appears quite 
probable that the absorptive power of water 
vapor up to 30 km for long-wave radiation 
is considerably greater than that of ozone in 
the same region, in view of the greater per- 
centage of water vapor present. 

The fact that the water vapor pressure and 
the mixing ratio may be higher in the middle 
stratosphere than lower down in the at- 
mosphere is a matter of interest from the 
viewpoint of atmospheric circulation. It is, of 
course, impossible to formulate any rigorous 
hypotheses on the mechanisms which carry 
water vapor to great heights on the basis of 
only three isolated soundings; however a 
few tentative statements may be advanced. 
Possible mechanisms by which water vapor 
may arrive at a given region are: (a) vertical 
mixing through a combination of molecular 
and small-scale eddy diffusion; (b) large-scale 
vertical eddy mixing; the water vapor rising 
in cyclonic regions and diverging aloft; (c) 
quasi-horizontal advection along isentropic 
surfaces with lateral eddy mixing; and (d) 
vertical transport in convection cells. The 
distributions of water vapor observed in 
these soundings would seem to eliminate (a) 
at the outset; the mixing ratio would be con- 
stant or monotonically decreasing with height 
if this mechanism were valid. Mechanism (b) 
appears doubtful since it demands that con- 
vergence and lifting in cyclones should extend 
up to heights of over 20 km. Mechanism (c) 
appears reasonable but actually begs the ques- 
tion: since none of the isentropic surfaces 
in the 30-km region ever descends into the 
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source region of moisture in the lower levels 
one is still left with the need for an explanation 
of the means by which the water vapor 
ascends initially to the region of the isentropic 
surfaces involved. 

Mechanism (d) may provide the means by 
which the water vapor can reach the required 
isentropic surfaces. Byers and BRAHAM (1949) 
have published data showing that radar echoes 
from thunderstorm cells frequently reach alti- 
tudes of 15 km and have been observed to 
attain altitudes of nearly 20 km on occasion. 
The ice crystals which are thus lifted into the 
stratosphere provide a rich source of water 
vapor which may be distributed by lateral 
mixing. This mechanism, suggested to the 
authors by V. E. Suomi in a private communi- 
cation, appears reasonable but must be checked 
by rough calculations of the total mass of 
water required to produce the observed 
mixing ratios and a comparison of this figure 
with an estimate of the mass of water ‘carried 
to such heights by the total number of strong 
convection cells existing on the earth at a 
given time. 

By way of a final summary, a few qualita- 
tive correlations deducible from the flight 
data are presented. The summer stratosphere 
in anticyclonic situations is characterized by 
an increase of temperature with height, 
easterly zonal winds increasing with height 
and mixing ratio increasing with height. The 
winter situation is characterized by isothermal 
conditions, westerly zonal winds and mixing 
ratio decreasing with height. The present data, 
taken in conjunction with the wind data of 
Crary (1950) and others, indicates that in the 
20 to 30-km region a circulation pattern 
exists which is antisymmetric with respect to 
the equatorial plane; the summer hemisphere 
having a warm anticyclone centered at the 
pole and easterly zonal winds first increasing 
as the latitude decreases, then decreasing and 
becoming westerly in the winter hemisphere. 
The winter pole in this model is the seat of a 
cold cyclone. Such a general pattern, already 


proposed by SCHERHAG (1948) and called to 
the author’s attention by H. Rızur and C.-G. 
RossBy, appears to be in agreement with all 
of the existing data. 

Any attempt to examine the hypotheses 
offered in this paper, or to formulate others 


more satisfactory, can only be made after the 


EARL W. BARRETT, LEE R. HERNDON, Jr, AND HOWARD JE GARTER 


acquisition of much more data on the tempera- 
ture, wind and moisture distributions in the 
stratosphere at many locations in all latitudes. 
The greatest value of the experiments reported 
here is the demonstration of the feasibility of 
making measurements of the moisture distri- 
bution at stratospheric heights by means of 
automatic balloon-sounding instruments and 
radio telemetering. 
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Note added in proof 


Since the preparation of the original manu- 
script of this paper, the authors’ attention has 
been called by A. B. Meinel of Yerkes Ob- 
servatory, the University of Chicago, to a paper 
by D. R. Bates and Marcel Nicolet, 1950 (The 
photochemistry of atmospheric water vapor, J. 
Geophys. Res. 55, 301—327), in which these 
authors attempt to account for the intensity of 
the emission spectrum of excited hydroxyl (OH) 
radicles in the light of the night sky (first iden- 
tified by Meinel, 1950, Astrophysical Journal 111, 
555—564) by calculations involving the rates of 
a number of possible chemical reactions. 

In carrying out their calculations, these authors 
assume that the concentration of water mole- 
cules in the stratosphere is constant and equal 
to one part in 105 of the total number of 
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molecules present. This concentration corresponds 
to a mixing ratio of 6.2 X 10-3 gm/kg in the 
stratosphere. Comparison of this figure with the 
data of Tables 9.1, 9.2 and 9.3 indicates that this 
estimate is too low by nearly one order of 
magnitude in summer, and slightly too low in 
winter. Since Bates and Nicolet are just able to 
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account for the observed intensity of the Meinel 
bands on the basis of their assumptions, it ap- 
pears that their calculations could be placed on 
a more secure footing by an upward revision 
of the average mixing ratio in the stratosphere 
in accordance with the data given in the present 


paper. 
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Some Aspects of Argentine Rainfall 


By M. J. OLASCOAGA, Army Meteorological Service, Argentina : 


(Manuscript received 25 April, 1950) 


Abstract 


This report investigates at first what percentage of the number of days with rain accounts 
for most of the rainfall in widely separated climatic provinces of Argentina. It turns out that this 
percentage, which may be termed the relative precipitation intensity, is a constant in space 


and time. 


Then the frequency distribution of natural rainfall periods (rainstorms) is analysed in relation 
to the rainfall derived from these disturbances. It is found that the same percentage of the rain- 
storms accounts for the same large fraction of total rainfall regardless of season or geographical 


effects. 


It is well known that average annual or 
monthly precipitation values at any station 
fail to give some of the most important in- 
formation wanted from such data. In particular, 
frequency distributions of rainfall amount over 
any time interval seldom, if ever, follow the 
Gaussian law of normal distribution but are 
skewed. The number of days (or months) with 
rainfall below average exceeds those with 
rainfall above average. Large seasonal and 
regional differences of skewness, however, 
are to be expected. It is of practical importance 
to know these differences and to determine 
the number of days that account for the bulk 
of the rainfall at any location. 

This report gives the result of computations 
made with this objective for different climatic 
provinces of Argentina. In order to hold the 
statistical labor within feasible limits ten char- 
acteristic areas were selected (fig. 1). Within 
these areas, the rainfall was analyzed by means 
of the daily rainfall values given in Carta del 
Tiempo published by the Argentine Meteoro- 
logical Service. The study covers five years, 
1944—1948, a sufficient statistical sample for 
the present purposes. 

The selection of the ten areas shown in 
fig. 1 was partly dependent on the following 
consideration. Analysis of rainfall data taken 
at a particular rain-gage often is of restricted 
value only. Such data on many days may not 
be representative even for the surroundings 


of the gage, especially when the clouds are 
cumuliform. Integration over several closely 
adjoining gages is preferable. For this reason, 
it was attempted to locate groups of ten sta- 
tions situated near each other and climatically 
as uniform as possible, especially with regard 
to monthly and annual rainfall. This objective 
was readily realized in such provinces as 
Buenos Aires, but necessarily could not be 
satisfied fully in the mountainous and sparsely 
settled regions. In the far south, it was not 
possible to form a group. There, five. stations 
were treated separately.! 

Fig. 2 shows the mean annual rainfall of 
Argentina. It is readily seen that in spite of 
the difficulties in selecting material it was 
possible to locate nine zones consisting of 
10 stations each in climatic regimes that vary 
widely as concerns latitude, proximity to high 
mountains or coast, dynamic regime and total 


. precipitation. Table ı gives the average seasonal 


precipitation measured in each zone during 
the s-year interval, and Table 2 gives the 
corresponding number of days with rain. 


Rainfall analysis 


The analysis procedure was as follows: At 
first, the average precipitation was determined 
for each day of the s-year period in each 
zone. Then these daily values were written 


ı The names of all stations utilized in this report 
appear in an appendix at the end of the report. 
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Fig. 1. Principal topographic features of Argentina, and 
location of the zones used in study (mean elevation of 
each zone in m). 


in ascending order of amount on a monthly 
basis (for 5 Januaries, 5 Februaries, etc.) and 
both rainfall and number of days with rain 
were summed on a cumulative basis. Fin- 


Table 1. Mean seasonal and annual rainfall (mm) 

for 1944—48 in each zone and at the southern 

stations El Turbio (52° S, 72° W) and Ushuania 
(55° S, 68° W). 
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Fig. 2. Mean annual rainfall of Argentina (mm). 


ally, all cumulative subtotals were expressed in 
percent of the totals. Following is an illustra- 
tion of this computation scheme (not based on 
actual data). 


Table 2. Mean seasonal and annual number of 
days with rain, 1944—48, in each zone and at 
El Turbio and Ushuaia 


an.— | Apr.—| July— | Oct.— 
Zone ee ne Sept. | Dec. ace 
| 

Tucumän....| 350 76 II 194 631 
Monte Caseros| 319 335 250 341 | 1,245 
Santa Ferse eso 143 77 230 782 
Gordobar..n 226 66 36 213 541 
San Rafael... 2120 47 24 69 260 
Mendoza..... 89 13 24 48 174 
Mar del Plata..| 200 197 157 153 707 
Bahia Blanca.. 97 78 94 II4 383 
PuertoDeseado| 19 | 22 Ser 81 
El Turbio.... 33 43 so 35 161 
Ushualar. en. 142 127 100 123 492 


Jan.— | Apr.—| July— |Oct.— 

Zone Mar. ae sh Dec. Yea 
Mucumonr. ere 34 17 8 25 84 
Monte Caseros 33 26 26 24 109 
Santa Reis sae 31 21 16 29 97 
@oördobapan. 30 13 9 29 81 
San) Rafaeli... 22 8 7 17 54 
Mendoza..... 18 7 7 iS 47 
Mar del Plata.. 35 34 Fly 34 140 
Bahia Blanca. 25 17 17 31 90 
PuertoDeseado| 14 16 19 II 60 
El’ Wurbior... 8 7 7; 7 29 
Wish tala penser 36 30 26 31 123 


314 
4 — — 
Cum % days| No days RR Cum Cum 
with RR cum Free RR (mm) RR 
8 is ik I I 
17 2 I 2 2 
25 3 2 4 5 
33 4 3 7 8 
42 5 3 10 12 
so 6 4 14 17 
58 ws 6 20 2 
67 8 7 2 32 
75 9 9 36 43 
84 10 12 48 57 
92 II 16 64 76 
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Following completion of the computation 
sheets, the results were plotted on graphs with 
cumulative percent of the number of days with 
rain and cumulative percent of rainfall as 
axes (cf. figs. 4,11). From these graphs it is 
possible to determine what percentage of the 
total rain falls during what percentage of the 
number of days with rain. In order to study 
seasonal and regional variations, is was con- 
venient to investigate certain significant per- 
centages only. From the appearance of all 
curves it was clear that a study of the levels 
50 % and 75 % of the precipitation (measured 
downward from 100 %) would bring out the 
significant features of the distribution. Figs. 
3a, b show the seasonal trend of these Exo 
levels in two of the zones, which are widely 
separated geographically and situated in en- 
tirely different rainfall regimes. The graphs say 
that 50% of the rainfall is accounted for by 
less than 20 %, of the number of days with rain, 
and 75 % of the rainfall is accounted for by 
less than 30% of the number of days with 
rain. Neither zone has a regular scasonal trend 
and exhibits only minor irregular fluctuations. 
As figs. 3a, b are typical for all zones, it 
follows that the relative concentration of rainfall 
is independent of season in Argentina. 

Table 3 shows the annual means and mean 
deviations of the two precipitation levels for 
all zones. For this purpose, the results for the 
five stations of the southernmost zone could 
be combined, since they were practically 
identical. If turns out that the percentages are 
practically independent of geography, rainfall 
tagime and total rainfall in addition to season, 
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Fig. 3. Seasonal distribution of the percent number of 
days with rain (ordinate) that account for the indicated 
percentages of the monthly precipitation for two zones. 


especially for the 50% rainfall. The extreme 
rainfall regime of southern Patagonia and 
Tierra del Fuego proves the only ‘exception, 
and even there the departure is small at the 
50 % level. 

Fig. 4 gives the composite curve for all zones 
combined. This curve can be closely approxi- 
mated by the analytical expression. 


y = axe™, (1) 


Table 3. The percentages of the number of days 
with rain that account for 50% and 75% of the 
annual precipitation in each zone and in all zones 
combined. The mean deviations refer to the depar- 
tures of the monthly values (as shown in fig. 3) 
from the annual mean in each zone. 


| so % of total 75% of total 
precipitation precipitation 
Zone Mean Mean 
Mean | devia- | Mean | devia- 
tion tion 
HLUCUMAR SS er 12 3 25 6 
Monte Caseros 13 2 27 4 
Santa: Recon. 2 13 2 27 4 
Cordobe cc kon 13 2 28 4 
San Rafael... 1$ 3 32 4 
Mendoza...... IS 5 29 7 
Mar del Plata... 16) 2 25 4 
Bahia Blanca... t2 4 26 5 
Puerto Deseado. 12 5 27 S 
SUR Ever one 20 2 | 41 | 3 
Mean for all 
BONES AEs LE | 14 | 3 | 29 | 5 
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where y is the rainfall expressed in cumulative 
percent of the total, x the number of days with 
rain expressed in cumulative percent of the 
total, and a and b are constants (a= 2.30 : 1072, 
b = 3.78: 102). This relation has just been 
shown to hold over almost all of Argentina, 
from the subtropics to the subantarctic, from 
the Atlantic Coast to the Andes Mountains, 
in regions where the mean annual rainfall 
ranges from 80 mm to 1,250 mm and the 
annual number of days with rain from 
30 to 150 days. It is therefore plausible to sug- 
gest that equation (1) represents a general 
statistical law that is valid for most regions of the 
globe. At this time it is not possible to prove 
this contention because of the labor involved, 
but reference is made to the rainfall of San 
Juan, Puerto Rico, (RıekL, 1950) where a 
similar distribution prevails. The practical 
usefulness of the distribution law (1) for 
climatological and hydrological purposes — if 
proven generally valid — requires no elabora- 
tion. 


Rainstorm analysis 


It is of interest now to turn to some dy- 
namical features of Argentine rainfall. The 
measurement of precipitation in 24-hourly 
amounts is a statistical device that does not 
give any information about the frequency 
of natural rainfall periods (rainstorms). Re- 
cently, RIEHL (1949) made an application of 
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Fig. 4. Percent cumulative rainfall (abscissa) against 
percent cumulative number of days with rain (ordinate) 
for the annual precipitation of all zones combined. 
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Fig. 5. The percentage of the mean annual rainfall that 
is not due to organized rainstorms. 


the rainstorm concept developed in hydrome- 
teorology to the rainfall of Hawaii. The 
following treatment for Argentina will follow 
the method used for Hawaii. 

Rainfall can occur in organized form pro- 
duced by general atmospheric disturbances, 
or in form of scattered local precipitation. 
It is reasonable to expect that during passage 
of a synoptic disturbance more than half of 
the stations located in a zone will report some 
precipitation. Therefore it will be assumed 
that a rainstorm has occurred when more than 
five stations of a zone record rainfall on a 
given day. The duration of a rainstorm is the 
number of consecutive days with precipita- 
tion measured at the majority of the stations. 

Based on the foregoing definition, it is 
possible to calculate the rainfall due to rain- 
storms in each zone. The difference between 
this amount and total rainfall is the scattered 
local precipitation. Fig. 5 gives the percentage 
of annual rainfall that is derived from such local 
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Fig. 6. The percentage of the annual number of days with 
rain that is due to days with unorganized precipitation. 


precipitation. The percentages are small every- 
where and in no part of Argentina approach 
the amounts noted in Hawaii. Nevertheless, 
a pronounced gradient is evident, generally 
directed from the areas with high toward those 
with low annual rainfall. This situation corre- 
sponds to that found in Hawaii. 

If we consider the percentage of the number 
of days with rain that is due to days with 
local precipitation (fig. 6), we encounter much 
larger numbers, again in correspondence to 
Hawaii. The gradient is the same as in case of 
fig. 5. The areas with lowest total rainfall have 
the highest frequency of scattered rain, which 
usually is the least effective form of precipita- 
tion for accumulation of soil moisture. 

Fig. 7 shows the annual number of rain- 
storms, fig. 8 the annual number of rain- 
storm days and fig. 9 the average durations of 
rainstorms. It is seen that the storm frequency 
increases rapidly from the mountains toward 
the coast, but that the duration of the disturb- 
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Fig. 7. Mean annual rainstorm frequency. 


ances varies within relatively narrow limits, 
roughly from two to three days. Fig. toa 
shows the storm frequency as a function of 
longitude, therefore distance from the Andes 
Mountains, in cross-section form. Many of the 
provinces experience months without any 
rainstorms, especially in winter. The number 
of months without storms also is mainly 
dependent on proximity to the Andes (fig. 
to b). 

We can now analyze the rainstorms in the 
same way as the daily precipitation and ask 
what percentage of the rainfall attributable to 
disturbances is accounted for by what per- 
centage of the number of rainstorms. The re- 
sult of the analysis (table 4) is quite analogous 
to that previously encountered in table 3. 
There is a remarkable uniformity at both 50% 
and 75% levels of precipitation, and the mean 
deviations are only slightly larger than in case of 
the daily rainfall. The southernmost zone (Sur) 
is not included in table 4 since, as stated 
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Fig. 8. Annual number of days with rain due to rainstorms. 


initially, it was not possible to form a group of 
10 stations in that region. Because of the great 
uniformity, both in space and time, it is again 
possible to combine all zones and construct a 
graph representative for all of Argentina (fig. 11). 
The curve shown in this figure also can be 
approximated by the distribution law (x). In 
ths case the constants have the values a = 
=1.73- 1071 and b = 1.74 : 10 2. General va- 
lidity of these constants over a major portion of the 
globe is again suggested. 
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Fig. 10a. East-west profile of annual rainstorm frequency. 
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Fig. 9. Mean duration of rainstorms (days). 


Within Argentina, we obtain a clear picture 
of the meaning of fig. 11 through comparison 
with figs. 7—8. Since the upper 50% of the 
rainfall stem from slightly more than 20 % of 
the number of rainstorms, the region near the 
mouth of the Rio de la Plata will have about 
20—25 days with appreciable rainfall in the 
annual average and the Andes foothills s—6 
days. The same so % of the precipitation fall 
on eight to ten occasions per year in the 
extreme east, but are concentrated in 2—3 
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Fig. 10 b. East-west profile of the number of months 
per year that are without rainstorms, 
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Table 4. The percentages of the number of rain- 
storms that account for 50 % and 75 % of the annual 
precipitation due to rainstorms in each zone and 
. in all zones combined. Mean deviation as indi- 
cated in Table 3. Months wichout any rainstorms 
(cf. fig. 10 b) are not included in the deviations. 


| 50% of ‘total 75% of total 


rainstorm rainstorm 
Zone precipitation precipitation 

Mean Mean 

Mean | devia- | Mean | devia- 

tion tion 
Tucuman | 25 | 3 48 4 
Monte Caseros. 22 3 43 6 
Sata Fer ner 22 5 43 7 
Gordobar -- DB 5 45 4 
Sans Rafael. 21 6 45 7 
Mendoza...... 25 6 48 8 
Mar del Plata.. 22 3 42 4 
Bahia Blanca... 2 4 45 3 
Puerto Deseado. 22 YU Kr MR 

Mean for all | | | 
ZONES Gyo ain 23 4 45 5 | 


disturbances in the west. This is a precarious 
situation that demonstrates clearly how small 
changes of the general circulation from year to year 
(small displacements of mean troughs and ridges 
aloft) can produce sharp fluctuations in the annual 
rainfall. It has already been shown by Biel 
(for convenient reference, see Trewartha, 1943, 
p. 360) that a high annual variability actually 
exists. Biel’s isolines of equal rainfall variability 
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Fig. 11. Percent cumulative rainfall due to rainstorms 

(abscissa) against percent cumulative rainstorm frequency 

(ordinate) for the annual precipitation derived from rain- 
storms of all zones combined. 


closely parallel the lines of annual storm fre- 
quency given in fig. 7. 
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List of stations used in the computation of the rainfall 
indices. 


Tucuman: Bella Vista, Concepciön, Famailla, Lules, Mon- 
tero, Padilla, Pala Pala, Rio Colorado, Rio Seco, Tucumän. 


Monte Caseros: Acuña, Baibiene, Bonpland, Curuzü 
Cuatiä, La Cruz, Libertad, Mercedes; Monte Caseros, 
Paso de los Libres, Tapebicua. 

Santa Fé: Canada Rosquin, Coronda, Esperanza, Gal- 
vez, Las Rosas, Pilar, Santa F6, Totoras, Crespo, Parana. 

Cördoba: Alta Gracia, Cördoba, Cosquin, Jesüs 
Maria, La Calera, La Cumbre, La Falda, Oncativo, 
Pilar, Villa Rosario. 


San Rafael: Corral de Lorca, Cuadro Nacional, Gaspar 


Campos, Goico, Goudge, Media Luna, Monte Coman, 
Ovejeria, Pampa del Tigre, San Rafael. 


Mendoza: Alto Verde, Buena Nueva, Lujän de Cuyo 
Medrano, Mendoza, Palmira, Rivadavia, Santa Rosa 
Tres Portenas, Ugarteche. 

Mar del Plata: Ayacucho, Balcarce, Calfucura, Gral. 
Madariaga, Juancho, Las Armas, Maipü, Mar del Plata, 
Miramar, Vivoratä. 


Bahia Blanca: Algarrobo, Argerich, Bahia Blanca, 
Buratovich, Fortin Mercedes, N. Levalle, Medanos, 
Montes de Oca, Ombucta, Pedro Luro. 

Puerto Deseado: A. De Biedma, Cerro Blanco, Faro 
Cabo Blanco, Fitz Roy, Jaramillo, Pico Truncado, 
Puerto Deseado, Ramén Lista, Tehuelches, Tellier. 


Sur: El Turbio, Faro Cabo Virgenes, Rio Gallegos, 
Rio Grande, Ushuaia. 


Letters to 


Variations in the Earth’s Angular Velocity 
Resulting from Fluctuations in Atmospheric 
and Oceanic Circulation 


By F. H. van den DUNGEN, 
J. F. COX and J. van MIEGHEM 


University of Brussels 
November 6th, 1950 


1. It is generally accepted that available aerological 
information concerning the distribution, in space 
and time, of the zonal component of the atmospheric 
circulation, is not yet complete enough to permit a 
reliable calculation of the angular momentum of the 
atmosphere. 

If it is remembered that the zonal circulation 
varies strongly with longitude (Namıas and Crapp, 
1949), even probably more strongly in the northern 
than in the southern hemisphere, it will appear that 
a calculation based on the information available 
about the distribution of the zonal velocity in one 
or two meridional cross sections only, does not seem 
very satisfactory. 

W. H. Munk and R. L. Mitier (1950) have 
calculated the angular momentum of the southern 
atmosphere using winter and summer cross sections 
along the 150°E meridian, as established by F. 
LoEweE and U. RADOK (1950 a). These cross sections 
show that the wind speed in the center of the west- 
erly jet exceeds 120 mi/h (at 25° S), in winter, and 
50 mi/h (at 40° S) in summer. 

More complete and recent information show 
(Lorwe and Rapox, 19506) that the summer 
southern “jet stream” would in fact be much more 
intense, exceeding 100 mi/h in the center. The dis- 
similarity or assymmetry of the westerlies in the two 
hemispheres is, in summer, much more marked 
than previously assumed. 

We are of the opinion that the annual variation 
of the momentum transfer between the earth and the 
atmosphere offers a better basis for calculating the 
annual variation of the angular momentum of the 
atmosphere. 

This transfer results from the friction between the 
earth surface and the atmosphere (including the 
pressure effect due to mountains). 

Utilising tables calculated by V. P. Starr and his 
collaborators (1949) on the basis of numerous sur- 
face information, we believe to have followed a 
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sound method in our calculations (vAN DEN DUNGEN 
et al. 1950). r 

2. The flux of west momentum across the earth’s 
surface is an ascending one where easterlies dominate 
(low latitudes) and a descending one where wester- 
lies dominate (middle latitudes). In the mean, as 
easterlies are better developed in the northern than 
in the southern hemisphere, the ascending flux of 
the low latitudes is not compensated in the northern 
hemisphere by the descending flux of the middle 
latitudes. A part of the ascending flux of the tropical 
regions of the northern hemisphere passes across the 
equator (WIDGER, 1949) and is compensated in the 
middle latitudes of the southern hemisphere, where 
more intense westerlies bring a descending flux of 
west momentum, more considerable than in the 
northern hemisphere. 

Assuming that for the entire surface of the earth, 
at a uniform rate throughout the year, a complete 
compensation of the ascending and descending west 
momentum happens, we obtain (vAN DEN DUNGEN 
et al. 1950) a scheme of the fluctuations of the earth’s 
rotation in good accordance with the published 
values (Stoyko, Finch, Uhink). 

When numerical tables of the kind calculated by 
V. P. Starr and his collaborators and related to the 
southern hemisphere are available, the assumption 
about yearly compensation may be tested. It seems 
probable that a complete yearly compensation at 
uniform rate does not occur and this may open the 
way for a meteorological explanation of apparently 
erratic changes. In that connection it must be re- 
membered that the reduction of the clock rates 
remains intricate and that future work may pro- 
vide further valuable information in that field. 


3. Column (4) of Table 2 in W. H. Munk and 
R. L. Miller cited paper needs, in our opinion, to 
be carefully examined. Let Ap represent the varia- 
tion from January to July of the mean pressure p 
(on the ground) not reduced to sea level, along a 
parallel of latitude; obviously, the total mass of the 
atmosphere remaining constant, it 1s necessary that 
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In other words, the mean value of the A p in the 
northern hemisphere must be equal and opposite 
in sign to the mean value of the A p in the southern 
hemisphere. Dividing the surface of the spherical 
globe into zones of 5° and 10° latitude, let us 
denote by [A p] the mean values of A p for every of 
these zones. 

Equation (1) shows that the [4 p] brought into a 
(sing, [A p]) diagram must appear such (stairs 
curve) that the algebraic sum of the areas comprised 
between it and the “sin y”’ axis be zero. 

It is easy to check that the [A p] values tabulated 
by W. H. Munk and R. L. Miller do not fulfill 
this condition. As a special case, one finds with 
those values 


80° | + 0.85 mb for the northern 
Ap = : hemisphere ; 
Bas 170% | — 2.19 mb for the southern 


hemisphere; 
it follows that 


Ap = — 0.67 mb, 


0723 
representing the mean value of A p be- 
tween latitudes qm; and q@2. 


Ap 
4 Pr 


It seems quite unlikely that this mean negative 
pressure difference between 80° S and 80° N could 
be balanced by a reasonable positive difference in 
the two small polar caps of 10° of latitude. 

4. In their above quoted paper W. H. Munk and 
R.L. Miter (1950) report about our paper (VAN DEN 
DunGen et al. 1949) in terms such that it could be 
understood that we were not aware of the importance 
for the problem stated of the correction from sea 
level pressure to local pressure. 

We believe that every reader who will have read 
through our paper will be aware that this was not 
the case. Our scheme was fully developed by one 
of our students, Miss M. LAURENT (1949, 1950). 
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Reply 


By W. H. MUNK and R. L. MILLER 
Scripps Institution of Oceanography 
La Jolla, November ısth, 1950 


1. The effect of the wind circulation on the earth’s 
rotation can be computed either from (a) changes in 
the angular momentum of the atmosphere as a 
whole, or (b) momentum transfer between atmos- 
phere and earth. We have utilized the former 
method, and arrived at results consistent with 
astronomic observations; as described in the above 
note, VAN DEN Duncen et al. (1950): utilized the 
latter method, with similar results. 

It is true that for our method it would have been 
highly desirable to have had available additional 
meridional sections in the Southern Hemisphere, 
and we have said so in our article. In van den Dun- 
gen’s method, the pressure effect due to moun- 
tains and a choice for the frictional coefficient be- 
tween air and land, and air and water, is also not 
without uncertainties. Since the transfer of angular 
momentum over six months from the westerlies or 
the easterlies to the earth is large compared to the 


angular momentum of the atmosphere at any instant, 


{ This paper was apparently published at about the 
same time as our paper, but no copy has as yet been 
available to us. In our article we refer to van den 
Dungen’s earlier paper in which the variations in the 
earth’s rotation are ascribed to latitudinal changes of 
air masses. 
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van den Dungen’s method must depend on the 
small difference between two large and somewhat 
uncertain terms, whereas in our case the difference 
is of the same magnitude as the terms from which it 
is formed. 

2. The pressure data in Table 2 go back largely to 
Napier Shaw, the source used also by van den 
Dungen in his computations. The fact that the 
pressure changes do not integrate to zero over the 
whole globe can be ascribed to uncertainties in these 
observations, and in the conversion of sea level 
pressure to local pressure. We have accordingly 
stated that “it would be desirable to recompute 
Table 2 using the original data for the uncorrected 
local pressure”. In any event it can be shown that 
our estimate for the magnitude of the inertia term 
(which is small compared to the wind term) will 
hardly be effected by allowing for this discrepancy. 

3. We regret to have given the impression that van 
den Dungen et al were unaware of the importance 
of reverting the sea level pressures to local pressures. 
The complete corrections involve the elevation and 
the temperature of the weather stations. The correc- 
tions, as entered in Table 1 of van den Dungen’s 
paper, take into account only the elevations. This 
part of the correction, which is of course the same 
throughout the year, will drop out from computa- 
tions of the seasonal differences with which our prob- 
lem is concerned. This explains why van den 
Dungen’s February minus August values (their 
Table 4) are nearly the same as our uncorrected January 
minus July values (our Table 2, column 2). 

The complete corrections depend partly on fem- 
perature, and these corrections do not drop out when 
differences between seasons are computed. The 
corrections are remarkably large. North of 30° N 
latitude the average local pressure variation between 
January and July is about minus 1.5 mb, whereas 
van den Dungen et al. obtained roughly plus 4 mb. 


Comments on the Behavior of Jet Streams 
over Eastern North America 


Bye. Ba CLAPP and J. S. WINSTON 
U. S. Department of Commerce, Weather Bureau 


October ııth, 1950 


Mr Phillips’ article in the May 19 so issue of TEL- 
LUS (Vol. 2, No. 2) touches on an interesting point 
regarding confluence. It is quite true that the mere 
confluence of streamlines on upper level charts does 


not necessarily lead to increased thermal gradients 
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and consequent increases in high level wind speeds. 
Confluence of the streamflow must be accompanied 
by differential temperature advection at most levels 
of the troposphere before it can result in accelerating 
wind speeds in the upper troposphere. However, 
Mr Phillips fails to consider that lack of strong 
advection at 300 mb may not represent conditions 
at all levels. It is very likely that there is cold air 
advection in the northwesterly flow and warm air 
advection in the southwesterly flow at middle and 
lower levels of the troposphere (i.e. up to about 
400 mb). Hydrostatically, the mean temperature 
field of the layer below 300 mb must be similar in 
appearance to the 300 mb height field, but it is well- 
known that the flow at lower levels in middle lati- 
tudes is generally less zonal and more meridional 
than at 300 mb. Thus, even in the case Mr Phillips 
presents the streamline-isotherm relationship pre- 
scribed by the confluence theory of Namias and 
Crapp (J. Meteor., 1949, Vol. 6, pp. 330—336) 
probably does exist at most levels of the troposphere. 
In fact examination of the pertinent sea level and 
700 mb charts indicates that this is so in the layer 
from sea level to 10,000 feet. Our opinion, contrary 
to Mr Phillips’ implication, is that the situation 
described by Namias and Clapp is the frequent 
rather than the rare occurrence in cases of confluence. 


Note on recent study of stability by R. 
Fjortoft 


By VICTOR P. STARR, 
Massachusetts Institute of Technology 
August 1950 


The purpose of this note is to call attention to 
a method for investigating the dynamic stability 
of certain flowsystems recently presented by 
R. FJ&RTORT (1950), and to show that his results 
may be expressed in a clear and concise form 
through the use of the quasi-Lagrangian system 
of hydrodynamical equations discussed by the 
writer (STARR 1945) on a previous occasion. In 
order to achieve clarity, the line of reasoning 
employed by Fjortoft is here repeated in re- 
vised form. 

Let us consider the two-dimensional horizontal 
motion of an inviscid incompressible fluid of 
uniform density on the surface of a sphere of radius 
r, in the region between two fixed latitudes ®,, and 
D, at which walls are assumed to be present so that 
a complete zonal belt extending through all longi- 
tudes A is formed. Following Fjortoft we examine 
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the area integral I of the absolute vorticity & multi- 
plied by sin ® over the entire area A of the belt. 
Either by introducing the proper definition of & in 
terms of the particle velocity and integrating, or by 
the following application of Green’s theorem, I can 
be shown to remain invariant for arbitrary motions 
of the fluid. Since a stream function W must exist 
such that ru = —0W/d®; rv cos D = 0W/0); 
¢ = V2W, where u and v are the eastward and 
northward components of particle velocity and 
V2 is the Laplacian operator for the curved sur- 
face, we have that 


I= ff sin @-V?¥dA =— (j sin a 
— [[r2 cos D aA 


Here ds is an element of the boundary of A and n 
is the outward normal to it. Because the region is 
doubly-connected, the line integral involves the 
use of a cut, but due to the fact that the integrand is 
single-valued the cut may be eliminated and the inte- 
gral separated into two cyclic line integrals along ®,, 
and ®,. Upon substituting the particle velocities for 
the derivatives of Y it may be observed that 1 is 
cqual to a linear combination with constant coefli- 
eients of (a) the circulation around @,, (b) the 
circulation around ®,, and (c) the total angular 
momentum of the fluid (contained in the last 
integral). In view of the fact that the fluid chains 
next to the walls must remain there as the motion 
progresses and that no solenoids or friction are 
present, each of these three quantities and therefore 
also I must remain constant with time f. 

Examining I otherwise, we may first change the 
variables of integration from ®, A to y, A where 

= sin ®. Furthermore we may then in turn change 
these variables from y, A to yy, A where y, is the 
average value of y with respect to longitude for 
a given material chain of particles encircling the 
sphere, according to the quasi-Langrangian scheme 
of hydrodynamical equations discussed elsewhere 
(STARR 1945). Thus 


T=? Jf sin DE cos Ddbd} = r? | | yédydà = 


2 0) 
tr | RA ae dy da. 
0 


If this is applied to motion which represents a 
breakdown from purely zonal flow, we may choose 


dd 


the material chains in such a fashion that &, —- 


dy, 
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are functions of y, alone. In view of the fact 
that under the circumstances assumed d£/dt = 0, 
this functional dependence on 7, is then preserved 
in the subsequent development of the motion. 
Supposing further that we define e by the rela- 
tionship y = y, + & it is simple to show that 
from continuity considerations G edd = o for each 
chain y, = constant. Upon making use of these 
several conditions it follows that 


. Ô 
I=? LG + dE (s +) di = 


9 


5 Be 
= 179 J yobdyydA — me 


where the first integral on the right is the value 
of I corresponding to zonal flow, i.e., when ¢ = 0, 


and y,=y, and / is given by 


A=| a Os etd) | dyo- 
dy : 


Due to the invariant character of J demonstrated 
before, it is clear that A must vanish for all 
motions developing from a state of zonal flow 
and for all motions which develop into zonal 
flow as time progresses. Several cases may be dis- 
tinguished. 

(x) If d£/dy, > o everywhere from y, = yj, to 
Yo = Yq or if dê/dy, <o everywhere in the same 
interval of yg, then /\ cannot vanish unless ¢ = o. 
In these cases zonal flow cannot develop into dis- 
turbed motion. 

(2) If d£/dy, changes sign or becomes zero within 
the interval Yo =Y1 tO Yo = Yo, we may have À =o 
and ¢=40. In this case zonal flow may develop 
into disturbed motion or vice versa, but no in- 
formation is given by the present analysis as to 
whether it will actually do so. 

(3) Preserving the conditions given under (1) let 
us suppose that a disturbance is introduced forcibly 
by some unspecified external agency, after which 
the system is left to its own devices. It follows that 
the finite value of A so introduced cannot of its 
own accord disappear, so that the motion will then 
continue to be disturbed. 

It is interesting to note the following. Under (1) 
stability is implied whatever be the distribution of 
¢ itself so long as the derivative of £ has the properties 
stated. This indicates that criteria which relate 
the presence or absence of stability to the value 


of Ë or to the gradient of angular momentum 
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are not applicable to problems of the present 
type where viscosity is absent and the mo- 
tions are two-dimensional. On the other hand, 
the results recapitulated in this note are in 
agreement with those obtained by Lin (1945) and 
Kuo (1949). 

For certain disturbances covered under (3), ac- 
cording to Lin (1945) and Kuo (1949), the in- 
troduction of a very small viscosity is effective in 
producing damping without at the same time 
leading to a sensible dissipation of the total kinetic 
energy. The introduction of such viscosity cannot, 
however, produce instability for the cases covered 
under (r). 

The invariance of I is in essence a spherical adapta- 
tion of the principle of invariance of the centroid 
of a vortex system as demonstrated, for example, 
in the textbook by Lams (1932). The extremal 
properties of I (as here demonstrated by the character 
of /\) for certain circumstances, are mathematically 
similar to those of the geopotential energy integral 
for properly selected stratified fluids as pointed out 
by Fyorrorr (1950). 
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Floating ice islands in the Arctic Ocean 


BIS Ow Ech RCHER US AE. 
Eilson Air Force Base, Alaska 


November 30th, 1950. 


During his recent visit to Alaska, Dr Sverre 
Petterssen suggested that the readers of Tellus might 
be interested in information on the recent discov- 
eries of large floating ice islands in the Arctic Ocean. 

The undersigned, having reason to suspect the 
existence of such islands, initiated visual and radar 
search in Spring 1950 by the 375th Reconnaissance 
Squadron in connection with the routine weather 
reconnaissance flights to the North Pole (the Ptar- 
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Fig. r. At a distance of 65 miles, T-2 could be 
distinguished clearly from the surrounding pack. This 
photograph was taken at a distance of 45 miles. 


Fig. 2. The high edge of T-2. Note the washboard 
pattern, the drainage system and the steep cliff. 


migan Flight). Although no deviation from the 
standard track could be made in search of such 
islands, the suspicion was soon corroborated. On 
June 3, 1950 the radarscope showed indications of 
a large floating island (called T-1) at 73°15’N and 
159°05’W. This, or a similar ice mass was seen by 
radar ten days later at 74°30’N and 159’05’W. 
Whether or not these radar returns were from 
the same or from two different ice masses could 
not be determined. The relatively large distance 
between them might, however, suggest that they 
were two different islands. 

The first visual identification of floating islands 
was made at 86°40’N and 167°00’E on 20 July 
1950 when three photographs were taken; during 
the following two weeks several additional pictures 
were obtained. This island has been called T-2. 

Island T-3 was discovered toward the end of 
July 1950 when several radar oscilloscope photo- 
graphs were taken of arather small but well-defined 
island at approximately 75°24’N and 173°00’W. 
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On 24 August 1950 this island was identified 
visually and one photograph obtained. 

The islands observed so far have several common 
features which clearly distinguish them from the 
surrounding pack. The surface shows a regular 
corrugated pattern, like that of a washboard, with 
troughs and ridges 300 to 600 feet apart. Superim- 
posed on this pattern is a drainage system of rivers 
about 100 to 200 feet wide, running or meandering 
across the island from one side bounded by a steep 
cliff to the other side where no appreciable cliff is 
apparent, and suggesting a wedge-shaped cross- 
section of the ice mass. Some of these features are 
apparent from the attached photograph. 

The projection of the ice mass above the water 
was variously estimated by the observers from 30 
to 100 feet or more, suggesting a total thickness of 
240 to 800 feet. Experienced drainage engineers 
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who have examined the photographs of T-2 have 
estimated the relief to more than 200 feet, implying 
a minimum thickness of 1600 feet at the high 
edge. 

Although the size of T-1 could not be determin- 
ed, fairly accurate information on the other islands 
is available. The area T-2 was found to be about 
300 square miles; the shape was slightly elliptical 
without any sharp corners, suggesting that the 
island had drifted in the pack for a long period of 
time. On the other hand, T-2, which covered an 
area of about 50 square miles, had a pronounced 
angular shape, which might be taken to indicate a 
relatively short history as a floating island. 

A full report on the information collected dur- 
ing last summer’s exploration is being prepared 
for publication by the American Meteorological 
Society. 
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